洛谷P2472 [SCOI2007]蜥蜴 题解

题目链接:

https://www.luogu.org/problemnew/show/P2472

分析:

这道题用最大流解决。

首先构建模型。

一根柱子可以跳入和跳出,于是拆成两个点:入点和出点。

每一根柱子的入点和出点连一条流量为高度的边,来限制蜥蜴跳入的次数。

当柱子a可以调到柱子b时,就从a的出点向b的入点连边,流量inf。

S向所有有蜥蜴的柱子的入点连边,流量为1

T表示地图外一点,当一根柱子能跳到地图外时,则出点向T连流量为inf的边。

然后跑最大流即可。

这里要注意数组的范围以及拆点。

代码:

#include<cstdio>
#include<cstring>
#include<queue>
#include<vector>
#include<iostream> 
#define inf 0x7fffffff
using namespace std;
int s,t,ans;
int d[1005];
struct edge
{
	int to,val,rev;
	edge(int _to,int _val,int _rev)
	{
		to=_to;
		val=_val;
		rev=_rev;
	}
};
vector<edge>e[1005]; 
void add(int x,int y,int w)
{
	e[x].push_back(edge(y,w,e[y].size()));
	e[y].push_back(edge(x,0,e[x].size()-1)); 
}
bool bfs()
{
    memset(d, -1, sizeof(d));
    queue<int> q;
    q.push(s);
    d[s]=0;
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
        for(int i=0;i<e[x].size();i++)
        {
            int y=e[x][i].to;
            if(d[y]==-1 && e[x][i].val)
            {
                q.push(y);
                d[y]=d[x]+1;
            }
        }
    }
    if(d[t]==-1)
        return 0;
    else
        return 1;
}

int dfs(int x,int low)
{
    if(x==t || low==0)
        return low;
    int totflow=0; 
    for(int i=0;i<e[x].size();i++)
    {
        int y=e[x][i].to;
        int rev=e[x][i].rev;
        if(d[y]==d[x]+1 && e[x][i].val) 
        {
            int a=dfs(y,min(low,e[x][i].val));
            e[x][i].val-=a;
            e[y][rev].val+=a;
            low-=a;
            totflow+=a;
            if(low==0)
                return totflow;
        }
    }
    if(low!=0)
        d[x]=-1;
    return totflow;
}

void dinic()
{
	while(bfs())
	{
		ans+=dfs(s,inf);
	}
}
int main()
{
	int n,m,c,cnt=0;
	char ss;
	scanf("%d%d%d",&n,&m,&c);
	s=0,t=n*m*2+1;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			cin>>ss;
			if(ss>'0')
			{
				add((i-1)*m+j,(i-1)*m+j+n*m,ss-'0');
			} 
		}
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			cin>>ss;
			if(ss=='L')
			{
				add(s,(i-1)*m+j,1);
				cnt++;
			}
		}
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			if(i+c>n||i-c<1||j+c>m||j-c<1)
			{
				add((i-1)*m+j+n*m,t,inf);
			}
		}
	} 
	for(int x1=1;x1<=n;x1++)
	{
		for(int y1=1;y1<=m;y1++)
		{
			for(int x2=1;x2<=n;x2++)
			{
				for(int y2=1;y2<=m;y2++)
				{
					if((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)<=c*c)
					{
						add((x1-1)*m+y1+n*m,(x2-1)*m+y2,inf);
					}
				}
			}
		}
	}
	dinic(); 
	printf("%d",cnt-ans); 
	return 0;
}
全部评论

相关推荐

不愿透露姓名的神秘牛友
11-21 17:16
科大讯飞 算法工程师 28.0k*14.0, 百分之三十是绩效,惯例只发0.9
点赞 评论 收藏
分享
一颗宏心:华为HR晚上过了十二点后还给我法消息。
点赞 评论 收藏
分享
点赞 收藏 评论
分享
牛客网
牛客企业服务