一些代码最大权闭合图

HDU49712014多校赛一共有n个任务，m个技术，完成一个任务可一些钱，学个#define_CRT_SECURE_NO_WARNINGS#include<cstdio>#include<cstring>#include<vector>#include<algorithm>#definesz1000#defineinf0x3fffffffusingnamespacestd;intn,m;inta[sz],b[sz];vector<int>p[sz];vector<int>g[sz],intlow[sz],dfn[sz],stack[sz],instack[sz],belong[sz],cost[sz];intdeep,top,kinds;voidtarjan(ints)//tarjan{low[s]=dfn[s]=deep++;stack[++top]=s;instack[s]=1;for(inti=0;i<g[s].size();{intt=if(dfn[t]==-{tarjanlow[s]=min(low[t],}elseif{low[s]=min(dfn[t],}}if(low[s]=={intintct=0;{t=stack[top--belong[t]=kinds;//每个属于哪个缩点后的点instack[t]=0;ct+=b[t];}while(t!=cost[kinds]ct;//cost记录缩点后的点集中包含的几个技术的花费和}}structnode{ints,t,v,nxt;}e[sz*inthd[sz],dis[sz],cur[sz],que[sz],cnt,S,T;voidinsert(ints,intt,intv){e[cnt].s=s,e[cnt].t=t,e[cnt].v=v,e[cnt].nxt=hd[s],hd[s]=e[cnt].s=t,e[cnt].t=s,e[cnt].v=0,e[cnt].nxt=hd[t],hd[t]=}boolbfs(ints,intt,intn){inthead=0,tail=0,i;memset(dis,-1,sizeof(dis[0])*(n+1)),dis[s]=0;que[tail++]=s;whilefor(i=hd[que[head++]];i!=-1;i=if(e[i].v>0&&dis[e[i].t]==-dis[e[i].t]=dis[e[i].s]+1;if(e[i].t==t)return1;que[tail++]=e[i].t;}}return}intdinic(ints,intt,intn){//最大流模板,这个无视，暂时还只会套模板，有点看不懂inti,mf=0,tp;while(bfs(s,t,for(i=0;i<n;i++)cur[i]=hd[i];intu=s,tail=0;while(cur[s]!=-if(u!=t&&cur[u]!=-1&&e[cur[u]].v>0&&dis[u]+1==que[tail++]=cur[u];u=}elseif(u==for(tp=inf,i=tail-1;i>=0;i--)tp=min(tp,e[que[i]].v);for(mf+=tp,i=tail-1;i>=0;i--){e[que[i]].v-=tp;e[que[i]^1].v+=tp;if(e[que[i]].v==0)tail=i;}u=}elsewhile(u!=s&&cur[u]==-1)u=e[que[--tail]].s;cur[u]=e[cur[u]].nxt;}}}return}voidmake_map{for(inti=0;i<kinds+n;i++)ng[i].clear();for(inti=0;i<m;for(intj=0;j<g[i].size();{intt=if(belong[i]!={ng[belong[i]].push_back}}for(inti0;i<n;i++)//下标[0,kinds)是技术的缩点,[kinds,kinds+n){intid=i+for(intj=0;j<p[i].size();ng[id].push_back(belong[p[i][j]]);//项目和技术之间建边，belong[p[i][j]]技}for(inti=0;i<kinds+n;{sort(ng[i].begin(),ng[i].endng[i].resize(unique(ng[i].beginng[i].endng[i].begin一个项}memset(hd,-1,sizeof(hd));cnt=0;S=kinds+T=S+for(inti0;i<n;i++)//源点s往所有的项目建流量为收益的边insert(S,i+kinds,a[i]);for(inti=0;i<n+kinds;for(intj0;j<ng[i].sizej++)//中间的点的流量全是infinsert(i,ng[i][j],inf);for(inti0;i<kinds;i++)//所有技术往汇点tinsert(i,T,cost[i]);}intmain{intscanf("%d",for(intcas=1;cas<=cs;{scanf("%d%d",&n,&m);intsum=0;for(inti=0;i<n;scanf("%d",&a[i]),sum+=a[i];for(inti=0;i<m;i++)scanf("%d",&b[i]);intnum,c;for(inti=0;i<n;{scanf("%d",p[i].clearfor(intj=0;j<num;{scanf("%d",}}for(inti=0;i<m;{g[i].clearfor(intj=0;j<m;{scanf("%d",&c);if(c)}}memset(dfn,-1,sizeof(dfn));top=deep=kinds=0;for(inti=0;i<m;i++)if(