bzoj4000傳送門
我一上來打了個傻乎乎的狀壓。。成功tle 50%(不要阻止我裝sb。。)
其實這道題敘述有點問題,給的那個3*p的矩陣,第一行是第0行。。。那麼就發現轉移只跟自己上一行的狀態有關,但n太大了,而狀態很少,少到能寫進乙個矩陣,快速冪get。
只要構造出f[i][j],i狀態能轉移到j狀態,則f[i][j]=1;
把這個矩陣自乘n次即可。
#include
#include
#include
#include
#include
#define ll unsigned long long
#define mod (ll)4294967296
using
namespace
std;
int n,m,p,k,tot,w[4],xp[10],g[65];
ll sum;
struct node
return c;
}}ans,h;
inline
void check(int x,int y)
else
}for(int i=1;i<=m;i++)if(xp[i-1]&g[y])
else
}h.f[x][y]=1;
}inline
void get(int x)
else
}g[++tot]=x;
}int yjn()
for(int i=0;ifor(int i=1;i<=tot;i++)
for(int j=1;j<=tot;j++)
check(i,j);
for(int i=1;i<=tot;i++)ans.f[i][i]=1;
while(n)
for(int i=1;i<=tot;i++)sum=(sum+ans.f[1][i])%mod;
cout
int main()
快速冪(矩陣快速冪)
求 3 0 3 1 3 n mod 1000000007 input 輸入乙個數n 0 n 10 9 output 輸出 計算結果 sample input 3sample output 40 分析 利用等比數列的求和公式得所求和是 3 n 1 1 2,如果暴力求3 n 1 會超時,這裡引入快速冪來...
快速冪 矩陣快速冪
快速冪 正常情況下求乙個數的冪時間複雜度為o n 而快速冪能把時間複雜度降到o logn 舉個例子 求5的13次方 思想首先把13化為二進位制 1101,即13 1101 8 1 4 1 2 0 1 1 即5 13 58 1 54 1 52 0 5 1 15 5 8 1 5 4 1 5 2 0 5 ...
快速冪 矩陣快速冪
快速冪 我們求a ba b ab最直接的方法就是把a乘b次這樣的話複雜度就是o n o n o n 但是在比賽時面對1e9的資料時還是會輕鬆超時的,此時就需要一種更快的乘法來幫助我們 我們把b拆成二進位制的形式得到a ba b ab a 10.01 a a1 0.01此時對b分解得到的序列10.01...