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

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 >

【DP】方格计数(nowcoder 20107-B)

發布時間:2023/12/3 30 豆豆
生活随笔 收集整理的這篇文章主要介紹了 【DP】方格计数(nowcoder 20107-B) 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

正題

nowcoder 20107-B


題目大意

在一個H*W的平面上,選擇N個點(坐標為非負整數),使得每個點之間的距離大于D,問你有多少種方案


解題思路

fi,j,kf_{i,j,k}fi,j,k?為有i個格子,選擇k個格子,距離最少為j的方案數,預處理f

然后枚舉斜率,得出以各個點為起點的最大長度,然后乘上方案數


code

#include<cmath> #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #define ll long long #define N 510 #define wyc 1000000007 using namespace std; ll T,n,m,L,s,nn,g,len,ans,d[N][N],f[N][N][100]; int main() {for(ll i=1;i<=500;++i)d[i][0]=d[0][i]=i;for(ll i=1;i<=500;++i)for(ll j=1;j<=i;++j)d[i][j]=d[j][i]=d[j][i%j];for(ll j=1;j<=500;++j){f[0][j][0]=1;for(ll i=1;i<=500;++i){for(ll k=0;k<=50;++k)f[i][j][k]=f[i-1][j][k];if(i>=j)for(ll k=1;k<=50;++k)(f[i][j][k]+=f[i-j][j][k-1])%=wyc;}}scanf("%lld",&T);while(T--){scanf("%lld%lld%lld%lld",&s,&n,&m,&L);ans=0;if(s==1){printf("%lld\n",(n+1)*(m+1)%wyc);continue;}for(ll i=0;i<=n;++i)for(ll j=0;j<=m;++j){if(d[i][j]!=1)continue;if(i==0||j==0)g=1;else g=2;//水平翻轉也計算一遍len=ceil(L/sqrt(i*i+j*j));for(ll k=1;k;++k){if(i*k>n||j*k>m)break;(ans+=(n-i*k+1)*(m-j*k+1)%wyc*(f[k][len][s-1]-f[k-1][len][s-1])%wyc*g%wyc)%=wyc;//(n-i*k+1)*(m-j*k+1)為有多少點可以以此為起點,f[k][len][s-1]-f[k-1][len][s-1]為每個點求最長的,如果不是最長要減去}}printf("%lld\n",(ans+wyc)%wyc);}return 0; }

總結

以上是生活随笔為你收集整理的【DP】方格计数(nowcoder 20107-B)的全部內容,希望文章能夠幫你解決所遇到的問題。

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