woc这题目的输入格式和输出格式真的恶心
首先我们就着样例讲一下闭合图
如图所示，第一层是两个实验节点，带来正收益；第二层是三个仪器节点，带来负收益；问讲道理到终点可以获得多大收益。
闭合图是什么呢？闭合图是一个点集，这个点集中所有点的出边所指向的点都必须在闭合图中。我们用点权来表示点的话，比如点集{10，-5，-6，终点}就是个闭合图，然而{10,25}就不是。
然后这题显而易见的叫我们求出最大权闭合图。
于是我们可以把它变成这样一个图：
如图，从源点到每个正点权点连一条容量为点权的边，从正点权点向负点权点连一条容量为INF的边，再从负点权点向汇点连一条容量是该点点权绝对值的边。
这样原来的求最大权闭合图变成求这个图的最小割。证明在这边
这样就完啦
#include<cstdio> #include<cstring> #include<cctype> #include<cstdlib> #include<algorithm> #include<iostream> #include<queue> #define maxn 1000 #define maxm 1000 using namespace std; inline long long read(){ long long num=0,f=1; char ch=getchar(); while(!isdigit(ch)){ if(ch==‘-‘) f=-1; ch=getchar(); } while(isdigit(ch)){ num=num*10+ch-‘0‘; ch=getchar(); } return num*f; } inline int count(int i){ return i&1?i+1:i-1; } struct Edge{ int next,to,val; }edge[maxm*4]; int head[maxn*3],num; inline void addedge(int from,int to,int val){ edge[++num]=(Edge){head[from],to,val}; head[from]=num; } inline void add(int from,int to,int val){ addedge(from,to,val); addedge(to,from,0); } bool vis[maxn*3]; int dfn[maxn]; int dis[maxn]; bool ext[maxn]; int list[maxn]; int Start,End; bool bfs(){ memset(vis,0,sizeof(vis)); memset(dis,0,sizeof(dis)); queue<int>f; f.push(Start);vis[Start]=1;dfn[Start]=1; while(!f.empty()){ int from=f.front();f.pop(); for(int i=head[from];i;i=edge[i].next){ int to=edge[i].to; if(edge[i].val<=0||vis[to]) continue; dfn[to]=dfn[from]+1; vis[to]=1; dis[to]=dis[from]+1; f.push(to); } } return vis[End]; } int dfs(int x,int val){ //printf("%d\n",x); if(x==End||val==0) return val; int flow=0;vis[x]=1; for(int &i=list[x];i;i=edge[i].next){ int to=edge[i].to; if(edge[i].val<=0||dfn[to]!=dfn[x]+1||vis[to]) continue; int now=dfs(to,min(val,edge[i].val)); edge[i].val-=now;edge[count(i)].val+=now;flow+=now;val-=now; if(val<=0) break; } if(val!=flow) dfn[x]=-1; return flow; } int maxflow(){ int ans=0; while(bfs()){ memset(vis,0,sizeof(vis)); for(int i=0;i<=End;++i) list[i]=head[i]; int now=dfs(Start,0x7fffffff); if(!now) break; ans+=now; } return ans; } int cnt; bool mst[maxn+maxm+10]; int c[maxn][maxm]; int tot[maxm]; int main(){ int m=read(),n=read();End=n+m+1; string S; for(int i=1;i<=m;++i){ int x=read(); cnt+=x; add(Start,i,x); getline(cin,S); int len=S.length(); for(int j=0;j<len;++j){ if(S[j]==‘ ‘) continue; x=0; while(isdigit(S[j])){ x=x*10+S[j]-‘0‘; j++; } c[i][++tot[i]]=x; } for(int j=1;j<=tot[i];++j) add(i,c[i][j]+m,0x7fffffff); } for(int i=1;i<=n;++i){ int x=read(); add(i+m,End,x); } int ans=maxflow(); for(int i=1;i<=m;++i) if(dis[i]^0) printf("%d ",i); printf("\n"); for(int i=1;i<=n;++i) if(dis[i+m]^0) printf("%d ",i); printf("\n%d",cnt-ans); return 0; }
