题意
左边有n个点标号0n-1,右边有n个点标号0n-1,左边每个点连出两条边向右边。问图中是否存在完美匹配,若能,将左边点的匹配点标号按照左边点的标号大小依次写下。如果有多组方案,输入字典序最小的方案。
$ n leqslant 10000$
题解
一道加深对匈牙利算法理解的题。
首先最简单的判是否有完美匹配很好办。但是这题要求字典序最小。一开始我是想按标号从小到大枚举左边的点,并且邻接表也按标号从小到大排序,然后当前枚举到的点如果尝试匹配一个点,而这个点已经被匹配了,那就暂时不选它。但是这样似乎有点复杂,并且不确定正确性。于是我又想啊想,因为匈牙利的方法是枚举到一个匹配点,就尽力匹配这个点,如果实在不行再放弃。所以我们直接把枚举左边的点改成标号从大到小,那么晚枚举的点,会优先匹配标号小的点。
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=1e4;
int n,Time;
int map[maxn+8][2],match[maxn+8],vis[maxn+8],ans[maxn+8];
int read()
{
int x=0,f=1;char ch=getchar();
for (;ch<'0'||ch>'9';ch=getchar()) if (ch=='-') f=-1;
for (;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0';
return x*f;
}
bool extra(int x)
{
for (int i=0;i<2;i++)
{
int child=map[x][i];
if (vis[child]==Time) continue;
vis[child]=Time;
if ((!match[child])||(extra(match[child])))
{
match[child]=x;
return 1;
}
}
return 0;
}
int main()
{
n=read();
for (int i=1;i<=n;i++)
{
int d=read(),x=((i-d)+n-1)%n+1,y=((i+d)-1)%n+1;
if (x<y) swap(x,y);
map[i][0]=y,map[i][1]=x;
}
for (int i=n;i;i--)
{
Time++;
if (!extra(i))
{
puts("No Answer");
return 0;
}
}
for (int i=1;i<=n;i++) ans[match[i]]=i;
for (int i=1;i<=n;i++) printf("%d ",ans[i]-1);
return 0;
}