这是道好题呀,果然NOI里出好题,先跪一下题解。
https://www.byvoid.com/blog/noi-2009-transform/
神犇说得非常详细了,其实就是倒着匹配,然后贪心的去最小的,因为如果不可以的话,前面的会就行调整,所以这个算法是正确的。趁机接着道题学习了一下匈牙利算法,真的是精炼巧妙呀。
匈牙利算法的主体是寻找增广路,这里的增广路不同于网络流中的增广路,从一个未匹配的顶点出发,到匹配的顶点,再到未匹配的顶点,……,最后以未匹配的顶点结尾。这样的一条路径,显然将它翻转是更优的,然后就没有了。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#define maxn 10010
using namespace std;
int n,m,num;
int to[2*maxn],next[2*maxn],head[maxn];
int ans[maxn],lk[maxn];
bool vis[maxn];
int dis(int x,int y)
{
return min(abs(x-y),n-abs(x-y));
}
void addedge(int x,int y)
{
num++;to[num]=y;next[num]=head[x];head[x]=num;
}
bool find(int x)
{
for (int p=head[x];p;p=next[p])
if (!vis[to[p]])
{
vis[to[p]]=1;
if (lk[to[p]]==0 || find(lk[to[p]]))
{
lk[to[p]]=x;
return 1;
}
}
return 0;
}
int main()
{
scanf("%d",&n);
for (int i=0;i<n;i++)
{
int x;
scanf("%d",&x);
int a=i-x+n,b=i+x;
a%=n;b%=n;
if (a<b) swap(a,b);
if (dis(a,i)==x) addedge(i,a);
if (dis(b,i)==x) addedge(i,b);
}
for (int i=n-1;i>=0;i--)
{
memset(vis,0,sizeof(vis));
if (!find(i))
{
printf("No Answer\n");
return 0;
}
}
for (int i=0;i<=n-1;i++) ans[lk[i]]=i;
for (int i=0;i<n-1;i++) printf("%d ",ans[i]);printf("%d\n",ans[n-1]);
return 0;
}