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),再遇到 ...
随机推荐
- controller,sevices层,java初步了解
一.controller层 二.service层 1.接口 2.接口的实现 转换 ClearingAccountArgument对象
- bash基础——管道符、通配符
1.多命令顺序执行 多命令顺序执行 格式 作用 ; 命令1 ; 命令2 多个命令之间没有任何逻辑联系 && 命令1&&命令2 逻辑与 当命令1正确执行,则命令2才会执行 ...
- ajax 传参数 数组 会加上中括号
解决办法 1,由于版本过高导致 我用的是1.9版本 2, 有三种选择. 一种是JSON.stringify([1,2,3]),到后端再解析. 另外一种是后端的接受的contentType改成appli ...
- idea搭建SSM的maven项目(tomcat容器)
一.创建maven的web项目 (1)选择项目的骨架 (2)写项目的坐标 (3)maven的设置 设置maven的本地仓库,以及配置文件的位置,同时点击+号,填入archetypeCatalog和in ...
- UVA 221 城市化地图(离散化思想)
题意: 给出若干个栋楼俯视图的坐标和面积,求从俯视图的南面(可以视为正视图)看过去到底能看到多少栋楼. 输入第一个n说明有n栋楼,然后输入5个实数(注意是实数),分别是楼的左下角坐标(x,y), 然后 ...
- Django——分页功能Paginator
Django分页功能----Paginator Paginator所需参数: Paginator(object_list,per_page) Paginator常用属性: per_page: 每页显示 ...
- HDU 4945 (dp+组合数学)
2048 Problem Description Teacher Mai is addicted to game 2048. But finally he finds it's too hard to ...
- noip模拟赛 补兵
分析:比较难想的一道dp题.要想补兵的数量最多,最后每个小兵的血量肯定是呈一个阶梯状的:i,i+1,i+2......i+k.那么记录一下每个血量i离它最近的小兵的血量是多少,记作cur[i].那么把 ...
- 食物(bzoj 3280)
Description 明明这次又要出去旅游了,和上次不同的是,他这次要去宇宙探险! 我们暂且不讨论他有多么NC,他又幻想了他应该带一些什么东西.理所当然的,你当然要帮他计算携带N件物品的方案数. 他 ...
- msp430入门编程01
msp430单片机最小系统 msp430入门学习 msp430入门编程