[NOI2009]变换序列

来源:互联网 发布:php开源crm 编辑:程序博客网 时间:2024/05/16 01:05

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

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

首先我觉得

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

那么很明显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#definerep(i,a,b)for(int i=a,tt=b;i<=tt;++i)#definedrep(i,a,b)for(int i=a,tt=b;i>=tt;--i)#defineerep(i,e,x)for(int i=x;i;i=e[i].next)#defineirep(i,x)for(typedef(x.begin()) i=x.begin();i!=x.end();i++)#defineread()(strtol(ipos,&ipos,10))#definesqr(x)((x)*(x))#definepbpush_back#definePSsystem("pause");typedeflong longll;typedefpair<int,int>pii;const int oo=~0U>>1;const double inf=1e100;const double eps=1e-6;string name="transform", in=".in", out=".out";//Varstruct 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;}


原创粉丝点击