欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页  >  IT编程

bzoj3438 小M的作物

程序员文章站 2022-07-11 11:18:07
题目链接 思路 先考虑没有额外收益的时候怎么做。 从$S$向第$i$点连一条容量为$A_i$边,表示种在$A$中的收益。 从第$i$个点向$T$连一条容量为$B_i$的边,表示种在$B$中的收益。 然后求出来最小割,用总收益减去即可。 完成之后如下图: 然后考虑如何处理额外收益 对于每一个额外的收益 ......

题目链接

思路

先考虑没有额外收益的时候怎么做。

\(s\)向第\(i\)点连一条容量为\(a_i\)边,表示种在\(a\)中的收益。

从第\(i\)个点向\(t\)连一条容量为\(b_i\)的边,表示种在\(b\)中的收益。

然后求出来最小割,用总收益减去即可。

完成之后如下图:

bzoj3438 小M的作物

然后考虑如何处理额外收益

对于每一个额外的收益,我们先新建一个点\(x\),表示全部建在\(a\)的收益。

需要保证如果不割掉这条边,那么就说明这些点全部建在\(a\).所以这些点向\(t\)的连边必须全部割掉。

那么就从\(s\)\(x\)连边,从\(x\)向该组合中所有点连边。

对于全部建在b处的额外收益,同理。

加入上图中含有组合\(\{1,2\}\),建图之后如下:

bzoj3438 小M的作物

然后用总收益减去最小割就是答案。

代码

/*
* @author: wxyww
* @date:   2019-06-10 08:33:40
* @last modified time: 2019-06-10 09:10:44
*/
#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include<ctime>
using namespace std;
typedef long long ll;
const int n = 10100,m = 5000100,inf = 2e9;
ll read() {
    ll x=0,f=1;char c=getchar();
    while(c<'0'||c>'9') {
        if(c=='-') f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9') {
        x=x*10+c-'0';
        c=getchar();
    }
    return x*f;
}
struct node {
    int v,nxt,w;
}e[m];
queue<int>q;
int dep[n],head[n],ejs = 1;
void add(int u,int v,int w) {
    e[++ejs].v = v;e[ejs].w = w;e[ejs].nxt = head[u];head[u] = ejs;
    e[++ejs].v = u;e[ejs].w = 0;e[ejs].nxt = head[v];head[v] = ejs;
}
int s,t,cur[n];
bool bfs() {
    memset(dep,0,sizeof(dep));
    while(!q.empty()) q.pop();
    dep[s] = 1;q.push(s);
    while(!q.empty()) {
        int u = q.front();q.pop();
        for(int i = head[u];i;i = e[i].nxt) {
            int v = e[i].v;
            if(!dep[v] && e[i].w) {
                dep[v] = dep[u] + 1;
                q.push(v);
                if(v == t) return true;
            }
        }
    }
    return false;
}
int dfs(int u,int now) {
    if(u == t) return now;
    int ret = 0;
    for(int &i = cur[u];i;i = e[i].nxt) {
        int v = e[i].v;
        if(dep[v] == dep[u] + 1 && e[i].w) {
            int k = dfs(v,min(now - ret,e[i].w));
            e[i].w -= k;e[i ^ 1].w += k;
            ret += k;
            if(ret == now) return ret;
        }
    }
    return ret;
}
int dinic() {
    int ret = 0;
    while(bfs()) {

        memcpy(cur,head,sizeof(cur));
        ret += dfs(s,inf);
    }
    return ret;
}
int ans = 0;
int main() {
    int n = read();
    s = n + 1,t = s + 1;
    for(int i = 1;i <= n;++i) {
        int k = read();ans += k;add(s,i,k);
    }
    for(int i = 1;i <= n;++i) {
        int k = read();ans += k;add(i,t,k);
    }
    int m = read();
    int cnt = n + 2;
    for(int i = 1;i <= m;++i) {
        int k = read(),c1 = read(),c2 = read();
        int x = ++cnt,x_ = ++cnt;
        ans += c1,ans += c2;
        add(s,x,c1);add(x_,t,c2);
        for(int j = 1;j <= k;++j) {
            int x = read();
            add(x,x,inf);add(x,x_,inf);
        }
    }
    cout<<ans - dinic();
    return 0;
}