scu oj 4442 Party(2015年四川省acm程序设计竞赛)
Party
n frogs
are invited to a tea party. Frogs are conveniently numbered by 1,2,…,n.
The tea party has black and green tea in service. Each frog has its own preference. He or she may drink only black/green tea or accept both.
There are m pairs
of frogs who dislike each other. They fight when they are serving the same type of tea.
Luckily, frogs can be divided into 2 groups
such that no two frogs in the same group dislike each other.
Frogs like gems. If the i-th
frog can be paid wi gems
instead of serving tea, it will not fight with others anymore.
The party manager has to dicide how to serve tea/pay gems to avoid fights, minimizing the total gems paid.
Input
The input consists of multiple tests. For each test:
The first line contains 2 integers n,m (1≤n≤103,0≤m≤104).
The second line contains n integers w1,w2,…,wn.
(1≤wi≤106).
The third line contains n integers p1,p2,…,pn.
(1≤pi≤3). pi=1 means
the i-th
frog drinks only black tea. pi=2 means
it drinks only green one, while pi=3means
it accepts both.
Each of the following m lines
contains 2 integers ai,bi,
which denotes frog ai and bi dislike
each other. (1≤ai,bi≤n)
Output
For each test, write 1 integer
which denotes the minimum total gems paid.
Sample Input
2 1
1 1
3 3
1 2
2 1
1 1
2 2
1 2
3 2
2 1 2
1 3 2
1 2
2 3
Sample Output
0
1
1
然后,yy了几个晚上。最终发现果然是神一样的建图。
如今m对青蛙之间有矛盾。有矛盾的青蛙他们不能喝一样的茶,对于每一仅仅青蛙。能够给他w[i]金币。让他不喝茶。他就不会和不论什么青蛙矛盾了。最少须要给多少金币让他们之间没有矛盾存在。这题里面另一个非常重要的信息,我一開始没有看出来。
这句话————“Luckily,
frogs can be divided into 2 groups
such that no two frogs in the same group dislike each other.”事实上这句话非常关键,它的意思就是这个矛盾关系图不存在奇环,是一个2分图。
我们如今考虑仅仅能喝一种茶的青蛙,红和绿之间即使存在矛盾也没有影响,那么仅仅用考虑红和红,绿和绿之间的矛盾关系。如果不考虑2种都能喝的青蛙。那么红和红,绿和绿都是单独的求出最大点权独立集。
因为是2分图,能够非常easy求。
如今考虑一仅仅2种茶都能够喝的青蛙。这里是关键,我的解法是把这个青蛙拆成2个点,一个点代表他喝红茶,一个点代表他喝绿茶,首先这2个点是不能同一时候存在。那么必定有边,然后如果这仅仅青蛙和其它仅仅喝红茶的青蛙有矛盾关系。那么相应就是这个青蛙喝红茶和那个青蛙喝红茶不能同一时候存在,其它情况事实上是类似的分析。这样拆点把矛盾关系图建立出来以后,实际上还是一个2分图。如今就是选出最大点权的点让这些点之间不存在矛盾关系。
实际上就是点权最大独立集。转换成求点权最小覆盖的模型,就能非常轻松的求解了。
。。。。想这么久,果然是神建图。事实上主要还是拆点。
- //#pragma comment(linker, "/STACK:1024000000,1024000000")
- #include<cstdio>
- #include<cstring>
- #include<algorithm>
- #include<iostream>
- #include<cmath>
- #include<queue>
- using namespace std;
- const int mmax = 2010;
- const int inf = 0x3fffffff;
- struct node
- {
- int flow;
- int en;
- int next;
- }E[100*mmax];
- int p[mmax];
- int num;
- void init()
- {
- memset(p,-1,sizeof p);
- num=0;
- }
- void add(int st,int en,int flow)
- {
- E[num].en=en;
- E[num].flow=flow;
- E[num].next=p[st];
- p[st]=num++;
- E[num].en=st;
- E[num].flow=0;
- E[num].next=p[en];
- p[en]=num++;
- }
- int d[mmax];
- int cur[mmax];
- bool vis[mmax];
- bool bfs(int st,int en)
- {
- memset(vis,0,sizeof vis);
- d[st]=0;
- vis[st]=1;
- queue<int>q;
- q.push(st);
- while(!q.empty())
- {
- int x=q.front();
- q.pop();
- for(int i=p[x]; i+1; i=E[i].next)
- {
- int v=E[i].en;
- if(!vis[v]&&E[i].flow)
- {
- vis[v]=1;
- d[v]=d[x]+1;
- q.push(v);
- }
- }
- }
- return vis[en];
- }
- int dfs(int st,int en,int flow)
- {
- if(st==en||flow==0)
- return flow;
- int f=0,dd;
- for(int &i=cur[st]; i+1;i=E[i].next)
- {
- int v=E[i].en;
- if(d[st]+1==d[v]&&(dd=dfs(v,en,min(flow,E[i].flow)))>0)
- {
- E[i].flow-=dd;
- E[i^1].flow+=dd;
- flow-=dd;
- f+=dd;
- if(flow==0)
- break;
- }
- }
- return f;
- }
- int dinic(int st,int en,int n)
- {
- int flow=0;
- while(bfs(st,en))
- {
- for(int i=0;i<=n;i++)
- cur[i]=p[i];
- flow+=dfs(st,en,inf);
- }
- return flow;
- }
- vector<int>e[mmax];
- int fg[mmax];
- int w[mmax];
- void dfs(int u)
- {
- int SZ=e[u].size();
- for(int i=0;i<SZ;i++)
- {
- int v=e[u][i];
- if(d[v]==-1)
- {
- d[v]=d[u]^1;
- dfs(v);
- }
- }
- }
- int main()
- {
- int n,m;
- while(~scanf("%d %d",&n,&m))
- {
- for(int i=0;i<mmax;i++)
- e[i].clear();
- for(int i=1;i<=n;i++)
- scanf("%d",&w[i]);
- int sum=0;
- for(int i=1;i<=n;i++)
- {
- scanf("%d",&fg[i]);
- if(fg[i]==3)
- sum+=w[i];
- }
- for(int i=0;i<m;i++)
- {
- int u,v;
- scanf("%d %d",&u,&v);
- e[u].push_back(v);
- e[v].push_back(u);
- }
- init();
- memset(d,-1,sizeof d);
- for(int i=1;i<=n;i++)
- {
- if(d[i]==-1)
- {
- d[i]=0;
- dfs(i);
- }
- }
- int st=0,en=2*n+1;
- for(int i=1;i<=n;i++)
- {
- if(d[i]==0)
- {
- if(fg[i]==1)
- add(st,i,w[i]);
- if(fg[i]==2)
- add(i,en,w[i]);
- if(fg[i]==3)
- {
- add(st,i,w[i]);
- add(i,i+n,inf);
- add(i+n,en,w[i]);
- }
- }
- else
- {
- if(fg[i]==1)
- add(i,en,w[i]);
- if(fg[i]==2)
- add(st,i,w[i]);
- if(fg[i]==3)
- {
- add(i,en,w[i]);
- add(st,i+n,w[i]);
- add(i+n,i,inf);
- }
- }
- }
- for(int u=1;u<=n;u++)
- {
- int Sz=e[u].size();
- for(int i=0;i<Sz;i++)
- {
- int v=e[u][i];
- if(d[u]==0 && d[v]==1)
- {
- if( (fg[u]==fg[v]) && (fg[u]!=3) )
- {
- if(fg[u]==1)
- add(u,v,inf);
- else
- add(v,u,inf);
- }
- else if( (fg[u]==fg[v]) && (fg[u]==3) )
- {
- add(u,v,inf);
- add(v+n,u+n,inf);
- }
- else if(fg[u]==3 || fg[v]==3 )
- {
- if(fg[u]==3)
- {
- if(fg[v]==1)
- add(u,v,inf);
- else
- add(v,u+n,inf);
- }
- if(fg[v]==3)
- {
- if(fg[u]==1)
- add(u,v,inf);
- else
- add(v+n,u,inf);
- }
- }
- }
- }
- }
- printf("%d\n",dinic(st,en,en)-sum);
- }
- return 0;
- }
scu oj 4442 Party(2015年四川省acm程序设计竞赛)的更多相关文章
- 华南师大 2017 年 ACM 程序设计竞赛新生初赛题解
题解 被你们虐了千百遍的题目和 OJ 也很累的,也想要休息,所以你们别想了,行行好放过它们,我们来看题解吧... A. 诡异的计数法 Description cgy 太喜欢质数了以至于他计数也需要用质 ...
- SCU 4436 Easy Math 2015年四川省赛题
题目链接:http://acm.scu.edu.cn/soj/problem/4436/ 题意:给你n个整数,求这n个数的平方根和是否是一个整数: 解题思路:如果这题每个数给他算出来,必然费时间,可能 ...
- 2010山东省第一届ACM程序设计竞赛
休眠了2月了 要振作起来了!!... http://acm.sdut.edu.cn/sdutoj/problem.php?action=showproblem&problemid=2155 因 ...
- UPC 2224 Boring Counting ★(山东省第四届ACM程序设计竞赛 tag:线段树)
[题意]给定一个长度为N的数列,M个询问区间[L,R]内大于等于A小于等于B的数的个数. [题目链接]http://acm.upc.edu.cn/problem.php?id=2224 省赛的时候脑抽 ...
- 哈尔滨工业大学(威海)第九届ACM程序设计竞赛 Virtual Youtuber
链接 [https://ac.nowcoder.com/acm/contest/624/G] 题意 其实题意说的辣鸡死了,没有说明确. y is the subsequences that its s ...
- 桂林电子科技大学第三届ACM程序设计竞赛 G 路径
链接:https://ac.nowcoder.com/acm/contest/558/G来源:牛客网 小猫在研究树. 小猫在研究路径. 给定一棵N个点的树,每条边有边权,请你求出最长的一条路径,满足经 ...
- (java保留n位小数)precise math function 北京信息科技大学第十届ACM程序设计竞赛 第2题
precise math function Time Limit : 3000/1000ms (Java/Other) Memory Limit : 65535/32768K (Java/Othe ...
- 重复 桂林电子科技大学第三届ACM程序设计竞赛
题目链接:https://ac.nowcoder.com/acm/contest/558/B import java.util.HashSet; import java.util.Scanner; p ...
- 相聚 桂林电子科技大学第三届ACM程序设计竞赛
题目链接:https://ac.nowcoder.com/acm/contest/558/D 就是求有多少块区域,用DFS就可以解决,一遇到一个1就从其开始深搜,将其所在的区域块覆灭(变为0),再遇到 ...
随机推荐
- 使用cmd查看Tensorboard的生成图
代码中"graph_practice7/"文件存放的是生成的神经网络图 现在要查看这个图,使用cmd 进入文件指定目录(我的文件存储地址是E:\workspace\PycharmP ...
- $("[lay-id='"+this.id+"']")
$("[lay-id='"+this.id+"']") $("[lay-id='"+this.id+"'] .layui-tabl ...
- B2. Concurrent 线程池(Executor)
[概述] 与数据库连接管理类似,线程的创建和销毁会耗费较大的开销,使用 “池化技术” 来更好地利用当前线程资源,减少因线程创建和销毁带来的开销,这就是线程池产生的原因. [无限创建线程的不足] 在生产 ...
- Myeclipse下载
地址:http://www.jb51.net/softs/593889.html#download
- Java中创建对象的内存图
所有人都知道面向对象思想,Java中的对象的创建在内存中是如何创建的,传智播客的视频看了一遍,把一些讲解的比较清晰的内容记录下来,方便记忆的更加深刻,Java中创建对象的过程,首先要理解JVM中栈.堆 ...
- python爬取网页图片
# html:网页地址 def getImg2(html): soup = BeautifulSoup(html, 'html.parser') href_regex = re.compile(r'^ ...
- 创建Tensor
目录 创建Tensor numpy, list numpy list zeros, ones, fill zeros ones fill random 打乱idx后,a和b的索引不变 constant ...
- Maximun product
Given a sequence of integers S = {S1, S2, ..., Sn}, you shoulddetermine what is the value of the max ...
- UVA 140 Bandwidth (dfs 剪枝 映射)
题意: 给定一个n个结点的图G和一个结点的排列, 定义结点i的带宽b(i)为i和相邻结点在排列中的最远距离, 所有b(i)的最大值就是这个图的带宽, 给定G, 求让带宽最小的结点排列. 给定的图 n ...
- Codeforces Round #261 (Div. 2) E (DP)
E. Pashmak and Graph Pashmak's homework is a problem about graphs. Although he always tries to do hi ...