[NOI2009 变换序列]

[关键字]:图论 二分图匹配

[题目大意]:太麻烦,就不说了

//=========================================================================

[分析]:因为每个距离都对应着两个数字——就是+d或-d,这样0~n-1和它所转换的序列可以看成一个二分图的模型,而一个完备匹配就对应着一个序列。只要求一遍最大匹配看是否等于n。而字典序的问题,因为用匈牙利算法找增广路就是找到第一个没有被匹配或和它匹配的点可以找到另外一条增广路的点,这样每次都是将已经匹配过的点所匹配的点编号变大,如果倒着找匹配就恰好能保证序号越小就越可以找到编号小的点匹配。

[代码]:

View Code
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;

const int MAXN=20000;

int n;
int m1[MAXN],m2[MAXN];
bool v[MAXN];
vector<int> e[MAXN];

void Init()
{
int t1,t2,d;
scanf("%d",&n);
for (int i=0;i<n;++i)
{
scanf("%d",&d);
t1=i+d;
if (t1>=n) t1%=n;
t2=i-d;
if (t2<0) t2+=n;
if (t1>t2) swap(t1,t2);
e[i].push_back(t1),e[i].push_back(t2);
}
}

bool DFS(int u)
{
for (vector<int>::iterator i=e[u].begin();i<e[u].end();++i)
if (!v[*i])
{
v[*i]=1;
if (m2[*i]==-1 || DFS(m2[*i]))
{
m1[u]=*i;
m2[*i]=u;
return 1;
}
}
return 0;
}

void Solve()
{
memset(m1,255,sizeof(m1));
memset(m2,255,sizeof(m2));
int ans=0;
//for (int i=0;i<n;++i)
for (int i=n-1;i>=0;--i)
{
memset(v,0,sizeof(v));
if (DFS(i)) ++ans; else break;
}
if (ans!=n) printf("No Answer\n");
else
{
printf("%d",m1[0]);
for (int i=1;i<n;++i) printf(" %d",m1[i]);
}
}

int main()
{
Init();
Solve();
system("pause");
return 0;
}



你可能感兴趣的:(序列)