bzoj2893(费用流)
先缩点,然后拆点,其实是很经典的一种操作,把不好做的点拆成边,然后我一开始想的是网络流,答案当然是增广次数,
但可以发现跑网络流的话不同的跑法增广次数不一样,不太好找最小的。我们可以换一种神奇的思路,跑最大费用流,
这样根据费用流每次都在最长路上增广保证每次都跑了尽量多的点,根据贪心原理可知这样是正确的。
详见http://blog.csdn.net/iamzky/article/details/41846687;
(一开始我一直错,发现自己增广用的dinic,模板打惯了就直接打上了。。。。忘了dinic是可以一次dfs实现多次增广的,没法统计增广次数。。我好蠢啊)
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdio>
using namespace std;
const int maxn=,inf=1e9;
struct edg{
int nxt,to,f,c;
}e[maxn],g[maxn];
int res,last[maxn],H[maxn],t1=,t2=,cas,dis[maxn],q[maxn],head,tail;
void add1(int x,int y){
++t1;e[t1].nxt=last[x];last[x]=t1;e[t1].to=y;
}
void add2(int x,int y,int z,int zz){
++t2;g[t2].nxt=H[x];H[x]=t2;g[t2].to=y;g[t2].f=z;g[t2].c=zz;
++t2;g[t2].nxt=H[y];H[y]=t2;g[t2].to=x;g[t2].f=;g[t2].c=-zz;
}
int scc,low[maxn],dfn[maxn],bel[maxn],cnt,sta[maxn],top,in[maxn];
void tarjan(int x){
low[x]=dfn[x]=++cnt;sta[++top]=x;in[x]=;
for(int i=last[x];i;i=e[i].nxt){
int v=e[i].to;
if(!dfn[v]){
tarjan(v);low[x]=min(low[x],low[v]);
}
else{
if(in[v])low[x]=min(low[x],dfn[v]);
} }
//cout<<"orz"<<endl;
if(low[x]==dfn[x]){
int now=;scc++;
while(now!=x){
now=sta[top--];in[now]=;
bel[now]=scc;
}
}
}
struct dui{
int from,to;
}tmp[maxn];
int pe[maxn],pv[maxn],mp[][],ans,vis[maxn],N,n,m,a,b,x,y,S,T,A[maxn],B[maxn],a1[maxn],b1[maxn];
void solve(){
while(){
memset(dis,-,sizeof(dis));
head=tail=;q[++tail]=S;dis[S]=;vis[S]=;
while(head!=tail){
head=(head+)%maxn;
int u=q[head];
for(int i=H[u];i;i=g[i].nxt){
int v=g[i].to;
if(g[i].f&&dis[v]<dis[u]+g[i].c){
dis[v]=dis[u]+g[i].c;
pe[v]=i;pv[v]=u;
if(!vis[v]){
vis[v]=;tail=(tail+)%maxn;
q[tail]=v;
}
}
}
vis[u]=;
}
if(dis[T]<)break;
int mm=inf;
for(int u=T;u!=S;u=pv[u])
mm=min(mm,g[pe[u]].f);
for(int u=T;u!=S;u=pv[u])
g[pe[u]].f-=mm,g[pe[u]^].f+=mm;
res+=dis[T]*mm;
if(dis[T]>)
++ans;
}
}
int main(){
cin>>cas;
while(cas--){
memset(mp,,sizeof(mp));
memset(last,,sizeof(last));
memset(H,,sizeof(H));
memset(dfn,,sizeof(dfn));
cin>>n>>m>>a>>b;
t1=t2=;scc=cnt=;
for(int i=;i<=a;++i){scanf("%d",&A[i]);}
for(int i=;i<=b;++i)scanf("%d",&B[i]);
for(int i=;i<=m;++i){
scanf("%d%d",&tmp[i].from,&tmp[i].to);
add1(tmp[i].from,tmp[i].to);
}
for(int i=;i<=n;++i){
if(!dfn[i])tarjan(i);
}
S=;T=scc*+;
for(int i=;i<=n;++i){
for(int j=last[i];j;j=e[j].nxt){
int u=bel[i],v=bel[e[j].to];
if(mp[u][v]||u==v)continue;
add2(u+scc+,v+,inf,);
mp[u][v]=mp[v][u]=;
}
}
for(int i=;i<=a;++i){
add2(S,+bel[A[i]],inf,);
}
for(int i=;i<=b;++i){
add2(+bel[B[i]]+scc,T,inf,);
}
for(int i=;i<=scc;++i){
add2(+i,+scc+i,,);add2(+i,+scc+i,inf,);
}
res=ans=;solve();
if(res!=scc)puts("no solution");
else{printf("%d\n",ans);}
}
//system("pause");
return ;
}
bzoj2893(费用流)的更多相关文章
- hdu-5988 Coding Contest(费用流)
题目链接: Coding Contest Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Ot ...
- POJ2195 Going Home[费用流|二分图最大权匹配]
Going Home Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 22088 Accepted: 11155 Desc ...
- BZOJ3130: [Sdoi2013]费用流[最大流 实数二分]
3130: [Sdoi2013]费用流 Time Limit: 10 Sec Memory Limit: 128 MBSec Special JudgeSubmit: 960 Solved: 5 ...
- 洛谷 1004 dp或最大费用流
思路: dp方法: 设dp[i][j][k][l]为两条没有交叉的路径分别走到(i,j)和(k,l)处最大价值. 则转移方程为 dp[i][j][k][l]=max(dp[i-1][j][k-1][l ...
- Codeforces 730I [费用流]
/* 不要低头,不要放弃,不要气馁,不要慌张 题意: 给两行n个数,要求从第一行选取a个数,第二行选取b个数使得这些数加起来和最大. 限制条件是第一行选取了某个数的条件下,第二行不能选取对应位置的数. ...
- zkw费用流+当前弧优化
zkw费用流+当前弧优化 var o,v:..] of boolean; f,s,d,dis:..] of longint; next,p,c,w:..] of longint; i,j,k,l,y, ...
- 【BZOJ-4213】贪吃蛇 有上下界的费用流
4213: 贪吃蛇 Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 58 Solved: 24[Submit][Status][Discuss] Desc ...
- 【BZOJ-3638&3272&3267&3502】k-Maximum Subsequence Sum 费用流构图 + 线段树手动增广
3638: Cf172 k-Maximum Subsequence Sum Time Limit: 50 Sec Memory Limit: 256 MBSubmit: 174 Solved: 9 ...
- [bzoj4514]数字配对[费用流]
今年SDOI的题,看到他们在做,看到过了一百多个人,然后就被虐惨啦... 果然考试的时候还是打不了高端算法,调了...几天 默默地yy了一个费用流构图: 源连所有点,配对的点连啊,所有点连汇... 后 ...
随机推荐
- C#执行javascript代码,执行复杂的javascript代码新方式
1. 使用nuget 包"Jurassic", 注意,如果 nuget上的包 用起来出现错误,请自行下载 github代码,自行编译最新代码成dll,再引用. 官方的nuget包 ...
- python pyMysql 自定义异常 函数重载
# encoding='utf8'# auth:yanxiatingyu#2018.7.24 import pymysql __all__ = ['Mymysql'] class MyExcept(E ...
- 史上最全的MSSQL复习笔记
1.什么是SQL语句 SQL语言,结构化的查询语言(Structured Query Language),是关系数据库管理系统的标准语言.它是一种解释语言,写一句执行一句,不需要整体编译执行. 语法特 ...
- hihoCoder1159 扑克牌
一道记忆化搜索 原题链接 和着色方案很像,这里就不详细阐述,可以去我博客里的着色方案里看. 但要注意本题不一样的是同种面值的牌花色不同,所以在转移时还需要乘上同种面值的牌的个数. #include&l ...
- tableView中cell的复用机制
TableView的重用机制,为了做到显示和数据分离,IOS tableView的实现并且不是为每个数据项创建一个tableCell.而是只创建屏幕可显示最大个数的cell,然后重复使用这些cell, ...
- Python编程笔记(第三篇)【补充】三元运算、文件处理、检测文件编码、递归、斐波那契数列、名称空间、作用域、生成器
一.三元运算 三元运算又称三目运算,是对简单的条件语句的简写,如: 简单条件处理: if 条件成立: val = 1 else: val = 2 改成三元运算 val = 1 if 条件成立 else ...
- Https如何确保传输安全的
1.对称加密算法 加密和解密算法是公开的,那个密钥是保密的, 只有双方才知道, 这样生成的加密消息(密文) 别人就无法得知了. 2.非对称加密算法 RSA算法非常有意思,是有一对儿钥匙, 一个是保密的 ...
- Clover相关知识
-f 重建驱动缓存 darkwake=4 有深度睡眠有关的设置,不懂 kext-dev-mode=1 启用第三方驱动,比较重要. dart=0 修复因开启 VT-d 导致系统启动时SMC五国错误,系统 ...
- PSP(3.30——4.5)以及周记录
1.PSP 3.30 12:00 13:00 10 50 Account前端 A Y min 13:00 13:20 0 20 站立会议 A Y min 15:15 17:00 20 85 Accou ...
- Java第2章笔记
1.什么是变量:在程序运行过程中它的值是允许改变的量 2.java中常用的数据类型分为四类八种 第一类:整型 int(整数类型) byte(字节类型) short(短整形) ...