BZOJ 3438: 小M的作物 最大权闭合子图
2017-08-24 23:13
441 查看
3438: 小M的作物
Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 1403 Solved: 607
[Submit][Status][Discuss]
Description
小M在MC里开辟了两块巨大的耕地A和B(你可以认为容量是无穷),现在,小P有n中作物的种子,每种作物的种子有1个(就是可以种一棵作物)(用1...n编号),现在,第i种作物种植在A中种植可以获得ai的收益,在B中种植可以获得bi的收益,而且,现在还有这么一种神奇的现象,就是某些作物共同种在一块耕地中可以获得额外的收益,小M找到了规则中共有m种作物组合,第i个组合中的作物共同种在A中可以获得c1i的额外收益,共同总在B中可以获得c2i的额外收益,所以,小M很快的算出了种植的最大收益,但是他想要考考你,你能回答他这个问题么?Input
第一行包括一个整数n第二行包括n个整数,表示ai第三行包括n个整数,表示bi第四行包括一个整数m接下来m行,
对于接下来的第i行:第一个整数ki,表示第i个作物组合中共有ki种作物,
接下来两个整数c1i,c2i,接下来ki个整数,表示该组合中的作物编号。输出格式
Output
只有一行,包括一个整数,表示最大收益Sample Input
34 2 1
2 3 2
1
2 3 2 1 2
Sample Output
11样例解释A耕地种1,2,B耕地种3,收益4+2+3+2=11。
1<=k< n<= 1000,0 < m < = 1000 保证所有数据及结果不超过2*10^9。
最大权闭合子图
#include<cmath> #include<ctime> #include<cstdio> #include<cstring> #include<cstdlib> #include<complex> #include<iostream> #include<algorithm> #include<iomanip> #include<vector> #include<string> #include<bitset> #include<queue> #include<map> #include<set> using namespace std; typedef long long ll; inline int read() { int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();} return x*f; } void print(int x) {if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');} const int N=3010,inf=0X3f3f3f3f; int n,D,F,S=0,T=N-1; int ecnt=1,last ; struct EDGE{int to,nt,val;}e[N*N]; inline void readd(int u,int v,int val) {e[++ecnt]=(EDGE){v,last[u],val};last[u]=ecnt;} inline void add(int u,int v,int val) {readd(u,v,val);readd(v,u,0);} int d ,q ; bool bfs() { memset(d,0,sizeof(d)); register int head=0,tail=1,u,i; q[head]=0;d[0]=1; while(head<tail) { u=q[head++]; for(i=last[u];i;i=e[i].nt)if(e[i].val&&!d[e[i].to]) {d[e[i].to]=d[u]+1;q[tail++]=e[i].to;} } return d[T]; } int dfs(int u,int lim) { if(u==T|!lim)return lim; register int i,res=0,tmp; for(i=last[u];i;i=e[i].nt) if(e[i].val&&d[u]+1==d[e[i].to]) { tmp=dfs(e[i].to,min(e[i].val,lim)); lim-=tmp;res+=tmp;e[i].val-=tmp;e[i^1].val+=tmp; if(!tmp)d[e[i].to]=-1;if(!lim)break; } return res; } ll flow; void dinic() {while(bfs())flow+=dfs(S,inf);} int main() { n=read(); register int u,val,c1,c2,k,i;ll ans=0; for(i=1;i<=n;++i){val=read();add(S,i,val);ans+=val;} for(i=1;i<=n;++i){val=read();add(i,T,val);ans+=val;} int m=read(); for(i=1;i<=m;++i) { k=read();c1=read();c2=read();ans+=c1+c2; add(i+m+n,T,c2);add(S,i+n,c1); while(k--){u=read();add(u,i+m+n,inf);add(i+n,u,inf);} } dinic(); print(int(ans-flow));puts(""); return 0; } /* 3 4 2 1 2 3 2 1 2 3 2 1 2 11 */
相关文章推荐
- 【BZOJ 3438】小M的作物 最大权闭合子图
- bzoj 3438: 小M的作物 (最大权闭合子图)
- 【最大权闭合子图/最小割】BZOJ3438-小M的作物【待填】
- 【BZOJ-3438】小M的作物 最小割 + 最大权闭合图
- BZOJ 3438 小M的作物 最大权闭合图
- BZOJ 3438 最大权闭合子图
- bzoj 1497 [NOI2006]最大获利【最大权闭合子图+最小割】
- BZOJ1565 [NOI2009]植物大战僵尸 【最大权闭合子图 + tarjan缩点(或拓扑)】
- bzoj 1565: [NOI2009]植物大战僵尸 (最大权闭合子图+tarjan+拓扑序)
- [BZOJ1497] 最大权闭合子图+最小割+最大流
- [最大权闭合子图] BZOJ 4873 [Shoi2017]寿司餐厅
- hdu 3879 Base Station bzoj 1497 最大获利问题 最大权闭合子图
- BZOJ 1497 [NOI2006]最大获利 最大权闭合子图
- BZOJ 1497: [NOI2006]最大获利 最大权闭合子图
- bzoj1565 [NOI2009]植物大战僵尸(拓扑序+最大权闭合子图,最小割)
- [省选前题目整理][BZOJ 1565][NOI 2009]植物大战僵尸(最小割+最大权闭合子图建模)
- [BZOJ4873][六省联考2017]寿司餐厅(最大权闭合子图)
- bzoj 1565: [NOI2009]植物大战僵尸 最大权闭合子图+拓扑排序
- bzoj 4823: [Cqoi2017]老C的方块【最大权闭合子图】
- 【BZOJ 1497】 [NOI2006]最大获利 最大权闭合子图