[补档][COGS 426]血帆海盗

[COGS 426]血帆海盗

题目

传送门：http://cogs.pro/cogs/problem/problem.php?pid=426

随着资本的扩大，藏宝海湾贸易亲王在卡利姆多和东部王国大陆各建立了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

OUTPUT

1

解题报告

好难啊= =，在编译器炸了的情况下生交了3遍，然后各种编译错误加不过样例= =

最小割的唯一性判定

在残余网络上跑tarjan求出所有SCC，记id[u]为点u所在SCC的编号。显然有id[s]!=id[t]（否则s到t有通路，能继续增广）。

①对于任意一条满流边(u,v)，(u,v)能够出现在某个最小割集中，当且仅当id[u]!=id[v]；

②对于任意一条满流边(u,v)，(u,v)必定出现在最小割集中，当且仅当id[u]==id[s]且id[v]==id[t]。

①<==将每个SCC缩成一个点，得到的新图就只含有满流边了。那么新图的任一s-t割都对应原图的某个最小割，从中任取一个把id[u]和id[v]割开的割即可证明。

②<==：假设将(u,v)的边权增大，那么残余网络中会出现s->u->v->t的通路，从而能继续增广，于是最大流流量（也就是最小割容量）会增大。这即说明(u,v)是最小割集中必须出现的边。

然后就可以跑了= =

  1 #include<iostream>
  2 #include<cstring>
  3 #include<cstdio>
  4 #include<queue>
  5 using namespace std;
  6 inline int read(){
  7     int sum(0);
  8     char ch(getchar());
  9     for(;ch<'0'||ch>'9';ch=getchar());
 10     for(;ch>='0'&&ch<='9';sum=sum*10+(ch^48),ch=getchar());
 11     return sum;
 12 }
 13 struct edge{
 14     int e,n,w;
 15 }a[400005];
 16 int pre[100005],tot;
 17 inline void insert(int s,int e,int w){
 18     a[tot].e=e;
 19     a[tot].w=w;
 20     a[tot].n=pre[s];
 21     pre[s]=tot++;
 22 }
 23 int n,m;
 24 int ed;
 25 int S(0),T;
 26 int ans(0),inf(0x7fffffff);
 27 int dis[100005];
 28 inline bool bfs(int s,int t){
 29     memset(dis,0,sizeof(dis));
 30     dis[s]=1;
 31     queue<int>q;
 32     q.push(s);
 33     while(!q.empty()){
 34         int k(q.front());
 35         q.pop();
 36         for(int i=pre[k];i!=-1;i=a[i].n){
 37             int e(a[i].e);
 38             if(!dis[e]&&a[i].w){
 39                 dis[e]=dis[k]+1;
 40                 q.push(e);
 41                 if(e==t)
 42                     return true;
 43             }
 44         }
 45     }
 46     return false;
 47 }
 48 inline int my_min(int a,int b){
 49     return a<b?a:b;
 50 }
 51 inline int dfs(int now,int flow){
 52     if(now==T)
 53         return flow;
 54     int tmp(flow),f;
 55     for(int i=pre[now];i!=-1;i=a[i].n){
 56         int e(a[i].e);
 57         if(dis[e]==dis[now]+1&&tmp&&a[i].w){
 58             f=dfs(e,my_min(tmp,a[i].w));
 59             if(!f){
 60                 dis[e]=0;
 61                 continue;
 62             }
 63             a[i].w-=f;
 64             a[i^1].w+=f;
 65             tmp-=f;
 66         }
 67     }
 68     return flow-tmp;
 69 }
 70 inline void dinic(int s,int t){
 71     while(bfs(s,t))
 72         ans+=dfs(s,inf);
 73 }
 74 int dfn[100005],low[100005],bl[100005],stack[100005];
 75 int qlt,cnt,top;
 76 bool vis[100005];
 77 inline void tarjan(int u){
 78     dfn[u]=low[u]=++cnt;
 79     stack[++top]=u;
 80     vis[u]=1;
 81     for(int i=pre[u];i!=-1;i=a[i].n)
 82         if(a[i].w){
 83             int e(a[i].e);
 84             if(!dfn[e]){
 85                 tarjan(e);
 86                 low[u]=my_min(low[u],low[e]);
 87             }
 88             else
 89                 if(vis[e])
 90                     low[u]=my_min(low[u],dfn[e]);
 91         }
 92     if(low[u]==dfn[u]){
 93         qlt++;
 94         while(1){
 95             int tmp(stack[top--]);
 96             vis[tmp]=0;
 97             bl[tmp]=qlt;
 98             if(tmp==u)
 99                 break;
100         }
101     }
102 }
103 inline int gg(){
104     freopen("bloodsail.in","r",stdin);
105     freopen("bloodsail.out","w",stdout);
106     memset(pre,-1,sizeof(pre));
107     n=read(),m=read();
108     ed=n>>1;
109     T=n+1;
110     for(int i=1;i<=m;i++){
111         int x(read()),y(read());
112         insert(x,y,1),insert(y,x,0);
113     }
114     for(int i=1;i<=ed;i++)
115         insert(S,i,1),insert(i,S,0),insert(i+ed,T,1),insert(T,i+ed,0);
116     dinic(S,T);//cout<<ans<<endl;
117     for(int i=0;i<=n+1;i++)
118         if(!dfn[i])
119             tarjan(i);
120     for(int i=1;i<=ed;i++)
121         for(int j=pre[i];j!=-1;j=a[j].n)
122             if(!a[j].w&&bl[i]==bl[a[j].e]&&a[j].e)
123                 ans--;
124     printf("%d",ans);
125     return 0;
126 }
127 int K(gg());
128 int main(){;}

View Code