【bzoj1562】【[NOI2009]变换序列】匈牙利算法的性质利用

（上不了p站我要死了，侵权度娘背锅）

Description

Input

Output

Sample Input
5
1 1 2 2 1
Sample Output
1 2 4 0 3
HINT
30%的数据中N≤50；
60%的数据中N≤500；
100%的数据中N≤10000。

（要不是我看不懂它在变什么。。。）
实际上就是对其中每个元素进行加di mod n或者是减di加n后mod n的操作。

那么每个元素都对应两个元素。而新的序列要求是排列，即不能有相同的。想到二分图匹配，就可得到一一对应的关系。

至于如何输出最小字典序，就要看对匈牙利算法的熟悉度和理解深度了。首先是考虑越靠前的数字要尽量对应小的数，而匈牙利在跑的时候是优先匹配先扫描到的边。如果先让前面的数字匹配小的数，之后在“腾”的时候就可能让前面的数的匹配变大，这样就不优了。但如果让前面的数字先匹配大的数，之后也不一定让它与小的数配对。
那该怎么办呢。。。
于是又发现，越后匹配的其实优先级越高，在扫边的时候会优先用先扫到的边进行匹配。所以这道题就从后往前for进行匹配，每次优先for小的边（这个就和建边的顺序有关了）

AC（PE*1）代码（要注意这道题要求了行末无空格）

#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>#define ll long longusing namespace std;template <typename T>inline void read(T &res){    T k=1,x=0;char ch=0;    while(ch<'0'||ch>'9'){if(ch=='-')k=-1;ch=getchar();}    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}    res=x*k;}const int N=10000+5;int n,d[N];int head[N],to[N*2],nxt[N*2],hh=0;bool vis[N];int bl[N],has[N];void adde(int a,int b){    hh++;    to[hh]=b;    nxt[hh]=head[a];    head[a]=hh;}bool find(int u){    if(u==-1) return true;    for(int i=head[u];i;i=nxt[i]){        int v=to[i];        if(vis[v]) continue;        vis[v]=1;        if(bl[v]==-1||find(bl[v])){            bl[v]=u;            has[u]=v;            return true;        }    }    return false;}int main(){    read(n);    for(int i=0;i<n;i++) read(d[i]);    for(int i=0;i<n;i++){        int x=(i+d[i])%n,y=(i-d[i]+n)%n;        if(x>y) swap(x,y);        adde(i,y),adde(i,x);    }    int cnt=0;    memset(bl,-1,sizeof(bl));    for(int i=n-1;i>=0;i--){        memset(vis,0,sizeof(vis));        if(find(i)) cnt++;    }    if(cnt<n){        printf("No Answer\n");        return 0;    }    for(int i=0;i<n-1;i++){        printf("%d ",has[i]);    }    printf("%d",has[n-1]);    return 0;}