2015湖南省选集训DAY5——work(BZOJ4177)

原创 2015年07月08日 15:54:57

Description

Mike有一个农场,这个农场n个牲畜围栏,现在他想在每个牲畜围栏中养一只动物,每只动物可以是牛或羊,并且每个牲畜围栏中的饲养条件都不同,其中第i个牲畜围栏中的动物长大后,每只牛可以卖a[i]元,每只羊可以卖b[i]元,为了防止牛羊之间相互影响,Mike找到了m条规律,每条规律给出一个三元组(i, j, k)表示如果第i个围栏和第j个围栏养的是不同的动物,那么Mike就需要花费k的代价请人帮忙处理牛羊之间的影响。不过同时Mike也发现k条特殊的规则(S, a, b),表示如果S中所有牲畜围栏中都养的是动物a,那么Mike可以获得b的额外收入。现在Mike想知道他该在哪些围栏中饲养什么动物才能使得总收益最大,为了简化问题,你只需要输出最大收益。

Input

第一行三个整数n、m、k,表示一共有n个围栏,m条规律,k条规则。

第二行有n个整数,表示a[i]。

第三行有n个整数,表示b[i]。

接下来m行,每行有三个整数(i, j, k)表示一条规则。

再接下来k行,每行一开始有三个整数t、a和b,表示一条规则(S, a, b),其中S的大小为t,接下来

t个整数表示S中的元素(a为0表示全为牛,a为1表示全为羊)。

Output

输出一个整数ans,表示最大收益。

Sample Input

4 2 1

1 2 3 1

2 3 1 2

1 2 3

1 3 2

2 0 100 1 2

Sample Output

108

HINT

对于100的数据,n <= 5000, m <= 5000, k <= 5000, a = 0 or 1。

题解与吐槽:

傻逼网络流
然后听取了tkd的建议拆了点。。。
wa后发现拆点毫无意义而且会影响正确性

从源点向每个点连一条流量为ai的边,每个点向汇点连一条流量为bi的边。于是每一种割法对应一种选择。然后考虑第二种限制,我们向有关系的两点互连流量为k的边,意会一下。第三种的话我们再加上一个新点,如果要求是割掉集合中的点与汇点的边,于是从源点向这个点连流量为收益大小的边,在从这个点向集合中的点连流量无穷大的边;另一种情况同理。于是就完了。

请无视掉那个奇怪的typedef,其实它的含义是tkd神犇

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef double db;
typedef long long tkdsb;

const tkdsb inf=0x3f3f3f3f3f3f;

tkdsb getint()
{
    char c=getchar();
    tkdsb f=1,g=0;
    while(c>'9' || c<'0'){if(c=='-')f=-1;c=getchar();}
    while(c<='9' && c>='0')g=(g<<3)+(g<<1)+c-'0',c=getchar();
    return f*g;
}

const tkdsb maxn=40005;

tkdsb n,m,k;

struct edge{
    tkdsb from,to,cap;
};

tkdsb s,t;

vector<tkdsb> g[maxn];
vector<edge> eds;

void addedge(tkdsb from,tkdsb to,tkdsb cap)
{
    g[from].push_back(eds.size());
    eds.push_back((edge){from,to,cap});
    g[to].push_back(eds.size());
    eds.push_back((edge){to,from,0});
}

tkdsb d[maxn];
queue<tkdsb> q;

bool bfs()
{
    memset(d,-1,sizeof d);
    d[s]=0;
    q.push(s);
    while(!q.empty())
    {
        tkdsb x=q.front();q.pop();
        for(vector<tkdsb>::iterator it=g[x].begin();it!=g[x].end();it++)
        {
            edge e=eds[*it];
            if(d[e.to]==-1 && e.cap>0)
            {
                d[e.to]=d[x]+1;
                q.push(e.to);
            }
        }
    }
    return d[t]!=-1;
}

tkdsb cur[maxn];

tkdsb dfs(tkdsb x,tkdsb f)
{
    if(x==t || f==0)return f;
    tkdsb used=0;
    tkdsb temp;

    for(tkdsb i=cur[x];i<g[x].size();i++)
    {
        edge e=eds[g[x][i]];
        if(d[e.to]==d[x]+1 && e.cap)
        {
            temp=dfs(e.to,min(e.cap,f-used));
            cur[x]=i;
            eds[g[x][i]].cap-=temp;
            eds[(g[x][i])^1].cap+=temp;
            used+=temp;
            if(used==f)return used;
        }
    }
    if(used==0)d[x]=-1;
    return used;
}

tkdsb dinic()
{
    tkdsb res=0;
    while(bfs())
    {
        memset(cur,0,sizeof cur);
        res+=dfs(s,inf);
    }
    return res; 
}

//expected score 100

