【TOJ】1070. Ouroboros Snake --构造？DFS？欧拉回路？【TOJ】3381. DeBruijin

 

张贴各种代码，仅供我来膜拜。。

 

先贴[TOJ]1070 Ouroboros Snake 的

 

---

 

1.宝哥的构造法，0.00s，在TOJ所有人的代码中排行第三

#include<cstdio> #include<cstring> int va[32800]; bool hs[32800]; int ca(int ans,int n){ if(ans>n/2) ans-=n/2; if(!hs[ans*2]) {hs[ans*2]=1;return ans*2;} else { hs[ans*2+1]=1; return ans*2+1; } } int main(){ int n,k,i,j,ans,t; while(scanf("%d%d",&n,&k)&&(n+k)){ memset(va,-1,sizeof(va[0])*(1<<n)); memset(hs,0,sizeof(hs[0])*(1<<n)); va[0]=0; hs[0]=1; t=ans=1<<(n-1); for(i=(1<<n)-1,j=0;j<n;i--,j++){ va[i]=ans; hs[ans]=1; ans+=(t>>1); t>>=1; } if(va[k]>-1) {printf("%d/n",va[k]);continue;} ans=0; for(i=1;i<=k;i++) ans=ca(ans,(1<<n)); printf("%d/n",ans); } return 0; }

 

---

 

2.我的欧拉回路法，改自我[TOJ]1416 Catenyms的代码 ，1.99s，pku超时，排行最后一名。

#include<stdio.h> #include<string.h> #include<iostream> #include<algorithm> using namespace std; #define MAXN 32770 struct node{ int a,b,num; }our[MAXN]; int path[MAXN],n,t,k,temp; bool f[MAXN],flag[MAXN]; int ans[16][MAXN]; void find_path_d(int now,int& step){ for (int i=0;i<temp;i++) if(our[i].a==now&&!f[i]){ f[i]=1; find_path_d(our[i].b,step); path[step++]=i; } } void save(int cnt,int k){ for(int i=cnt-1,j=0;i>=0;i--,j++) {ans[k][j]=our[path[i]].num;}//printf("%d ",our[path[cnt]].num);} //printf("/n"); } void solve(){ int i,j,cnt; temp=1<<15; for(j=0;j<temp;j++) our[j].num=j; for(i=1;i<16;i++){ memset(f,0,sizeof(f)); temp=1<<i; for(j=0;j<temp;j++){ our[j].a=j>>1; our[j].b=j&(~(1<<(i-1))); //printf("num=%d/n",our[i].num); //printf(" a=%d ",our[i].a);print_bit(our[i].a); //printf(" b=%d ",our[i].b);print_bit(our[i].b); }cnt=0; find_path_d(0,cnt); save(cnt,i); } } /* void print_bit(int num){ int i=0,list[MAXN]; while(num){ list[i++]=num%2; num/=2; } while(i--) printf("%d",list[i]); printf("/n"); }*/ int main(){ int i,j; solve(); while(scanf("%d%d",&n,&k),n||k){ printf("%d/n",ans[n][k]); } //system("pause"); return 0; }

 

---

 

3.DFS法，下述强大解题报告所附代码改编而成，0.00s

#include<stdio.h> #include<string.h> #define MAXN 32770 int ans[16][MAXN],n; int f[MAXN],flag; void dfs(int k,int num){ int next; f[num]=1; ans[n][k]=num; next=(num<<1)%(1<<n); if(flag&&!f[next]) dfs(k+1,next); next=(num<<1)%(1<<n)+1; if(flag&&!f[next]) dfs(k+1,next); f[num]=0; if(k==(1<<n)-1) flag=0; } void solve(){ int i; memset(ans,0,sizeof(0)); for(i=1;i<16;i++){ n=i; memset(f,0,sizeof(f[0])*(1<<n)); flag=1; dfs(0,0); } } int main(){ int i,j,k; solve(); while(scanf("%d%d",&n,&k),n||k){ printf("%d/n",ans[n][k]); } return 0; }

 

---

