[NOI2009]变换序列

这道题九点开始看……一直写到现在……

事实证明混乱的代码不要接着写，最好是重新写一遍……

首先我觉得

每个点只能有对应的两个点，这个不是很明显的约束了开头后面的直接就推的出来吗？

那么很明显O(n^2)可以过……

然后我发现这是个每个点最多两条边的二分图

但是没有仔细想，依旧觉得它约束了第一个后面就都推的出来

立马写了一个……模拟二分图匹配……

然后wa……

仔细一想发现MS会出现当前局面没有可以约束的点的情况……比如图中存在两个连通分量……所以不行……

我写的程序无法像匈牙利那样修改增广路,so……

它必须是二分图匹配……

然后有一个字典序的问题……

因为每个点只有两个选择，从后往前暴力可过，假设先选字典序小的有没有完备匹配，有的话就可选……

优化一下就是从后往前，假设先选字典序小的有没有增广路

Linux下面不知道一个什么东西莫名其妙的定义了link……导致我临时replace了一下……

//Lib
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>

#include<iostream>
#include<algorithm>
#include<vector>
#include<string>
#include<queue>
#include<set>
#include<map>
#include<list>
using namespace std;
//Macro
#define	rep(i,a,b)	for(int i=a,tt=b;i<=tt;++i)
#define	drep(i,a,b)	for(int i=a,tt=b;i>=tt;--i)
#define	erep(i,e,x)	for(int i=x;i;i=e[i].next)
#define	irep(i,x)	for(typedef(x.begin()) i=x.begin();i!=x.end();i++)
#define	read()	(strtol(ipos,&ipos,10))
#define	sqr(x)	((x)*(x))
#define	pb	push_back
#define	PS	system("pause");
typedef	long long	ll;
typedef	pair<int,int>	pii;
const int oo=~0U>>1;
const double inf=1e100;
const double eps=1e-6;
string name="transform", in=".in", out=".out";
//Var
struct E
{
	int next,node;
}e[40008];
int n,tot,ans;
int s[10008],c[10008],h[20008];
int llink[20008],to[10008][2];
bool vis[20008];
void add(int a,int b)
{
	e[++tot].next=h[a];e[tot].node=b;h[a]=tot;
}
void Init()
{
	scanf("%d",&n);int x,y;
	rep(i,1,n)scanf("%d",s+i);
	rep(i,1,n)
	{
		x=i-s[i];if(x<=0)x+=n;
		y=i+s[i];if(y>n)y-=n;
		to[i][0]=min(x,y);to[i][1]=max(x,y);
		add(i,x+n);if(x!=y)add(i,y+n);
	}
}
bool DFS(int u)
{
	erep(i,e,h[u])
	{
		int v=e[i].node;
		if(vis[v]==0)
		{
			vis[v]=true;
			if(!llink[v]||DFS(llink[v]))
			{
				llink[v]=u;
				llink[u]=v;
				return true;
			}
		}
	}
	return false;
}
void Work()
{
	rep(i,1,n)
	{
		memset(vis,0,sizeof vis);
		if(DFS(i))ans++;
	}
	if(ans<n){printf("No Answer\n");return;}
	int x,y,tmp;
	drep(i,n,1)
	{
		if(to[i][0]+n!=llink[i])
		{
			x=to[i][0];y=to[i][1];tmp=llink[x+n];
			memset(vis,0,sizeof vis);
			vis[x+n]=true;llink[x+n]=i;llink[y+n]=0;
			if(DFS(tmp))llink[i]=x+n;
			else llink[x+n]=tmp,llink[y+n]=i;
		}
	}
	rep(i,1,n-1)printf("%d ",llink[i]-n-1);
	printf("%d\n",llink[n]-n-1);
}
int main()
{
//	freopen((name+in).c_str(),"r",stdin);
//	freopen((name+out).c_str(),"w",stdout);
	Init();
	Work();
//	PS;
	return 0;
}