hdu3277
题意:
女选男,结婚游戏,可以选自己喜欢的,也可以选朋友喜欢的,每次换男的结婚,还可以容忍没有选的男生k次 ,问最多能玩几次(sum)
分析
二分匹配的改编
二分匹配是源点到女孩是1,女孩到男孩是1,男孩到汇点是1
之前是源点到女孩是sum,女孩到男孩是1,男孩到汇点是sum(满流,每一个女孩都可以换至少sum次)
这是是源点到女孩是sum,女孩到男孩是1,女孩可以选k个不喜欢的男孩玩,新建一个源点,每个女孩到自己源点是k,自己源点到自己不喜欢男孩是1,男孩到汇点是sum(满流,每一个女孩都可以换至少sum次)
dinic(多路)
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll INF=1e18;
const int maxn=805;
const int maxm=maxn*maxn;
struct edge
{
int v,next;
ll w;
};
edge no[maxm];
set<int >se[305];
set<int >se_[305];
int vis[305];
int fa[305];
int head[maxn];
int level[maxn];
int cnt;
queue <int> q;
int find_(int x)
{
return x==fa[x] ?x: fa[x]=find_(fa[x]);
}
bool bfs(int s,int t,int n)
{
for(int i=0;i<=n;i++)
level[i]=-1;
while(!q.empty())
{
q.pop();
}
q.push(s);
level[s]=0;
while(!q.empty())
{
int u=q.front();
if(u==t)
return true;
q.pop();
for(int i=head[u];i;i=no[i].next)
{
int v=no[i].v;
ll w=no[i].w;
if(w>0&&level[v]==-1)
{
level[v]=level[u]+1;
q.push(v);
}
}
}
return false;
}
ll dfs(int u,int t,ll max_)
{
if(u==t)
return max_;
ll sum=0;
for(int i=head[u];i;i=no[i].next)
{
edge & e=no[i];
if(e.w>0&&level[e.v]==level[u]+1)
{
ll f=dfs(e.v,t,min(max_-sum,e.w));
sum+=f;
no[i].w-=f;
no[1+((i-1)^1)].w+=f;
if(sum==max_)
break;
}
}
if(!sum)
level[u]=-1;
return sum;
}
ll dinic(int s,int t,int n)
{
ll ans=0;
while(bfs(s,t,n))
{
ll f;
while((f=dfs(s,t,INF))>0)
ans+=f;
}
return ans;
}
void init()
{
for(int i=0;i<255;i++)
{
se[i].clear();
se_[i].clear();
fa[i]=i;
}
cnt=1;
for(int i=0;i<maxn;i++)
{
head[i]=0;
}
//memset(vis,0,sizeof(vis));
}
void init2()
{
cnt=1;
for(int i=0;i<maxn;i++)
{
head[i]=0;
}
}
void add(int x,int y,ll w)
{
no[cnt].v=y;
no[cnt].w=w;
no[cnt].next=head[x];
head[x]=cnt;
cnt++;
}
void solve()
{
int n,m,f,k;
cin>>n>>m>>k>>f;
int s=0,t=3*n+1;
for(int i=1;i<=m;i++)
{
int x,y;
cin>>x>>y;
se[x].insert(y);
}
for(int i=1;i<=f;i++)
{
int x,y;
cin>>x>>y;
int fax=find_(x);
int fay=find_(y);
if(fax!=fay)
{
fa[fax]=fay;
for(auto it:se[fax])
{
se[fay].insert(it);
}
}
}
for(int i=1;i<=n;i++)
{
memset(vis,0,sizeof(vis));
int fai=find_(i);
for(auto it:se[fai])
{
vis[it]=1;
}
for(int j=1;j<=n;j++)
{
if(!vis[j])
{
se_[i].insert(j);
}
}
}
int l=0,r=n;
while(l<r)
{
ll res=0;
int mid=(l+r+1)/2;
for(int i=1;i<=n;i++)
{
int fai=find_(i);
add(s,i,mid);
add(i,s,0);
res+=mid;
for(auto it:se[fai])
{
add(i,it+2*n,1);
add(it+2*n,i,0);
}
for(auto it:se_[i])
{
add(i+n,it+2*n,1);
add(it+2*n,i+n,0);
}
add(i,i+n,k);
add(i+n,i,0);
}
for(int i=1;i<=n;i++)
{
add(i+2*n,t,mid);
add(t,i+n*2,0);
}
ll ans=dinic(s,t,3*n+1);
if(ans==res)
{
l=mid;
}
else
{
r=mid-1;
}
init2();
}
cout<<l<<"\n";
}
int main()
{
ios::sync_with_stdio(0);
cin.tie(0),cout.tie(0);
int __=1;
cin>>__;
while(__--)
{
init();
solve();
}
return 0;
}
isap
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll INF=1e18;
const int maxn=255;
int fa[maxn];
struct Edge
{
int u,v;
};
Edge e[maxn*maxn];
struct edge
{
int v,next;
ll w;
};
int cnt,start,End,ans,head[maxn*3];
ll gap[maxn*3],dis[maxn*3];
edge no[maxn*9*maxn];
int map_[maxn][maxn];
void init()
{
memset(head,-1,sizeof(head));
cnt=0;
}
int find_(int x)
{
return x==fa[x]?x:fa[x]=find_(fa[x]);
}
void add(int x,int y,ll w)
{
no[cnt].v=y;
no[cnt].w=w;
no[cnt].next=head[x];
head[x]=cnt;
cnt++;
}
ll dfs(int u,ll minflow)
{
if(u==End)return minflow;
int i,v;
ll f,min_dis=ans-1,flow=0;
for(i=head[u];i!=-1;i=no[i].next)
{
v=no[i].v;
if(no[i].w<=0)continue;
if(dis[v]+1==dis[u])
{
f=dfs(v,no[i].w>minflow-flow?minflow-flow:no[i].w);
no[i].w-=f;
no[i^1].w+=f;
flow+=f;
if(flow==minflow)break;
if(dis[start]>=ans)return flow;
}
min_dis=min_dis>dis[v]?dis[v]:min_dis;
}
if(flow==0)
{
if(--gap[dis[u]]==0)
dis[start]=ans;
dis[u]=min_dis+1;
gap[dis[u]]++;
}
return flow;
}
ll isap()///ans等于End+1
{
ll maxflow=0;
memset(gap,0,sizeof(gap));
memset(dis,0,sizeof(dis));
gap[0]=ans;
while(dis[start]<ans)
maxflow+=dfs(start,INF);
return maxflow;
}
void solve()
{
int n,m,k,f;
cin>>n>>m>>k>>f;
for(int i=0;i<=n;i++)
{
fa[i]=i;
}
for(int i=1;i<=m;i++)
{
cin>>e[i].u>>e[i].v;
}
for(int i=1;i<=f;i++)
{
int x,y;
cin>>x>>y;
int fax=find_(x);
int fay=find_(y);
if(fax!=fay)
{
fa[fax]=fay;
}
}
for(int i=1;i<=n;i++)
{
fa[i]=find_(i);
}
for(int i=1;i<=m;i++)
{
int u=find_(e[i].u);
int v=e[i].v;
map_[u][v]=1;
}
int l=0,r=n;
start=0;End=3*n+1;
ans=End+1;
ll bns;
while(l<=r)
{
init();
int mid=(l+r)/2;
/*for(int i=1;i<=n;i++)
{
add(start,i,mid);
add(i,start,0);
}*/
for(int i=1;i<=n;i++)
{
add(start,i,mid);
add(i,start,0);
for(int j=1;j<=n;j++)
{
if(map_[fa[i]][j])
{
add(i,j+2*n,1);
add(j+2*n,i,0);
}
else
{
add(i+n,j+2*n,1);
add(j+2*n,i+n,0);
}
}
add(i,i+n,k);
add(i+n,i,0);
}
for(int i=1;i<=n;i++)
{
add(i+2*n,End,mid);
add(End,i+2*n,0);
}
if(isap()==mid*n)
{
bns=mid;
l=mid+1;
}
else
{
r=mid-1;
}
}
cout<<bns<<"\n";
}
int main()
{
ios::sync_with_stdio(0);
cin.tie(0),cout.tie(0);
int __=1;
cin>>__;
while(__--)
{
memset(map_,0,sizeof(map_));
//init();
solve();
}
return 0;
}