2017百度之星资格赛1002:度度熊的王国战略(最小割)
2017-08-06 22:14
525 查看
很大胆的用了最小割 可能因为数据太水 5000+ms过的
也根据这个道题仔细想了想最小割 其实最小割就是把一个连通块分为两部分
思路:
最小割其实就是把一个连通块分成两部分
只不过我们之前做的最小割 比如求从节点1号到节点n号的最小割 割完边后图肯定是分成了两部分 而且1号与n号肯定不在同一部分
这题没有具体要求哪两个点要分开 只要分成两部分就好
所以 只要图被分成了两部分 n号节点肯定在其中一部分 那么肯定存在某个点x在另一部分 这个最小割就是x号到n号的最小割
所以我们只要将源点从1一直赋值到n-1 进行n-1次循环去求最小的最小割就好了
也根据这个道题仔细想了想最小割 其实最小割就是把一个连通块分为两部分
思路:
最小割其实就是把一个连通块分成两部分
只不过我们之前做的最小割 比如求从节点1号到节点n号的最小割 割完边后图肯定是分成了两部分 而且1号与n号肯定不在同一部分
这题没有具体要求哪两个点要分开 只要分成两部分就好
所以 只要图被分成了两部分 n号节点肯定在其中一部分 那么肯定存在某个点x在另一部分 这个最小割就是x号到n号的最小割
所以我们只要将源点从1一直赋值到n-1 进行n-1次循环去求最小的最小割就好了
#include<algorithm> #include<iostream> #include<cstdlib> #include<cstring> #include<cstdio> #include<string> #include<stack> #include<queue> #include<cmath> #include<stack> #include<list> #include<map> #include<set> typedef long long ll; using namespace std; const int MAXN=20010; const int MAXM=880010; const int INF=0x3f3f3f3f; int bian[100005][3];//用来存m组边 int used[3000];//use[i]==0代表i这个节点从未出现过 ==1代表已经出现过了 struct Node { int from,to,next; int cap; }edge[MAXM]; int tol; int head[MAXN]; int dep[MAXN]; int gap[MAXN]; void init() { tol=0; memset(head,-1,sizeof(head)); } void addedge(int u,int v,int w) { edge[tol].from=u; edge[tol].to=v; edge[tol].cap=w; edge[tol].next=head[u]; head[u]=tol++; edge[tol].from=v; edge[tol].to=u; edge[tol].cap=w;//无向图 edge[tol].next=head[v]; head[v]=tol++; } void BFS(int start,int end) { memset(dep,-1,sizeof(dep)); memset(gap,0,sizeof(gap)); gap[0]=1; int que[MAXN]; int front,rear; front=rear=0; dep[end]=0; que[rear++]=end; while(front!=rear) { int u=que[front++]; if(front==MAXN)front=0; for(int i=head[u];i!=-1;i=edge[i].next) { int v=edge[i].to; if(dep[v]!=-1)continue; que[rear++]=v; if(rear==MAXN)rear=0; dep[v]=dep[u]+1; ++gap[dep[v]]; } } } int SAP(int start,int end,int n) //n 为节点的个数 { int res=0; BFS(start,end); int cur[MAXN]; int S[MAXN]; int top=0; memcpy(cur,head,sizeof(head)); int u=start; int i; while(dep[start]<n) { if(u==end) { int temp=INF; int inser; for(i=0;i<top;i++) if(temp>edge[S[i]].cap) { temp=edge[S[i]].cap; inser=i; } for(i=0;i<top;i++) { edge[S[i]].cap-=temp; edge[S[i]^1].cap+=temp; } res+=temp; top=inser; u=edge[S[top]].from; } if(u!=end&&gap[dep[u]-1]==0) break; for(i=cur[u];i!=-1;i=edge[i].next) if(edge[i].cap!=0&&dep[u]==dep[edge[i].to]+1) break; if(i!=-1) { cur[u]=i; S[top++]=i; u=edge[i].to; } else { int min=n; for(i=head[u];i!=-1;i=edge[i].next) { if(edge[i].cap==0)continue; if(min>dep[edge[i].to]) { min=dep[edge[i].to]; cur[u]=i; } } --gap[dep[u]]; dep[u]=min+1; ++gap[dep[u]]; if(u!=start)u=edge[S[--top]].from; } } return res; } void getmap(int m) { init(); int i; for(i=0;i<m;i++) { if(bian[i][0]==bian[i][1]) { continue; } addedge(bian[i][0], bian[i][1], bian[i][2]); } } int main() { int i,j,k; int n,m; while(scanf("%d%d",&n,&m)==2) { memset(used,0,sizeof used);//初始化为0 int cnt=0; for(i=0;i<m;i++) { scanf("%d%d%d",&bian[i][0],&bian[i][1],&bian[i][2]); if(used[bian[i][0]]==0) { used[bian[i][0]]=1; cnt++; } if(used[bian[i][1]]==0) { cnt++; used[bian[i][1]]=1; } } if(cnt!=n)// 如果点并没有全部出现 那么那些没出现的点肯定是孤立的 所建的图就根本不是一个连通块 直接输入0即可 { printf("0\n"); continue; } int source,sink=n; int minn=INF; for(i=1;i<=n-1;i++)//这里非常大胆 直接开了n-1次循环 每次循环都要重新建图 { source=i;//源点每次都赋值为i 汇点值不变 getmap(m);// 重新建图 int tt=SAP(source,sink,n); if(i==1) { minn=tt; } else { minn=min(minn,tt);//更新最小割的值 保存最小的最小割 } } printf("%d\n",minn); } return 0; }
相关文章推荐
- 2017百度之星资格赛 1002 度度熊的王国战略(并查集 or 无向图最小割)
- 2017百度之星资格赛 1002 度度熊的王国战略
- 2017 百度之星资格赛1002 度度熊的王国战略【图+贪心】
- 2017百度之星资格赛:1002. 度度熊的王国战略
- 2017百度之星资格赛 度度熊的王国战略
- 2017"百度之星"程序设计大赛-资格赛-1002-度度熊的王国战略
- 2017百度之星资格赛1002度度熊的王国战略
- hdu 6081 度度熊的王国战略 2017"百度之星"程序设计大赛 - 资格赛 1002
- 2017"百度之星"程序设计大赛 - 资格赛:1002 度度熊的王国战略
- 2017"百度之星"程序设计大赛 - 资格赛 度度熊的王国战略
- 2017Astar资格赛1002 度度熊的王国战略
- 2017百度之星资格赛 1002-图 1003-DP 1004-01背包(解选择方案)
- 【2017"百度之星"程序设计大赛 - 资格赛-1002度度熊的王国战略-并差集】
- hdu6081--度度熊的王国战略(百度之星资格赛)
- 2017百度之星资格赛:1001. 度度熊保护村庄(floyd最小环)
- 2017百度之星资格赛(A)1005 今夕何夕(日期计算)
- 2017百度之星资格赛1003 度度熊与邪恶大魔王(完全背包)
- 2017百度之星-1001-度度熊保护村庄(floyd求最小环)
- 2017百度之星初赛:A-1002. 数据分割(并查集+set)
- 【2017"百度之星"程序设计大赛 - 资格赛】 度度熊与邪恶大魔王