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

歡迎訪問 生活随笔!

生活随笔

當(dāng)前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

CF451E-Devu and Flowers【组合计数,容斥】

發(fā)布時(shí)間:2023/12/3 编程问答 32 豆豆
生活随笔 收集整理的這篇文章主要介紹了 CF451E-Devu and Flowers【组合计数,容斥】 小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.

正題

題目鏈接:https://www.luogu.org/problem/CF451E


題目大意

長度為nnn的序列aia_iai?,有如下要求

  • ai∈[0...xi]a_i\in[0...x_i]ai?[0...xi?]
  • ∑i=1nai=s\sum_{i=1}^na_i=si=1n?ai?=s
  • 求序列個(gè)數(shù)。


    解題思路

    考慮容斥,這樣我們就可以將條件轉(zhuǎn)換為ai>xia_i>x_iai?>xi?,然后我們可以讓n?xi?1n-x_i-1n?xi??1這樣就將條件變?yōu)榱?span id="ozvdkddzhkzd" class="katex--inline">ai≥0a_i\geq 0ai?0,然后可以用組合數(shù)插板法計(jì)算,也就是答案為
    ∑SCs?∑xi∈S(xi+1)?1n?1?(?1)∣S∣\sum_{S}C_{s-\sum_{x_i\in S}(x_i+1)-1}^{n-1}*(-1)^{|S|}S?Cs?xi?S?(xi?+1)?1n?1??(?1)S


    codecodecode

    #include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; const ll XJQ=1e9+7; ll n,s,inv[50],x[50],ans; bool v[50]; 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 C(ll n,ll m) {if(n<0||m<0||n<m)return 0;if(n%XJQ==0||m==0)return 1;ll x=1;for(ll i=n-m+1;i<=n;i++)x=x*(i%XJQ)%XJQ;return x*inv[m]%XJQ; } void dfs(ll dep,ll _,bool flag) {ll _n=s+n;for(ll i=1;i<=n;i++)if(v[i]) _n-=x[i]+1;if(flag)ans=((ans+_*C(_n-1,n-1))%XJQ+XJQ)%XJQ;if(dep>n) return;v[dep]=1;dfs(dep+1,-_,1);v[dep]=0;dfs(dep+1,_,0); } int main() {scanf("%lld%lld",&n,&s);inv[1]=1;for(ll i=2;i<=20;i++)inv[i]=inv[i-1]*power(i,XJQ-2)%XJQ;for(ll i=1;i<=n;i++)scanf("%lld",&x[i]);ans=C(s+n-1,n-1);dfs(1,1,0);printf("%lld",ans); }

    總結(jié)

    以上是生活随笔為你收集整理的CF451E-Devu and Flowers【组合计数,容斥】的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。