int main()
{
//  freopen("work.in","r",stdin);
//  freopen("work.out","w",stdout);

    n=getint();
    m=getint();
    k=getint();

    s=0;
    t=maxn-1;

    tkdsb x,y;
    tkdsb z;
    tkdsb res=0;

    for(tkdsb i=1;i<=n;i++)
    {
        x=getint();
        addedge(s,i,x);
        res+=x;
    }
    for(tkdsb i=1;i<=n;i++)
    {
        x=getint();
        addedge(i,t,x);
        res+=x;
    }
    for(tkdsb i=1;i<=m;i++)
    {
        x=getint();
        y=getint();
        z=getint();

        addedge(x,y,z);
        addedge(y,x,z);
    }

    tkdsb temp;

    for(tkdsb i=1;i<=k;i++)
    {
        temp=getint();
        x=getint();
        y=getint();
        res+=y;
        if(x==1)
        {
            addedge(i+n,t,y);
            for(tkdsb j=1;j<=temp;j++)
            {
                z=getint();
                addedge(z,i+n,inf);
            }
        }
        else
        {
            addedge(s,i+n,y);
            for(tkdsb j=1;j<=temp;j++)
            {
                z=getint();
                addedge(n+i,z,inf);
            }
        }       
    }
    res-=dinic();
    printf("%lld\n",res);

    return 0;
}

2015湖南省队集训DAY8——梦工厂

梦工厂 (yume.cpp/c/pas)Time Limit: 1 s Memory Limit: 128 M 问题描述「有时候用烂了的名字也会别有深意」 ——摘自EN 语录 “这里是...
  • Starry__Night
  • Starry__Night
  • 2015年07月11日 20:37
  • 591

湖南集训Day5

T1果然是回文自动机(树),zyh讲了一下然而还是听不懂。 下午的网络流看各路神犇秒题,分分钟口胡出来,真是6的不行啊,我这种弱渣只能看看,顺便围观一堆小号233333。 晚上好像整栋楼都在刷Go...
  • nlj1999
  • nlj1999
  • 2015年12月31日 00:00
  • 185

[清华集训2015]V

题目大意给出一个序列a,现在有五种操作 1、把[l,r]每个ai加上x 2、把[l,r]每个ai变成x 3、把[l,r]每个ai变成max(ai-x,0) 4、询问目前ai的值 5、询问ai...
  • WerKeyTom_FTD
  • WerKeyTom_FTD
  • 2016年10月31日 16:36
  • 744

#164. 【清华集训2015】V

线段树维护 ? 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 ...
  • sxb_201
  • sxb_201
  • 2016年03月13日 20:33
  • 311

uoj164. 【清华集训2015】V

Picks博士观察完金星凌日后,设计了一个复杂的电阻器。为了简化题目,题目中的常数与现实世界有所不同。 这个电阻器内有编号为 1∼n1∼n 的 nn 个独立水箱,水箱呈圆柱形,底面积为 1 m21 ...
  • CRZbulabula
  • CRZbulabula
  • 2017年01月09日 11:05
  • 456

数据结构优化连边

什么是线段树连边对于一类题目,它们建出来的图好大好大,这时候通常可以利用线段树(大多数时候是可持久化线段树)来进行优化连边,保证点的连通性不变,路径权值信息不变。例题1:湖南集训2014Jabberw...
  • WerKeyTom_FTD
  • WerKeyTom_FTD
  • 2016年01月15日 19:55
  • 548

[霍尔定理]「2017 山东一轮集训 Day2」LOJ 6062——PAIR

什么是霍尔定理这是关于二分图的一个非常重要的定理。 但是博主今天第一接触,真不愧是蒟蒻。 主要内容: 设M(U1)M(U1)为与U1U1中的点相连的点集,一个二分图U,V(|U|=|X| 感谢...
  • CHN_JZ
  • CHN_JZ
  • 2017年07月26日 22:15
  • 1228

福建省队集训被虐记——DAY3

昨天没写……今天补上吧
  • zhb1997
  • zhb1997
  • 2014年07月15日 22:55
  • 562

2015关于第十一届"蓝狐网络杯"湖南省大学生计算机程序设计竞赛的总结

------------------------------------------------------------------------人生总是需要我们不断受挫,不断总结。然后继续前行!   ...
  • u010579068
  • u010579068
  • 2015年08月31日 12:09
  • 1785

uoj#164. 【清华集训2015】V

传送门:http://uoj.ac/problem/164 思路:科学的题面: 请你写一个数据结构支持以下功能: 1:区间[l,r]加x 2:区间[l,r]减x并和0取max 3:区间覆盖 ...
  • thy_asdf
  • thy_asdf
  • 2016年04月26日 10:37
  • 540
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:2015湖南省选集训DAY5——work(BZOJ4177)
举报原因:
原因补充:

(最多只允许输入30个字)