太空飞行计划问题(最大流)

题目

---

最大权闭合图。

对于每个实验，连一条从s到实验，边权为实验赞助的边。

对于每个需要的仪器，连一条从实验到器材，边权为INF的边。

对于每个仪器，连一条从器材到t，边权为器材费用的边。

求最小割即最大流，然后用实验赞助总和减去最大流即为最大净收益

如图，我们求出最大流后，再拿总利润减去时=>可得将原利益减掉$<=>$不用这个实验

至于求路径，个人觉得数据太水，故dinic最后一次bfs的d数组正好可以用来判断这个
但可能只有暴力（由于n,m只有50，枚举删除器材和t之间的边，如果删去后的最大流和原来的最大流的差值等于这条边的边权，那么这条边就是必须满流的，也就是这个器材是必要的。再根据需要的器材我们容易知道要做的实验有哪些。）比较符合心意

#include<bits/stdc++.h>
#define re return
#define inc(i,l,r) for(int i=l;i<=r;++i)
const int maxn=105,maxm=100005,inf=1e9;
using namespace std;
int flag;
template<typename T>inline void rd(T&x)
{
    char c;bool f=0;
    while((c=getchar())<'0'||c>'9')if(c=='-')f=1;
    x=c^48;
    while((c=getchar())>='0'&&c<='9')x=x*10+(c^48);
    if(f)x=-x;
    if(c=='\n')flag=1;
}

int s,t,m,n,k=1,hd[maxn],deep[maxn],cur[maxn],use[maxn],vis[maxn];
struct node{
    int to,nt,flow;
}e[maxm];

inline void add(int x,int y,int z)
{
    e[++k].to=y;e[k].nt=hd[x];hd[x]=k;e[k].flow=z;
    e[++k].to=x;e[k].nt=hd[y];hd[y]=k;e[k].flow=0;
}

inline bool bfs()
{
    inc(i,1,t)deep[i]=0;
    deep[s]=1;
    queue<int>q;
    q.push(s);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int i=hd[u];i;i=e[i].nt)
        {
            int v=e[i].to;
            if(!deep[v]&&e[i].flow)
            {
                deep[v]=deep[u]+1;
                q.push(v);
            }
        }
    }
    re deep[t];
}
inline int dfs(int u,int flow)
{
    if(u==t) re flow;
    int delta=flow;
    for(int &i=cur[u];i;i=e[i].nt)
    {
        int v=e[i].to;
        if(deep[v]==deep[u]+1&&e[i].flow)
        {
            int d=dfs(v,min(delta,e[i].flow));
            e[i].flow-=d;
            e[i^1].flow+=d;
            delta-=d;
            if(!delta)break;
        }
    }
    re flow-delta;
}
void xfs(int u)
{
    for(int i=hd[u];i;i=e[i].nt)
    {
        int v=e[i].to;
        vis[v]=1;
    }
}
int main()
{
    //freopen("in.txt","r",stdin);
    int x,totans=0;
    rd(n),rd(m);
    s=n+m+2;
    t=s+1;
    inc(i,1,n)
    {
        rd(x);add(s,i+m,x);
        flag=0;
        totans+=x;
        while(flag!=1)
        {
            rd(x);
            add(i+m,x,inf);
        }
    }
    
    inc(i,1,m)
    {
        rd(x);
        add(i,t,x);
    }
    
    int ans=0; 
    while(bfs()){
        inc(i,1,t)cur[i]=hd[i];
        ans+=dfs(s,inf);
    } 
    
    inc(i,1,m)vis[i]=0;
    inc(i,1,n)
    if(deep[i+m])printf("%d ",i),xfs(i+m);
    printf("\n");
    
    inc(i,1,m)
    if(vis[i])printf("%d ",i);
    printf("\n");
    printf("%d",totans-ans);
    re 0;
}