P2762 太空飞行计划问题

本来是想做 NOI2009的植物大战僵尸来着，建了半天模型都没又成功。
于是看了看题解，发现这个问题叫做最大权闭合图
我TM好像没学过啊.
于是先看了看概念，做了下入门的题目，然后一做就是一晚上（惨~
我真不是颓了来着。

概念

这个题引入了最大权闭合图的概念。一个闭合图的定义为图中所有点的出边的终点都在这个图中。一个图的最大权闭合子图为其中所有闭合子图中权值和最大的。
本题的实质就是求最大权闭合图，试按照实验-仪器的关系建图，并将它们的收益（如果要花钱，则为负值）标注在点上。这个图应当形如以下示例。

要求本图上的最大权闭合图，就要先将原图中的正值点连向源点S，容量为点权，负值点连向汇点T，容量为点权的绝对值，例如下图示例。

需要用到的定理：
1.最大流==最小割
2.最小割分割成的两个图中，源点S所在的那个图是最大权闭合图。

建模：

在本题里，我们需要从源点向每个实验建一条容量为价值的边
从实验向所需器材建一条容量为INF的边
从器材向汇点连一条容量为花费的边
然后跑最大流。
答案即为所以实验的总价值-最大流。
为什么？？？
我们可以想象，要实现最小割，就要割掉最小的实验边和最小的器材边。
这些实验是没有实现的，而这些器材是用到的。
那么总的收入就等于总代价-没有实现的代价-器材带来的成本

输出方案：

可以看最后一次查出阻塞流的时候有哪些点被BFS遍历到了，如果被遍历到了说明跟源点S在一个子图内

坑点：

输入没有0作为最后一个数结束
那么我们就特判\n为结束符不就行了？？？
没错，在本机上是可以的
但是在OJ上死活就是不A，当我下载了第一个点的时候，发现这™就是样例
显示我输出的答案是12？？？（我怎么可能不过样例就交啊
然而我没有意识到这个问题，继续找错，一拖就是一晚上(大概1.30小时
最后，我发现，为啥人家特判的换行符是’\r’？？？
难道不是‘\n’吗？
我怀着试试的心态一交，A了？！
好吧，这出题人绝对是故意的
本傻B补充一下，其实是在Linux下进行测试，不是出题人故意的.

题目链接：

COGS
Luogu

Ac code ：

#include <cstdio>
#include <iostream>
#include <queue>
#include <cstring> 
#define il inline
using namespace std;
const int inf=0x7fffffff;
const int maxm=110000;
int head[maxm],to[maxm*2],cap[maxm*2],net[maxm*2],deep[maxm],cnt=1;
il void add(int x,int y,int c){cnt++,to[cnt]=y,cap[cnt]=c,net[cnt]=head[x],head[x]=cnt;}
queue <int> dl;
int vis[maxm];
il bool BFS(int s,int t)
{
    while(!dl.empty()) dl.pop();
    memset(deep,-1,sizeof(deep));
    dl.push(s),deep[s]=0;
    while(!dl.empty())
    {
        int x=dl.front();dl.pop();
        for(int i=head[x];i;i=net[i])
         if(cap[i]>0&&deep[to[i]]==-1)
          dl.push(to[i]),deep[to[i]]=deep[x]+1;
    }
    return deep[t]==-1?0:1;
}
int dfs(int now,int flow,int t)
{
    if(now==t) return flow;
    int w,used=0;
    for(int i=head[now];i;i=net[i])
    {
        int v=to[i];
        if(deep[v]==deep[now]+1&&cap[i])
        {
            w=dfs(v,min(flow-used,cap[i]),t);
            cap[i]-=w;
            cap[i^1]+=w;
            used+=w;
            if(used==flow) return flow;
        }
    }
    if(!used) deep[now]=-1;
    return used;
}
il int dinic(int s,int t)
{
    int maxflow=0;
    while(BFS(s,t)) maxflow+=dfs(s,inf,t);
    return maxflow;
}
inline void adx(int x,int y,int cax)
{
    add(x,y,cax),add(y,x,0);
}
inline bool read(int &res)
{
    res=0;char ch=getchar();
    while (!isdigit(ch)) { if (ch == '\r') return false; ch = getchar(); }
    while (isdigit(ch)) { res = res * 10 + ch - '0'; ch = getchar(); }
    if (ch == '\r') return false; else return true;
}
int main()
{
    //freopen("shuttle.in","r",stdin);
    //freopen("shuttle.out","w",stdout);
    int n,m;
    read(n),read(m); 
    int s=0,t=1000,sum=0;
    for(int i=1,c,x;i<=n;i++)
    {
        read(c);
        sum+=c;
        adx(s,i,c);
        while(read(x)) adx(i,x+n,inf);
        adx(i,x+n,inf);
    }
    for(int i=1,c;i<=m;i++) read(c),adx(i+n,t,c);
    int ans=dinic(s,t);
    for(int i=1;i<=n;i++)
     if(deep[i]!=-1) printf("%d ",i);

    puts("");

    for(int i=1;i<=m;i++)
     if(deep[i+n]!=-1) printf("%d ",i);
    puts(""); 
    printf("%d\n",sum-ans);
    return 0;
}