poj 1149 SAP 模板使用

构图想法：

 

 

假设某一个顾客能打开一定的猪圈，那么这个顾客将把所有全买了，然后，下一次，如果有人能打开和上一个顾客相同的猪圈，那么从上一个顾客到下一个顾客就有一个无穷大的边到这个顾客，这个顾客也将买走所有他能买的。

 

从顾客到汇点，容量就是，此顾客能买的猪的数目。

 

 

 

#include<stdio.h> #include<memory.h> #include<iostream> using namespace std; const int MAXNUM = 2100000000 ; //权值上限 const int MAXPT = 400 ; //顶点数上限 const int MAXEG = MAXPT*MAXPT*2 ; //边数上限 const int MAXQUE = 200000 ; /* s = 1 ; // 源点 t = n ; // 汇点 */ class MNF_SAP { private: //int m,n; int s,t; int d[MAXPT]; //距离标号 int pre[MAXPT]; //前置顶点 int flow[MAXPT]; //到当前点为止，所有弧的最小值 int cur[MAXPT]; //当前弧 int cnt[MAXPT]; //k标号出现次数 int que[MAXPT*2],ft,rt; bool vis[MAXPT]; void init_d() { int i,u; memset(vis,0,sizeof(vis)); ft=rt=0; d[t]=0; vis[t]=true; que[++rt]=t; while(ft!=rt) { u=que[(++ft)%MAXQUE]; for(i=gra[u]; i; i=E[i].next) if(!vis[E[i].v] && !E[i].w) { d[E[i].v]=d[u]+1; vis[E[i].v]=true; que[(++rt)%MAXQUE]=E[i].v; } } } public: struct { int v,w,next; } E[MAXEG]; int n; int k; int gra[MAXPT]; MNF_SAP(){} ~MNF_SAP(){} void adde(int x,int y,int z) { E[++k].v = y ; E[k].next = gra[x] ; E[k].w = z ; gra[x] = k ; E[++k].v = x ; E[k].next = gra[y] ; E[k].w = 0 ; gra[y] = k ; } int SAP() { int u,v,i,maxflow=0; //总最大流 s=1; t=n; u=s; flow[s]=MAXNUM; for(i=1; i<=n; i++) cur[i]=gra[i]; //当前弧初始化 init_d(); for(i=1; i<=n; i++) cnt[d[i]]++; cnt[0]=n; while(d[s]<n) { for(i=cur[u]; i; i=E[i].next) //找允许弧 if(E[i].w && d[E[i].v]+1==d[u]) break; if(i) //存在允许弧 { cur[u]=i; //设置当前弧 v=E[i].v; if(E[i].w<flow[u]) flow[v]=E[i].w; else flow[v]=flow[u]; //标记当前顶点为止经过的最小弧 u=v; pre[v]=i; //前置顶点边号 if(u==t) { do { E[pre[u]].w-=flow[t]; //正向弧减a[t] E[pre[u]^1].w+=flow[t]; //通过异或操作找反向弧 u=E[pre[u]^1].v; } while(u!=s); maxflow+=flow[t]; //memset(flow,0,sizeof(flow)); flow[s]=MAXNUM; } } else //不存在允许弧 { if(--cnt[d[u]]==0) break; //间隙优化 d[u]=n; cur[u]=gra[u]; for(i=gra[u]; i; i=E[i].next) if(E[i].w && d[E[i].v]+1<d[u]) d[u]=d[E[i].v]+1; //修改距离标号为 最小的非允许弧加1 cnt[d[u]]++; if(u!=s) u=E[pre[u]^1].v; } } return maxflow; } } ; int pigs[1005] ; int pbs[105][1005]; int bs(int a[],int f,int t,int x) // f is from 1 { int m ; while(f<=t) { m = (f+t)>>1 ; if(a[m]==x) return m; else if( a[m] < x ) { f = m + 1 ; } else { t = m - 1 ; } } return 0; } int main() { MNF_SAP obj ; int m,n; int i,j,k; int a,b; int tem; int sum ; while( cin>>m>>n ) { //init obj memset(obj.gra,0,sizeof(obj.gra)); obj.k = 1 ; obj.n = n+2; // input for(i=1;i<=m;i++) scanf("%d",&pigs[i]); for(i=1+1;i<=n+1;i++) { scanf("%d",&pbs[i][0]); for(j=1;j<=pbs[i][0];j++) { scanf("%d",&pbs[i][j]); for(k=1+1;k<i;k++) { if( bs(pbs[k],1,pbs[k][0],pbs[i][j]) ) { obj.adde(k,i,MAXNUM); break; } } } for(sum=0,j=1;j<=pbs[i][0];j++) { sum += pigs[pbs[i][j]] ; pigs[pbs[i][j]] = 0 ; } obj.adde(1,i,sum); scanf("%d",&tem); obj.adde(i,n+2,tem); } cout<<obj.SAP()<<endl; // working } return 0; }