好吧算增加了一个模板吧。。
过两天要校赛了,没过的题先放着,这两天把做过的整理整理。
构图思路:将点拆分成2组,一组代表牛,一组代表遮蔽点,然后构建源点和汇点;
用floyd求出点与点之间的最大距离,根据最大距离开始二分,使距离小于mid的点在最大流中开路,若得出最大流等于牛的个数,mid减少,否则增大。求出答案。
这题注意一下构图,先前的sap模板会改变原图的值,且构图的时候注意双向边。
自己的代码WA了。。。
///////////////////////////////
经过好几天的debug 在usaco上找到数据,结果证明我先前的那份模板不能处理大数据。
然后自己用这份模板重新做了一遍,处理floyd上有错 - - wa好几次。
模板:
/*Floyd+最大流判定(拆点)+二分答案*/ #include <iostream> #include <cstdio> #include <algorithm> #include <cstring> #include <queue> #define MAXN 405 #define INF 1e8 #define min(a,b) (a<b?a:b) #define max(a,b) (a>b?a:b) using namespace std; struct edge { int u,v,w,next; }E[200000]; int head[MAXN],ecnt; int gap[MAXN],cur[MAXN],pre[MAXN],dis[MAXN]; __int64 map[205][205]; int N,M,scr,sink,vn,num; int A[205],B[205]; void Insert(int u,int v,int w) { E[ecnt].u=u; E[ecnt].v=v; E[ecnt].w=w; E[ecnt].next=head[u]; head[u]=ecnt++; E[ecnt].u=v; E[ecnt].v=u; E[ecnt].w=0; E[ecnt].next=head[v]; head[v]=ecnt++; } void Init() { int i,u,v; __int64 c; memset(map,-1,sizeof(map)); scr=0;sink=2*N+1;vn=sink+1;num=0; for(i=1;i<=N;i++) { scanf("%d%d",&A[i],&B[i]); num+=A[i]; } for(i=1;i<=M;i++) { scanf("%d%d%I64d",&u,&v,&c); if(map[u][v]==-1||map[u][v]>c) {map[u][v]=map[v][u]=c;} } } void Floyd() { int i,j,k; for(k=1;k<=N;k++) { for(i=1;i<=N;i++) { if(i==k||map[i][k]==-1) continue; for(j=1;j<=N;j++) { if(j==k||i==j||map[j][k]==-1) continue; if(map[i][j]==-1||map[i][j]>map[i][k]+map[k][j]) map[i][j]=map[i][k]+map[k][j]; } } } } int Sap(int s,int t,int n)//核心代码(模版) { int ans=0,aug=INF;//aug表示增广路的流量 int i,v,u=pre[s]=s; for(i=0;i<=n;i++) { cur[i]=head[i]; dis[i]=gap[i]=0; } gap[s]=n; bool flag; while(dis[s]<n) { flag=false; for(int &j=cur[u];j!=-1;j=E[j].next)//一定要定义成int &j,why { v=E[j].v; if(E[j].w>0&&dis[u]==dis[v]+1) { flag=true;//找到容许边 aug=min(aug,E[j].w); pre[v]=u; u=v; if(u==t) { ans+=aug; while(u!=s) { u=pre[u]; E[cur[u]].w-=aug; E[cur[u]^1].w+=aug;//注意 } aug=INF; } break;//找到一条就退出 } } if(flag) continue; int mindis=n; for(i=head[u];i!=-1;i=E[i].next) { v=E[i].v; if(E[i].w>0&&dis[v]<mindis) { mindis=dis[v]; cur[u]=i; } } if((--gap[dis[u]])==0) break; gap[dis[u]=mindis+1]++; u=pre[u]; } return ans; } bool Judge(__int64 mid)//构图+判定 { int i,j; memset(head,-1,sizeof(head));ecnt=0; for(i=1;i<=N;i++) { Insert(scr,i,A[i]); Insert(i+N,sink,B[i]); Insert(i,i+N,INF); } for(i=1;i<=N;i++) { for(j=i+1;j<=N;j++) { if(map[i][j]==-1) continue; if(map[i][j]<=mid)//注意要连入点-->出点(INF) { Insert(i,j+N,INF); Insert(j,i+N,INF); } } } return Sap(scr,sink,vn)==num; } void Solve() { int i,j; __int64 l,r,mid,Max=0; Floyd(); /* for(i=1;i<=N;i++) { for(j=1;j<=N;j++) printf("%I64d ",map[i][j]); printf("\n"); } */ for(i=1;i<=N;i++) { for(j=i+1;j<=N;j++) Max=max(Max,map[i][j]); } l=0;r=Max; while(l!=r) { mid=(l+r)/2; if(Judge(mid)) r=mid; else l=mid+1; } if(Judge(l)) printf("%I64d\n",l); else printf("-1\n"); } int main() { while(scanf("%d%d",&N,&M)!=EOF) { Init(); Solve(); } return 0; }
自己根据这份模板写的代码:
#include<iostream> #include<cstring> #include<cstdio> #define MAXN 405 #define INF 1e8 #define min(a,b) (a<b?a:b) #define max(a,b) (a>b?a:b) using namespace std; struct edge { int u,v,w,next; }E[200000]; __int64 l;__int64 r;__int64 tmp;__int64 mid; int head[MAXN],ecnt; int gap[MAXN],cur[MAXN],pre[MAXN],dis[MAXN]; __int64 map[205][205]; int N,M,scr,sink,vn,num; int A[205],B[205],sum; int k,c1,m,nn; void Floyd() { for(int k=1;k<=N;k++) for(int i=1;i<=N;i++) if(i!=k&&map[i][k]!=-1) for(int j=1;j<=N;j++) { if(i!=j&&j!=k&&map[k][j]!=-1) { if(map[i][j]==-1) map[i][j]=map[i][k]+map[k][j]; else map[i][j]=map[i][j]>map[i][k]+map[k][j]?map[i][k]+map[k][j]:map[i][j]; } } } void Insert(int u,int v,int w) { E[ecnt].u=u; E[ecnt].v=v; E[ecnt].w=w; E[ecnt].next=head[u]; head[u]=ecnt++; E[ecnt].u=v; E[ecnt].v=u; E[ecnt].w=0; E[ecnt].next=head[v]; head[v]=ecnt++; } int Sap(int s,int t,int n)//核心代码(模版) { int ans=0,aug=INF;//aug表示增广路的流量 int i,v,u=pre[s]=s; for(i=0;i<=n;i++) { cur[i]=head[i]; dis[i]=gap[i]=0; } gap[s]=n; bool flag; while(dis[s]<n) { flag=false; for(int &j=cur[u];j!=-1;j=E[j].next)//一定要定义成int &j,why { v=E[j].v; if(E[j].w>0&&dis[u]==dis[v]+1) { flag=true;//找到容许边 aug=min(aug,E[j].w); pre[v]=u; u=v; if(u==t) { ans+=aug; while(u!=s) { u=pre[u]; E[cur[u]].w-=aug; E[cur[u]^1].w+=aug;//注意 } aug=INF; } break;//找到一条就退出 } } if(flag) continue; int mindis=n; for(i=head[u];i!=-1;i=E[i].next) { v=E[i].v; if(E[i].w>0&&dis[v]<mindis) { mindis=dis[v]; cur[u]=i; } } if((--gap[dis[u]])==0) break; gap[dis[u]=mindis+1]++; u=pre[u]; } return ans; } void build(__int64 mid) { memset(head,-1,sizeof(head));ecnt=0; for(int i=1;i<=N;i++)//源点 { Insert(scr,i,A[i]); Insert(i+N,sink,B[i]); Insert(i,i+N,INF); } for(int i=1;i<=N;i++) for(int j=i+1;j<=N;j++) if(map[i][j]!=-1&&map[i][j]<=mid) { Insert(i,j+N,INF); Insert(j,i+N,INF); } } bool judge(__int64 mid) { build(mid); //cout<<l<<' '<<r<<endl; int ans=Sap(scr,sink,vn); //cout<<ans<<"!!!!!!!!!!"<<endl; return ans==sum; } int main() { while(scanf("%d%d",&N,&M)!=EOF) { sum=0;scr=0;sink=2*N+1;vn=sink+1; memset(map,-1,sizeof(map)); for(int i=1;i<=N;i++) { scanf("%d%d",&A[i],&B[i]); sum+=A[i]; } for(int i=1;i<=M;i++) { int v,u; __int64 w; scanf("%d%d%I64d",&v,&u,&w); if(map[v][u]==-1||map[v][u]>w) { map[u][v]=w; map[v][u]=w; } } tmp=0;l=0;r=0;mid=0; Floyd(); for(int i=1;i<=N;i++) //找最短路的最大值 for(int j=i+1;j<=N;j++) { if(map[i][j]!=-1&&map[i][j]>tmp) tmp=map[i][j]; } r=tmp; //cout<<endl<<tmp<<endl; while(l!=r) { mid=(l+r)>>1; if(judge(mid)) { r=mid; } else l=mid+1; } if(judge(l)) printf("%I64d\n",l); else printf("-1\n"); } return 0; }
那份不能处理大数据的挫代码:
#include<cstdio> #include<cstring> #include<iostream> using namespace std; __int64 l,r,mid; const int maxn = 220; const int inf =1000000000;//不要开太大 int c[maxn][maxn]; int q[maxn], pre[maxn]; int level[maxn], gap[maxn]; __int64 map[maxn][maxn]; int A[maxn]; int B[maxn]; int INF =1000000000; int n, m; int s, t; int nn; int cnt,ans; void init_sap(){ memset(level, 1, sizeof (level)); memset(gap, 0, sizeof gap); memset(pre, -1, sizeof pre); int qs = 0, qe = 0; q[qe++] = t; level[t] = 0; gap[ level[t] ] ++; while(qs < qe){ int hd = q[qs++]; for(int i = 1; i <= n; i ++){ if(level[i] > n && c[i][hd] > 0){//level[i] >= n 也可以,why ? q[qe++] = i; level[i] = level[hd] + 1; gap[ level[i] ] ++; } } } } int find_path(int u){ for(int i = 1; i <= n; i ++) if(c[u][i] > 0 && level[u] == level[i] + 1) return i; return -1; } int relabel(int u){ int tmp = inf; for(int i = 1; i <= n; i ++) if(c[u][i] > 0 && tmp > level[i] + 1) tmp = level[i] + 1; if(tmp == inf) tmp = n; return tmp; } int sap(){
init_sap(); int flow = 0, u = s; while(level[s] <= n){ int v = find_path(u); if(v > 0){ pre[v] = u; u = v; if(u == t){ int min_flow = inf; for(int i = t; i != s; i = pre[i]) if(min_flow > c[ pre[i]][i]) min_flow = c[ pre[i]][i]; for(int i = t; i != s; i = pre[i]){ c[pre[i]][i] -= min_flow; c[i][pre[i]] += min_flow; } flow += min_flow; u = s; } }else{ if(-- gap[ level[u]] == 0) return flow; int v = relabel(u); gap[v] ++; level[u] = v; if(u != s) u = pre[u]; } } return flow; } void build() { memset(c, 0, sizeof c); for(int i=1;i<=nn;i++) { c[s][i]=A[i]; c[i+nn][t]=B[i]; c[i][i+nn]=inf; } for(int i=1;i<=nn;i++) //将map中小于mid的路打通c=inf for(int j=1;j<=nn;j++) { if(map[i][j]==-1)continue; if(map[i][j]<=mid) { c[i][j+nn]=inf; c[j][i+nn]=inf; } } /* for(int i=1;i<=n;i++) {cout<<endl; for(int j=1;j<=n;j++) { cout<<c[i][j]<<' '; } }*/ } int main(){ while(~scanf("%d%d", &n, &m)){ //n为点 m为边 nn=n;n=2*nn+2; memset(map,-1,sizeof(map)); for(int i=1;i<=nn;i++) //n的图 { for(int j=1;j<=nn;j++) { if(i==j)map[i][j]=0; } } s = 2*nn+1, t = 2*nn+2;cnt=0; //牛的数量 for(int i=1;i<=nn;i++) //加入源点和汇点的容量 { scanf("%d%d",&A[i],&B[i]); cnt+=A[i]; } for(int i = 1; i <= m; i ++){ int u, v ; __int64 w; scanf("%d%d%I64d", &u, &v, &w);//双向边 if(map[u][v]==-1||w<map[u][v]) map[u][v]=map[v][u]=w; }
for(int k=1;k<=nn;k++) //Floyd for(int i=1;i<=nn;i++) { if(i==k||map[i][k]==-1)continue; for(int j=1;j<=nn;j++) { if(j==k||j==i||map[k][j]==-1)continue; if(map[i][j]==-1||map[i][j]>map[i][k]+map[k][j]) map[i][j]=map[i][k]+map[k][j]; } } l=0,r=0; ans=0; for(int i=1;i<=nn;i++) { for(int j=i+1;j<=nn;j++) r=max(r,map[i][j]); } while(r!=l) //二分 { mid=(l+r)>>1; build(); ans = sap(); if(ans==cnt) { r=mid;
} else l=mid+1; } mid=l; build(); ans=sap(); //r=inf; // printf("%I64d\n",r); //printf("%I64d %I64d %d %d\n",l,r,ans,cnt); if(ans==cnt) printf("%I64d\n",l); else printf("-1\n"); } return 0; }