[FZYZOJ 1355] 8-2 太空飞行计划问题

P1355 -- 8-2 太空飞行计划问题

时间限制：1000MS

内存限制：131072KB

Description

 W 教授正在为国家航天中心计划一系列的太空飞行。每次太空飞行可进行一系列商业性实验而获取利润。现已确定了一个可供选择的实验集合E={E1，E2，…，Em}，和进行这些实验需要使用的全部仪器的集合I={I1，I2，…In}。实验Ej需要用到的仪器是I的子集RjÍI。配置仪器Ik的费用为ck美元。实验Ej的赞助商已同意为该实验结果支付pj美元。W教授的任务是找出一个有效算法，确定在一次太空飞行中要进行哪些实验并因此而配置哪些仪器才能使太空飞行的净收益最大。这里净收益是指进行实验所获得的全部收入与配置仪器的全部费用的差额。

对于给定的实验和仪器配置情况，编程找出净收益最大的试验计划。(n,m<50)

Input Format

第1行有2 个正整数m和n。m是实验数，n是仪器数。接下来的m 行，每行是一个实验的有关数据。第一个数赞助商同意支付该实验的费用；接着是该实验需要用到的若干仪器的编号。最后一行的n个数是配置每个仪器的费用。

Output Format

第1 行是实验编号；第2行是仪器编号；最后一行是净收益。

Sample Input

2 3
10 1 2
25 2 3
5 6 7

Sample Output

1 2
1 2 3
17

Hint

【题解+吐槽】

今天IOI day1，恭喜杜教成功AK,rank1，恭喜ChinaTeam总分比USATeam高，3Au1Ag，预祝day2:4AK！

听说T1签到题T2恶心题T3卡时题

听说尼日利亚队集体爆0为了引出大新闻- -

 

回到正事

这题最小割模型。

每个实验像S连边，费用为可得费用，累计进total。

每个仪器像T连边，费用为要花费的费用。

实验要用到的仪器 连边，费用为INF（表示不能割）。

跑最小割就好了

一直以为dinic写错了原来是建图脑洞太大一个地方错了TAT

#include <stdio.h>
#include <queue>
#include <string.h>
#include <iostream>
using namespace std;
const int V=210,E=100010;
int head[V],next[E],to[E],flow[E];
int m,n,S,T,tot=1,total;
bool vis[V];
inline int min(int a,int b) {
    return a<b?a:b;
}
inline void add(int u,int v,int fl) {
    tot++;
    to[tot]=v;
    next[tot]=head[u];
    head[u]=tot;
    flow[tot]=fl;
    tot++;
    to[tot]=u;
    next[tot]=head[v];
    head[v]=tot;
    flow[tot]=0;
}

int c[V];

queue<int> q;
inline void clean() {
    while(!q.empty())
        q.pop();
}

inline bool bfs() {
    memset(c,0,sizeof(c));
    clean();
    q.push(S);
    c[S]=1;
    while(!q.empty()) {
        int top=q.front();
        q.pop();
        for (int i=head[top];i;i=next[i]) {
            int _to=to[i],_flow=flow[i];
            if(_flow==0 || c[_to]>0) continue;
            c[_to]=c[top]+1;
            q.push(_to);
            if(_to==T) return 1;
        }
    }
    return 0;
}

int dfs(int now,int low) {
    if (now==T) return low;
    int flowx,r=low;
    for (int i=head[now];i;i=next[i]) {
        int _to=to[i],_flow=flow[i];
        if(c[_to]!=c[now]+1 || _flow==0) continue;
        flowx=dfs(_to,min(r,_flow));
        r-=flowx;
        flow[i]-=flowx, flow[i^1]+=flowx;
        if(!r) return low;
    }
    if(r==low) c[now]=-1;
    return low-r;
}

inline int dinic() {
    int ansx=0;
    while(bfs()) ansx+=dfs(S,210000000);
    return ansx;
}

int main() {
    scanf("%d%d",&m,&n);
    S=0;
    T=n+m+1;
    for (int i=1,cost;i<=m;++i) {
        char str;
        scanf("%d%c",&cost,&str);
        total+=cost;
        add(S,i,cost);
        while(str!='
') {
            scanf("%d%c",&cost,&str);
            add(i,cost+m,210000000);
        }
    }
    for (int i=1,cost;i<=n;++i) {
        scanf("%d",&cost);
        add(i+m,T,cost);
    }
    int answer=total-dinic();
    for (int i=head[S];i;i=next[i])
        if(flow[i])
            printf("%d ",to[i]);
    printf("
");
    for (int i=head[T];i;i=next[i]) {
        //cout<<i<<' '<<flow[i]<<endl;
        if(flow[i]) printf("%d ", to[i]-m);
    }
    printf("
%d
",answer);
    return 0;
}

网络流建图一直都那么坑TAT