jzoj4800-[GDOI2017模拟9.24]周末晚会【dp,循环重构】
生活随笔
收集整理的這篇文章主要介紹了
jzoj4800-[GDOI2017模拟9.24]周末晚会【dp,循环重构】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目大意
nnn個人的圓桌,可以放男可以放女,然后要求不能有連續kkk個女生坐在一起。求方案總數。
解題思路
先不考慮圓桌,設fi,jf_{i,j}fi,j?表示放了i個人,最前面有連續jjj個女生,動態轉移方程顯然
fi,j=fi?1,j?1(j!=0)f_{i,j}=f_{i-1,j-1}(j!=0)fi,j?=fi?1,j?1?(j!=0)
fi,0=∑j=0nfi?1,jf_{i,0}=\sum_{j=0}^n f_{i-1,j}fi,0?=j=0∑n?fi?1,j?
設rrr表示長度為lll不考慮循環重構的方案數那么有
ri=∑j=1kfi,j?(j+1)r_i=\sum_{j=1}^k f_{i,j}*(j+1)ri?=j=1∑k?fi,j??(j+1)
(j+1j+1j+1)表示圓桌旋轉的方法
然后我們定義qiq_iqi?表示長度為iii且不由循環節構成的方案總數那么有
qi=ri?∑d∣iqdq_i=r_i-\sum_{d|i}q_dqi?=ri??d∣i∑?qd?
然后定義sis_isi?表示長度為iii且包括循環重構的方案那么有
si=∑d∣iqdds_i=\sum_{d|i}\frac{q_d}ozvdkddzhkzdsi?=d∣i∑?dqd??
(循環重復了ddd次)
codecodecode
#include<cstdio> #include<algorithm> #include<cstring> #define ll long long using namespace std; const ll XJQ=1e8+7,N=2010; ll power(ll x,ll b){ll ans=1;while(b){if(b&1) ans=ans*x%XJQ;x=x*x%XJQ;b>>=1;}return ans; } ll t,p[N][N],r[N],q[N],ans,n,k; int main() {scanf("%lld",&t);while(t--){scanf("%lld%lld",&n,&k);memset(p,0,sizeof(p));p[1][0]=1;for(ll i=1;i<=2000;i++)for(ll j=0;j<=min(i-1,k);j++){p[i+1][0]=(p[i+1][0]+p[i][j])%XJQ;p[i+1][j+1]+=p[i][j];}memset(r,0,sizeof(r));for(ll i=1;i<=n;i++)for(ll j=0;j<=k;j++)(r[i]+=p[i][j]*(j+1)%XJQ)%=XJQ;ans=0;for(ll i=1;i<=n;i++){q[i]=r[i];for(ll j=1;j<i;j++)if(i%j==0)q[i]=(q[i]-q[j]+XJQ)%XJQ; }for(ll i=1;i<=n;i++)if(!(n%i)) (ans+=q[i]*power(i,XJQ-2)%XJQ)%=XJQ;if(n<=k) ans++;printf("%lld\n",ans);} }總結
以上是生活随笔為你收集整理的jzoj4800-[GDOI2017模拟9.24]周末晚会【dp,循环重构】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 什么是滤波电容 滤波电容是什么东西
- 下一篇: jzoj6191-[NOI2019模拟2