这题刚开始猜出来正确的思路和实现方法,但正是因为猜的,所以才并不能十分的确定,仅仅只是知其然而不知其所以然,
同时虽然思路对了,但原图多构了两排点,这是没有必要的,所以以后建完图后,看看是不是可以简化模型
将每一个点拆成两个点,一个点表示和喜欢的人连,另一个点表示和不喜欢的人连,流量为k,因为最多只能和k个不喜欢的人跳
二分舞曲数目mid,判定是否能举办mid次,将s与各个点连cap为mid的边,最后判断是否满流(是否flow==n*mid?),表示是否每个点都有与之对应的mid个点成为舞伴(贪心,怎么证)
感觉网络流好像就是通过点与边,对各种实际约束进行限制吧!
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<queue>
#include<algorithm>
using namespace std;
const int N=250;
const int inf=0x3f3f3f3f;
int n,id_b[55][2],id_g[55][2],s,t;
bool mp[55][55];
int head[N],cur[N],tot,k;
struct aa
{
int cap,flow,to,pre;
}edge[N*N];
void addedge(int x,int y,int z)
{
edge[++tot].to=y;edge[tot].pre=head[x];edge[tot].cap=z;head[x]=tot;
edge[++tot].to=x;edge[tot].pre=head[y];edge[tot].cap=0;head[y]=tot;
}
int lev[N];
bool bfs()
{
memset(lev,0,sizeof(lev));
lev[s]=1;queue<int> q;q.push(s);
while (!q.empty())
{
int u=q.front();q.pop();
for (int i=head[u];i;i=edge[i].pre)
if (edge[i].cap>edge[i].flow&&!lev[edge[i].to])
{
lev[edge[i].to]=lev[u]+1;
if (edge[i].to==t) return true;
q.push(edge[i].to);
}
}
return false;
}
int dfs(int u,int maxflow)
{
if (u==t||maxflow==0) return maxflow;
int ans=0;
for (int &i=cur[u];i;i=edge[i].pre)
if (lev[edge[i].to]==lev[u]+1)
{
int flow=dfs(edge[i].to,min(maxflow,edge[i].cap-edge[i].flow));
ans+=flow;
maxflow-=flow;
edge[i].flow+=flow;
edge[((i-1)^1)+1].flow-=flow;
if (maxflow==0) return ans;
}
return ans;
}
int work(int li)
{
int ans=0;
memset(head,0,sizeof(head));tot=0;
memset(edge,0,sizeof(edge));
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
if (mp[i][j]) addedge(id_b[i][0],id_g[j][0],1);
else addedge(id_b[i][1],id_g[j][1],1);
for (int i=1;i<=n;i++)
{
addedge(s,id_b[i][0],li);
addedge(id_b[i][0],id_b[i][1],k);
addedge(id_g[i][0],t,li);
addedge(id_g[i][1],id_g[i][0],k);
}
while (bfs())
{
for (int i=s;i<=t;i++) cur[i]=head[i];
ans+=dfs(s,inf);
}
return ans;
}
int main()
{
scanf("%d%d",&n,&k);
char ch[55];
int bj=0;
for (int i=1;i<=n;i++)
id_b[i][0]=++bj,id_b[i][1]=++bj,
id_g[i][0]=++bj,id_g[i][1]=++bj;
s=0,t=bj+1;//
for (int i=1;i<=n;i++)
{
scanf("%s",ch);
for (int j=0;j<strlen(ch);j++)
if (ch[j]=='Y') mp[i][j+1]=true;
}
int l=0,r=n,mid,ans=0;
while (l<=r)
{
mid=(l+r)>>1;
if (work(mid)==mid*n) {ans=mid;l=mid+1;}
else r=mid-1;
}
printf("%d",ans);
return 0;
}
好像有贪心方法
一下引自
BLADEVIL大神
另外还有一种更强大的贪心方法,首先对于二分图,假设每个点的度数都是k,那么这个图一定有k个不同的完备匹配,我们可以知道,假设二分图只有两个点,那么显然正确,4个点的情况,左2,右2这种情况也正确,那么假设对于2*n的二分图,我们新添加左右各一个节点,这样左面的节点可以和右面的前n个节点连接,原来前n各节点连接的节点可以给新添加的右面的点相连,所以有n+1个不同的完备匹配,根据数学归纳可证。那么假设没有喜欢的关系,每个人都可以和k个人跳舞,这样答案就是k,那么现在有喜欢的限制关系了,假设对于i这个人,他有a个喜欢的人,那么他一定可以和这a个人每个人跳一次舞,那么对于不喜欢的人,最多跳k次,那么这个人可以跳a+k次,我们就可以求出来所有人喜欢的人的数量,然后加上k,和n去min就好了。假设存在一种情况喜欢的人最少的内个人条的舞的次数不是最后的答案的话,这时这个人肯定可以和另一个人跳舞,那么除了他之外的所有人都至少有一个人可以跳舞,这样就可以再增加一次答案,因为喜欢是双向的,那么对于每一个人和不喜欢的人跳一次舞的时候,双方都会减少一次和不喜欢的人跳舞的次数。