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了一个费用流构图: 源连所有点,配对的点连啊,所有点连汇... 后 ...
随机推荐
- [z] Git pull 强制覆盖本地文件
git fetch --all git reset --hard origin/master git pull git fetch 只是下载远程的库的内容,不做任何的合并 git reset 把HEA ...
- TouchSlide 插件使用介绍
TouchSlide(PC端插件http://www.superslide2.com/demo.html#effect1) 可用于javascript触屏滑动特效插件,移动端滑动特效,触屏焦点图,触屏 ...
- power designer 从sqlserver数据库获取字段说明&导出rtf文档模板
具体的操作稍后在修改 附件下载:https://files.cnblogs.com/files/zinan/powerDesigner.rar
- hdu 4004 (二分加贪心) 青蛙过河
题目传送门:http://acm.hdu.edu.cn/showproblem.php?pid=4004 题目意思是青蛙要过河,现在给你河的宽度,河中石头的个数(青蛙要从石头上跳过河,这些石头都是在垂 ...
- HDU 2255.奔小康赚大钱 最大权匹配
奔小康赚大钱 Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Subm ...
- Ubuntu安装R及R包
安装R $sudo apt-get update $sudo apt-get install r-base $sudo apt-get install r-base-dev 安装一些可能的依赖包 $s ...
- 生活类App原型制作分享-AnyList
AnyList是一款可以帮你创建购物清单,并且帮助你整理食谱的生活工具App,前面引导页采用图片+文字的方式,介绍App的用法,登录注册采用选项卡切换的方式,减少了页面切换的繁琐操作,在Mockplu ...
- 新手必看,史上最全的iOS开发教程集锦,没有之一!
最近大火的iPhone XS Max和iPhone XS,不知道有没有同学已经下手了呢?一万三的价位确实让很多人望而却步啊.据说为了赢得中国的用户,专门出了双卡双待的,可想而知中国市场这块“肥肉”人人 ...
- python 字典遍历
dic1={"name":"kxb","age":28}for k,v in dic1.items(): print(k+",,, ...
- 手机(Android)资源
手机型号 API Android版本 Lenovo A238t API 10 2.3.5 华为 P7 API 19 4.4.2