1. 程式人生 > >bzoj1019題解

bzoj1019題解

static pen 暴力 一個 clas alt efi one 技術分享

【解題思路】

  對於一個hanoi,知道了各種移動操作的優先級,也就確定了方案。可以證明對於盤子數為N的hanoi,任意移動方案都等價於將數目為N-1的一疊盤子移動k次,並將最小的一個盤子經過b次後移動到目標柱頂端。這樣,hanoi的任一移動方案所需次數都滿足線性遞推式f[n]=k*f[n-1]+b。

  因此我們暴力算出數列f的前幾項(更確切地,前三項即可),然後遞推即可。加上矩乘優化後理論復雜度O(log2n)。

【參考程序】

(嗯,實際上N≤30根本不需要矩乘。。也懶得寫。。)

技術分享
 1 #include <bits/stdc++.h>
 2 #define
range(i,c,o) for(register int i=(c);i<(o);++i) 3 #define dange(i,c,o) for(register int i=(c);i>(o);--i) 4 using namespace std; 5 6 static int N; 7 long long f[35]; 8 string opt[6]; 9 stack<int> stick[3]; 10 11 inline long long solve(const int&x) 12 { 13 long long
ret=0; int las=3; 14 range(i,0,3) for(;!stick[i].empty();stick[i].pop()); 15 dange(i,x,0) stick[0].push(i); 16 for(;stick[1].size()<x&&stick[2].size()<x;++ret) 17 { 18 range(i,0,6) 19 { 20 int fr=opt[i][0]-A,to=opt[i][1]-A; 21 if
( 22 fr!=las&&!stick[fr].empty()&&( 23 stick[to].empty()|| 24 stick[fr].top()<stick[to].top() 25 ) 26 ) 27 { 28 stick[las=to].push(stick[fr].top()); 29 stick[fr].pop(); break; 30 } 31 } 32 } 33 return ret; 34 } 35 36 int main() 37 { 38 ios::sync_with_stdio(0); 39 cin>>N; range(i,0,6) cin>>opt[i]; 40 f[1]=1,f[2]=solve(2),f[3]=solve(3); 41 long long k=(f[3]-f[2])/(f[2]-1),b=f[2]-k; 42 range(i,4,N+1) f[i]=k*f[i-1]+b; 43 return cout<<f[N]<<endl,0; 44 }
View Code

bzoj1019題解