日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

P5004-专心OI - 跳房子【dp,矩阵乘法】

發布時間:2023/12/3 编程问答 48 豆豆
生活随笔 收集整理的這篇文章主要介紹了 P5004-专心OI - 跳房子【dp,矩阵乘法】 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

正題

題目鏈接:https://www.luogu.org/problemnew/show/P5004


題目大意

NNN個無色格子排成一行,可以把某些格子染成黑色,但兩個黑色格子之間必須至少有MMM個無色格子,求方案數


解題思路

首先很明顯
fn=∑i=0n?m?1fif_n=\sum_{i=0}^{n-m-1}f_ifn?=i=0n?m?1?fi?

然后我們發現∑i=0n?m?1fi\sum_{i=0}^{n-m-1}f_ii=0n?m?1?fi?就是前綴和,那我們在矩陣乘法中維護一個前綴和
A=[∑i=0n?m?1fifn?mfn?m+1...fn?1fn]A=\begin{bmatrix}\sum_{i=0}^{n-m-1}f_i \\ f_{n-m} \\ f_{n-m+1} \\ ... \\ f_{n-1} \\ f_{n} \end{bmatrix}A=?????????i=0n?m?1?fi?fn?m?fn?m+1?...fn?1?fn???????????
然后看著轉移就好了。


codecodecode

#include<cstdio> #include<algorithm> #include<cstring> #define ll long long using namespace std; const ll Size=20,XJQ=1e9+7; struct matrix{ll a[Size][Size]; }f; ll n,m; matrix operator *(matrix a,matrix b) {matrix c;memset(c.a,0,sizeof(c.a));for(ll i=0;i<=m;i++)for(ll j=0;j<=m;j++)for(ll k=0;k<=m;k++)(c.a[i][j]+=a.a[i][k]*b.a[k][j]%XJQ)%=XJQ;return c; } matrix power(matrix f,ll b) {b--;matrix ans;memset(ans.a,0,sizeof(ans.a));for(int i=0;i<=m;i++)ans.a[0][i]=1;ans.a[0][m]++;while(b){if(b&1) ans=ans*f;f=f*f;b>>=1;}return ans; } int main() {scanf("%lld%lld",&n,&m);for(ll i=1;i<=m;i++)f.a[i][i-1]=1;f.a[0][m]=1;f.a[m][m]=1;f=power(f,n);printf("%lld",f.a[0][m]); }

總結

以上是生活随笔為你收集整理的P5004-专心OI - 跳房子【dp,矩阵乘法】的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。