多校_BIT_Thieves hdoj3491

这个一看就知道是网络流，但是比赛的时候建图建了好久，我们这方面还有待加强。

就是拆点和最大流。

将每个城市拆成两个点，i*2-1，i*2，gra[i*2-1][i*2]=i 城市的点权，如果城市i和城市j有边，则有边

gra[i*2][j*2-1]和gra[j*2][i*2-1]，权值均为正无穷。

最后将起点和终点所拆出的两点间的边权设为正无穷，然后最大流。

赛后重写顺便复习了下EK和dinic。

 

EK

#include<cstdio> #include<cstring> #include<algorithm> #include<queue> #define inf 10000000 #define N 210 using namespace std; int n,m,sta,end,ans; int gra[N][N],fa[N]; bool s[N]; int EK() { int ans=0,mini,i,u; while (1){ memset(fa,-1,sizeof(fa)); memset(s,0,sizeof(s)); queue <int> que; que.push(sta); s[sta]=1; while (!que.empty()){ u=que.front(); que.pop(); for (i=1;i<=n;i++) if (!s[i]&&gra[u][i]>0){ que.push(i); s[i]=1; fa[i]=u; if (i==end) break; } if (s[end]) break; } if (!s[end]) break; mini=inf; for (i=end;i!=sta;i=fa[i]) mini=min(mini,gra[fa[i]][i]); for (i=end;i!=sta;i=fa[i]){ gra[fa[i]][i]-=mini; gra[i][fa[i]]+=mini; } // printf("%d/n",mini); ans+=mini; } return ans; } int main() { int i,t,a,b; scanf("%d",&t); while (t--){ scanf("%d%d%d%d",&n,&m,&sta,&end); memset(gra,0,sizeof(gra)); for (i=1;i<=n;i++){ scanf("%d",&a); gra[i*2-1][i*2]=a; } for (i=0;i<m;i++){ scanf("%d%d",&a,&b); gra[a*2][b*2-1]=inf; gra[b*2][a*2-1]=inf; } gra[sta*2-1][sta*2]=inf; gra[end*2-1][end*2]=inf; n=2*n; sta=sta*2-1; end=end*2; printf("%d/n",EK()); } return 0; }

 

dinic

#include<cstdio> #include<cstring> #include<algorithm> #define inf 100000000 #define N 210 using namespace std; int n,m,sta,end,ans; int gra[N][N],que[N],d[N]; bool bfs()//搜出层次图 { int tail,head,i,u; tail=head=1; memset(d,-1,sizeof(d)); que[tail++]=sta; d[sta]=0; while (head<tail){ u=que[head++]; for (i=1;i<=n;i++) if (gra[u][i]>0&&d[i]==-1){ d[i]=d[u]+1; que[tail++]=i; } } if (d[end]==-1) return 0; return 1; } int dfs(int k,int sum) { if (k==end) return sum; int os=sum,i,a; for (i=1;i<=n&∑i++) if (d[i]==d[k]+1&&gra[k][i]>0){ a=dfs(i,min(sum,gra[k][i])); gra[k][i]-=a; gra[i][k]+=a; sum-=a; } return os-sum; } int dinic() { int ans=0; while (bfs()) ans+=dfs(sta,inf); return ans; } int main() { int i,t,a,b; scanf("%d",&t); while (t--){ scanf("%d%d%d%d",&n,&m,&sta,&end); memset(gra,0,sizeof(gra)); for (i=1;i<=n;i++){ scanf("%d",&a); gra[i*2-1][i*2]=a; } for (i=0;i<m;i++){ scanf("%d%d",&a,&b); gra[a*2][b*2-1]=inf; gra[b*2][a*2-1]=inf; } gra[sta*2-1][sta*2]=inf; gra[end*2-1][end*2]=inf; n=2*n; sta=sta*2-1; end=end*2; printf("%d/n",dinic()); } return 0; }