正題

---

題目大意

$n?m$的矩陣，求有多少條路徑的乘積不小于$S$。

---

解題思路

我們用總路徑數減去乘積小于$S$的路徑數

我們很容易想到用$f_{i,j,k}$表示到$(i,j)$這個點，然后乘積之和為$k$的$dp$。但是時間復雜度$O(nmS)$顯然難以勝任本題。

我們考慮將$S?1$整除分塊，用$f_{i,j,k}$表示$(i,j)$這個點時，再乘上一個大于等于$k$的數就會大于等于$S$。

然后我們可以得到動態轉移方程
$$f_{i,j,k}=f_{i?1,j,z}+f_{i,j?1,z}(z=\frac{S?1}{?\frac{S?1}{k}??a_{i,j}})$$

然后$k$只有$2?\sqrt{S}$，所以時間復雜度$O(nm\sqrt{S})$

---

$code$

#include<cstdio> #include<cstring> #include<algorithm> #include<cmath> using namespace std; const int XJQ=1e9+7,N=310; int n,m,s,t,a[N][N],f[2][N][5000],ans,num[5000],v[1100000],c[N][N],cnt; int main() {freopen("mobitel.in","r",stdin);freopen("mobitel.out","w",stdout);scanf("%d%d%d",&n,&m,&s);c[1][0]=1;s--;for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)c[i][j]=(c[i][j-1]+c[i-1][j])%XJQ;for(int i=1,k;i<=s;i=k+1){k=s/(s/i);num[++cnt]=s/i;v[num[cnt]]=cnt;}/*for(int i=s;i>=1;i--)v[i]=v[i]?v[i]:v[i+1];*/for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)scanf("%d",&a[i][j]);f[1][1][v[s/a[1][1]]]=1;for(int i=1;i<=n;i++){memset(f[~i&1],0,sizeof(f[~i&1]));for(int j=1;j<=m;j++)for(int k=1;k<=cnt;k++){int z=num[k];if(!f[i&1][j][k]) continue;if(i<n&&z/a[i+1][j]>0) (f[~i&1][j][v[z/a[i+1][j]]]+=f[i&1][j][k])%=XJQ;if(j<m&&z/a[i][j+1]>0) (f[i&1][j+1][v[z/a[i][j+1]]]+=f[i&1][j][k])%=XJQ;}}for(int k=1;k<=cnt;k++)(ans+=f[n&1][m][k])%=XJQ;printf("%d",(c[n][m]-ans+XJQ)%XJQ); }

總結

以上是生活随笔為你收集整理的jzoj6311-Mobitel【dp,整除分块】的全部內容，希望文章能夠幫你解決所遇到的問題。

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