BZOJ1562: [NOI2009]变换序列

因为每个位置i给出了距离di，所以他最后的位置只有
$(i+di)Mod\ n,(i-di+n)Mod\ n$ 2（或1）种选择，因为是二分图，所以可以用匈牙利算法。
因为要求最后的序列字典序最小，所以按照升序建立邻接表，n~1倒序匹配Ti。
证明：根据匈牙利算法的匹配过程，这是显然的…..

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

inline void swap(int &x,int &y){x^=y;y^=x;x^=y;}
inline int read()
{
    char c; int x;
    while(!((c=getchar())>='0'&&c<='9'));
    x=c-'0';
    while((c=getchar())>='0'&&c<='9') (x*=10)+=c-'0';
    return x;
}
const int maxn = 11000;

struct edge
{
    int y,nex;
    edge(){}
    edge(const int &_y,const int &_nex){y=_y;nex=_nex;}
}a[maxn<<2]; int len,fir[maxn],n;
inline void ins(const int &x,const int &y){a[++len]=edge(y,fir[x]);fir[x]=len;}

int match[maxn],t[maxn],tp;
bool v[maxn];
bool find_(const int x)
{
    if(v[x]) return false;
    t[++tp]=x;
    v[x]=true;
    for(int k=fir[x];k;k=a[k].nex)
    {
        const int y=a[k].y;
        if(!match[y]||find_(match[y]))
        {
            match[y]=x;
            return true;
        }
    }
    return false;
}
int ans[maxn];

int main()
{
    memset(fir,0,sizeof fir); len=0;

    n=read();
    for(int i=1;i<=n;i++)
    {
        int x=read();
        int t1=i-x>0?i-x:n-(x-i);
        int t2=i+x>n?i+x-n:i+x;
        if(t1>t2) swap(t1,t2);
        ins(i,t2); ins(i,t1);
    }

    bool flag=true;
    for(int i=n;i>=1;i--)
    {
        tp=0; flag=find_(i); if(!flag) break;
        while(tp) v[t[tp--]]=false;
    }
    if(!flag) {printf("No Answer\n"); 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); printf("%d\n",ans[n]-1);

    return 0;
}