TJU 2944 Mussy Paper 最大权闭合子图

传送门

给你一些东西，  每个东西有一个值，有正有负。 在给一些关系， 选了其中一个物品， 和他有关系的也必须全都选上， 关系是单向的。 问最后的最大价值是多少， 如果小于0输出“   **** ”（记不得是什么了。

网络流， 如果一个东西的价值是正的， 源点向它连边， 权值为val； 如果是负的， 汇点向它连边， 权值为-val。 如果1->2， 就从1向2连边， 权值为inf。 最后就是求一个最小割。

 #include<bits/stdc++.h>
 using namespace std;
 #define pb(x) push_back(x)
 #define ll long long
 #define mk(x, y) make_pair(x, y)
 #define lson l, m, rt<<1
 #define mem(a) memset(a, 0, sizeof(a))
 #define rson m+1, r, rt<<1|1
 #define mem1(a) memset(a, -1, sizeof(a))
 #define mem2(a) memset(a, 0x3f, sizeof(a))
 #define rep(i, a, n) for(int i = a; i<n; i++)
 #define ull unsigned long long
 typedef pair<int, int> pll;
 const double PI = acos(-1.0);
 const double eps = 1e-;
 const int mod = 1e9+;
 const int inf = ;
 const int dir[][] = { {-, }, {, }, {, -}, {, } };
 const int maxn = ;
 int head[maxn*], s, t, num, q[maxn*], dis[maxn], cnt, vis[maxn];
 struct node
 {
     int to, nextt, c;
 }e[maxn*];
 void init() {
     mem1(head);
     num = cnt = ;
     mem(vis);
 }
 void add(int u, int v, int c) {
     e[num].to = v; e[num].nextt = head[u]; e[num].c = c; head[u] = num++;
     e[num].to = u; e[num].nextt = head[v]; e[num].c = ; head[v] = num++;
 }
 int bfs() {
     int u, v, st = , ed = ;
     mem(dis);
     dis[s] = ;
     q[ed++] = s;
     while(st<ed) {
         u = q[st++];
         for(int i = head[u]; ~i; i = e[i].nextt) {
             v = e[i].to;
             if(e[i].c&&!dis[v]) {
                 dis[v] = dis[u]+;
                 if(v == t)
                     return ;
                 q[ed++] = v;
             }
         }
     }
     return ;
 }
 int dfs(int u, int limit) {
     if(u == t)
         return limit;
     int cost = ;
     for(int i = head[u]; ~i; i = e[i].nextt) {
         int v = e[i].to;
         if(e[i].c&&dis[u] == dis[v]-) {
             int tmp = dfs(v, min(limit-cost, e[i].c));
             if(tmp>) {
                 e[i].c -= tmp;
                 e[i^].c += tmp;
                 cost += tmp;
                 if(cost == limit)
                     break;
             } else {
                 dis[v] = -;
             }
         }
     }
     return cost;
 }
 int dinic() {
     int ans = ;
     while(bfs()) {
         ans += dfs(s, inf);
     }
     return ans;
 }
 void dfs(int u) {
     vis[u] = ;
     cnt++;
     for(int i = head[u]; ~i; i = e[i].nextt) {
         int v = e[i].to;
         if(e[i].c&&!vis[v]) {
             dfs(v);
         }
     }
 }
 int main()
 {
     int n;
     while(scanf("%d", &n)&&n) {
         s = , t = n+;
         init();
         int sum = , val, p;
         for(int i = ; i<=n; i++) {
             scanf("%d%d", &val, &p);
             if(val>) {
                 add(s, i, val);
                 sum += val;
             }else if(val<) {
                 add(i, t, -val);
             }
             while(p--) {
                 scanf("%d", &val);
                 add(i, val, inf);
             }
         }
         int ans = dinic();
         if(ans == sum) {
             puts("Refused");
         } else {
             dfs();
             printf("%d %d\n", sum-ans, --cnt);
             for(int i = ; i<=n; i++) {
                 if(vis[i]) {
                     cnt--;
                     printf("%d", i);
                     if(cnt)
                         cout<<" ";
                     else
                         cout<<endl;
                 }
             }
         }
     }
 }