CJOJ 血帆海盗
Description
随着资本的扩大,藏宝海湾贸易亲王在卡利姆多和东部王 国大陆各建立了N/2 个港口。大灾变发生以后,这些港口之间失去了联系,相继脱离了藏宝海湾贸易亲王的管辖,各自为政。利益的驱动使得每个港口都想和对岸大陆的另一个港口建立 贸易合作关系,由于地理位置因素,只有存在直接到达的航线的两个港口才能建立合作,而且每个港口只与对岸一个港口建立合作,因此并不是所有的港口都能找到 合作伙伴。
血帆海盗得知这一消息以后,决定对其中一条航线进行干扰性的掠夺。经过分析,血帆海盗计算出最多能有W 对港口合作。如果两个港口之间只有一条航线,而且这条航线恰好是血帆海盗要掠夺的航线,这两个港口将不能建立合作关系。血帆海盗指挥官菲尔拉伦想知道他们 有几种选择,可以让地精们无法建立W 对港口。
Input
第1行,两个整数N,M,表示一共的港口个数和航线条数。
接下来M行,每行两个整数A,B,表示卡利姆多的港口A与东部王国的港口B之间有一条航线直接连接,其中1<=A<=N/2,N/2+1<=B<=N。
Output
一个整数,表示血帆海盗可以选择掠夺的航线条数。
解释:
如果掠夺一条航线以后,地精依然可以建立起最多的W个合作关系(可以有多种),那么这条航线是不值得掠夺的,否则就是掠夺方案之一。
Sample Input
8 5
1 5
1 6
2 7
3 7
4 8
Sample Output
1
Hint
样例说明
地精做多能建立起合作关系的数量为3,掠夺(4,8)这条航线后,最多能建立的合作关系的数量减少为2。
数据规模
40%的数据满足2<=N<=200,1<=M<=1000
100%的数据满足2<=N<=100000,1<=M<=100000,保证N为偶数
Source
by BYVoid
网络流 ,连通性,二分图
题目大意,给定一个二分图,问有多少中方案能使断掉一条边后的最大匹配数减小;
首先通过Dinic求出最大流,然后在残量网络上进行tarjan缩强连通分量;
1.首先非匹配边断掉没有卵用;
2.对于一条匹配边,如果他对应的匹配的两个点在一个强连通分量中,说明这条边可以被替换掉,也没啥卵用;
这是为什么呢?
对于残量网络上,向T方向流的边都是表示没有进行匹配的边;
而向S方向流的边都是表示这条边已经匹配了;
所以如果连成了环,也就说明这个环中的匹配边和非匹配可以互相转化
可以画画图加深理解...%%%战舰狗老师
#include<cstdio> #include<algorithm> #include<cmath> #include<iostream> #include<cstring> #include<vector> #define RG register using namespace std; typedef long long ll; const int N=1000000; const int Inf=19260817; int gi() { int x=0; char ch=getchar(); while(ch<'0'||ch>'9') ch=getchar(); while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar(); return x; } int head[N],nxt[N],to[N],s[N],cnt=1,S,T,n,m,sum,q[N],level[N],vis[N],F,c[N]; int dfn[N],low[N],zhan[N],tot,fr[N],tt,ans,vis2[N]; struct data{ int x,y,id; }edge[N]; inline void Addedge(int x,int y,int z) { to[++cnt]=y,s[cnt]=z,nxt[cnt]=head[x],head[x]=cnt; } inline void lnk(int x,int y,int z){ Addedge(x,y,z);Addedge(y,x,0); } inline bool bfs(){ for(RG int i=S;i<=T;i++) level[i]=0,vis[i]=0; int t=0,sum=1; q[0]=S,level[S]=1,vis[S]=1; while(t<sum){ int now=q[t++]; if(now==T) return 1; for(RG int i=head[now];i;i=nxt[i]){ int y=to[i]; if(level[y]==0&&s[i]){ level[y]=level[now]+1; q[sum++]=y; } } } return 0; } inline int dfs(int now,int maxf){ if(now==T) return maxf; int ret=0; for(RG int i=head[now];i;i=nxt[i]) { int y=to[i],f=s[i]; if(level[y]==level[now]+1&&f) { int minn=min(maxf-ret,f); f=dfs(y,minn); s[i]-=f;s[i^1]+=f;ret+=f; if(ret==maxf) break; } } return ret; } inline void Dinic(){ while(bfs()) F+=dfs(S,99999999); } void tarjan(int x){ dfn[x]=low[x]=++tt;zhan[++sum]=x; vis2[x]=1;int y; for(int i=head[x];i;i=nxt[i]){ y=to[i]; if(s[i]){ if(!dfn[y]){ tarjan(y); low[x]=min(low[x],low[y]); } else if(vis2[y]) low[x]=min(low[x],dfn[y]); } } if(low[x]==dfn[x]){ tot++; do { y=zhan[sum--]; vis2[y]=0;fr[y]=tot; } while(y!=x); } } int main(){ n=gi(),m=gi();S=0,T=n+1; for(int i=1;i<=m;i++){ int x=gi(),y=gi(); lnk(x,y,1);edge[i]=(data){x,y,cnt^1}; } for(int i=1;i<=n/2;i++) lnk(S,i,1); for(int i=n/2+1;i<=n;i++) lnk(i,T,1); Dinic();for(int i=S;i<=T;i++) if(!dfn[i]) tarjan(i); for(int i=1;i<=m;i++){ if(s[edge[i].id]==0&&fr[edge[i].x]==fr[edge[i].y]) ans++; else if(s[edge[i].id]==1) ans++; }printf("%d\n",m-ans); return 0; }