4.最最强大的DFS+DP+欧拉回路法，虽然内存用得有点多，可是代码中体现的思想令我折服。。改编自Roba的代码，0.00s。

 

#include<iostream> #include<cstring> #include<cstdio> #define MAXN 1 << 15 using namespace std; int pns; int ans[16][MAXN]; int mp[16][MAXN][2],tk[MAXN]; void build(int k) { int i,tmp; for(i=0;i< (1<<k);i++) { tmp=(i & ((1<< (k-1))-1)); mp[k][i][0]=tmp<<1; mp[k][i][1]=mp[k][i][0]+1; } } void dfs(int k,int p) { if(tk[mp[k][p][0]]==0){ tk[mp[k][p][0]]=1; dfs(k,mp[k][p][0]); ans[k][pns++]=mp[k][p][0]; } if(tk[mp[k][p][1]]==0){ tk[mp[k][p][1]]=1; dfs(k,mp[k][p][1]); ans[k][pns++]=mp[k][p][1]; } } void init() { int i; for(i=1;i<=15;i++) { memset(tk,0,sizeof(tk)); pns=0; tk[0]=1; build(i); dfs(i,0); ans[i][pns++]=0; } } int main() { int i; init(); int n,k; while(scanf("%d %d",&n,&k) && ( n || k)) { int kk=1<<n; printf("%d/n",ans[n][kk-k-1]); } return 0; }

 

---

 

这里有篇1070的强大的解题报告：http://hi.baidu.com/zab08/blog/item/914d6a89944952bb0f2444c4.html

 

---

 

然后是【TOJ】3381  DeBruijin

做法与1070类似的，或者就是从那改编而成？

1.我改编自上述代码3，00.02s

#include<stdio.h> #include<string.h> #include<iostream> #include<algorithm> using namespace std; #define MAXN 5000 int ans[12][MAXN],n; bool f[MAXN],flag; void dfs(int k,int num){ int next; f[num]=1; ans[n][k]=num; next=(num<<1)%(1<<n); if(flag&&!f[next]) dfs(k+1,next); next=(num<<1)%(1<<n)+1; if(flag&&!f[next]) dfs(k+1,next); f[num]=0; if(k==(1<<n)-1) flag=0; } void solve(){ int i; memset(ans,0,sizeof(0)); for(i=1;i<12;i++){ n=i; memset(f,0,sizeof(f[0])*(1<<n)); flag=1; dfs(0,0); } } int main(){ int i,j,k; solve(); while(scanf("%d",&k)!=EOF){ printf("%d ",1<<k); for(i=0;i<k;i++) printf("0"); for(i=1;i<(1<<k)-(k-1);i++) printf("%d",ans[k][i]&1); printf("/n"); } return 0; }

 

---

 

2.kobe改编自Roba的，与上面4类似 ，0.00s

#include<iostream> #include<cstring> #include<cstdio> #include<map> #define MAXN 1 << 11 using namespace std; int pns; int tk[MAXN]; int ans[MAXN]; int mat[MAXN][2]; void build(int n) { int i,j; int tmp; for(i=0;i<(1<<n);i++) { tmp=(i & ((1<<(n-1))-1)); mat[i][0]=tmp<<1; mat[i][1]=mat[i][0]+1; } } void dfs(int p) { if(tk[mat[p][0]]==0){ tk[mat[p][0]]=1; dfs(mat[p][0]); ans[pns++]=mat[p][0]; } if(tk[mat[p][1]]==0){ tk[mat[p][1]]=1; dfs(mat[p][1]); ans[pns++]=mat[p][1]; } } int main() { int n; while(scanf("%d",&n)!=EOF) { int k; int i; memset(tk,0,sizeof(tk)); build(n); pns=0; tk[0]=1; dfs(0); ans[pns++]=0; printf("%d ",1<<n); for(i=0;i<n;i++) printf("0"); for(i=pns-2;i>n-2;i--) { k=(1&(ans[i])); k=k+'0'; printf("%c",k); } printf("/n"); } return 0; }