[Luogu P1345] [USACO5.4]奶牛的电信Telecowmunication (最小割)

题面

传送门:https://www.luogu.org/problemnew/show/P1345

[Luogu P1345] [USACO5.4]奶牛的电信Telecowmunication (最小割)]

[Luogu P1345] [USACO5.4]奶牛的电信Telecowmunication (最小割)


Solution

这道题,需要一个小技巧了解决。

我相信很多像我这样接蒟蒻,看到这道题,不禁兴奋起来:“这道题是裸的割边,我会做!!!”

然后兴冲冲的打了个DINIC,交一发,80分。

所以说我们有时候还是太naive。

重新读题,会发现这题割的不是边,是点。这样还能80分,数据真水

所以说,我们需要一个割边转割点的小技巧。

我们可以考虑“拆点”,即把一个点拆成两个点,中间连一条边权为1的边

前一个点作为“入点”,别的点连边连入这里

后一个点作为“出点”,出去的边从这里出去

这样,只要我们切断中间那条边,就可以等效于除去这个点

如图:

[Luogu P1345] [USACO5.4]奶牛的电信Telecowmunication (最小割)

红色的边边权为1,黑色的边边权为inf。

原点和汇点的内部边权为inf,因为显然这两个点不能删除。

题面给的边删除没意义(因为我们要删点),所以也设为inf(事实上设为1也没问题,因为删除这条边的权值可以理解为删除了一个点)

至于怎么算割边,可以证明割边在数值上等于最大流(本蒟蒻不会证)

至于怎么求最大流........可以参考这个博客:http://www.cnblogs.com/SYCstudio/p/7260613.html

最后记得双倍空间

然后就OjbK了


Code

//Luogu P1345 [USACO5.4]奶牛的电信Telecowmunication
//June,3rd,2018
//割边转割点
#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
long long read()
{
long long x=0,f=1; char c=getchar();
while(!isdigit(c)){if(c=='-') f=-1; c=getchar();}
while(isdigit(c)){x=x*10+c-'0';c=getchar();}
return x*f;
}
const int N=200+10;
const int inf=0x3f3f3f3f;
struct road
{
int to,w,rev;
road (int A,int B,int C)
{
to=A,w=B,rev=C;
}
};
vector <road> e[N];
int n,m,c1,c2,depth[N];
queue <int> dl;
bool bfs()
{
memset(depth,0,sizeof depth);
depth[c1]=1;
dl.push(c1);
while(dl.empty()==false)
{
int now=dl.front();
dl.pop();
for(int i=0;i<int(e[now].size());i++)
if(e[now][i].w>0 and depth[e[now][i].to]==0)
{
depth[e[now][i].to]=depth[now]+1;
dl.push(e[now][i].to);
}
}
if(depth[c2]==0) return false;
return true;
}
int dfs(int now,int f)
{
if(now==c2) return f;
int ans=0;
for(int i=0;i<int(e[now].size());i++)
if(e[now][i].w>0 and depth[e[now][i].to]==depth[now]+1)
{
int temp=dfs(e[now][i].to,min(f,e[now][i].w));
e[now][i].w-=temp;
e[e[now][i].to][e[now][i].rev].w+=temp;
f-=temp,ans+=temp;
if(f==0) break;
}
return ans;
}
int Dinic()
{
int ans=0;
while(bfs()==true)
ans+=dfs(c1,inf);
return ans;
}
inline void AddLine(int s,int t,int w)
{
e[s].push_back(road(t,w,e[t].size()));
e[t].push_back(road(s,0,e[s].size()-1));
}
int main()
{
n=read(),m=read(),c1=read(),c2=read();
for(int i=1;i<=n;i++) e[i].reserve(8);
for(int i=1;i<=n;i++)
if(i==c1 or i==c2)
AddLine(i,i+n,inf);
else
AddLine(i,i+n,1);
for(int i=1;i<=m;i++)
{
int s=read(),t=read();
AddLine(s+n,t,inf);
AddLine(t+n,s,inf);
} printf("%d",Dinic());
return 0;
}

正解(C++)

上一篇:LK8382-删除排序链表的重复元素(保留一个重复节点+删除全部重复节点)


下一篇:洛谷——P1345 [USACO5.4]奶牛的电信Telecowmunication