bzoj1019題解
阿新 • • 發佈:2017-06-02
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 #defineView Coderange(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 longret=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 }
bzoj1019題解