BZOJ1856[Scoi2010]字符串——组合数学+容斥
生活随笔
收集整理的這篇文章主要介紹了
BZOJ1856[Scoi2010]字符串——组合数学+容斥
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
題目描述
lxhgww最近接到了一個生成字符串的任務,任務需要他把n個1和m個0組成字符串,但是任務還要求在組成的字符串中,在任意的前k個字符中,1的個數不能少于0的個數。現在lxhgww想要知道滿足要求的字符串共有多少個,聰明的程序員們,你們能幫助他嗎?輸入
輸入數據是一行,包括2個數字n和m輸出
輸出數據是一行,包括1個數字,表示滿足要求的字符串數目,這個數可能會很大,只需輸出這個數除以20100403的余數樣例輸入
2 2樣例輸出
2提示
【數據范圍】
對于30%的數據,保證1<=m<=n<=1000
對于100%的數據,保證1<=m<=n<=1000000
?
將選$1$看成往右走,選$0$看成往上走,那么要求的就是從$n*m$的網格的左下角走到右上角且不能穿過$y=x$的方案數。
將不能穿過$y=x$看成不能走到$y=x+1$,答案就是總方案數(即沒有不能穿過$y=x$限制的方案數)-走到$y=x+1$的方案數。
將起點關于$y=x+1$對稱到$(-1,1)$,那么走到$y=x+1$的方案數就是從$(-1,1)$走到$(n,m)$只能往右和往上走的方案數。
最終答案就是$C_{n+m}^{n}-C_{n+m}^{n+1}$,注意當$n<m$時答案為$0$。
#include<set> #include<map> #include<queue> #include<stack> #include<cmath> #include<vector> #include<bitset> #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> using namespace std; const int mod=20100403; int n,m; int fac[2000010]; int inv[2000010]; int C(int n,int m) {return 1ll*fac[n]*inv[m]%mod*inv[n-m]%mod; } int main() {scanf("%d%d",&n,&m);inv[0]=inv[1]=fac[0]=fac[1]=1;for(int i=2;i<=n+m;i++){fac[i]=1ll*fac[i-1]*i%mod;inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;}for(int i=2;i<=n+m;i++){inv[i]=1ll*inv[i-1]*inv[i]%mod;}if(n>=m){printf("%d",(C(n+m,n)-C(n+m,n+1)+mod)%mod);}else{printf("0");} }轉載于:https://www.cnblogs.com/Khada-Jhin/p/10954692.html
總結
以上是生活随笔為你收集整理的BZOJ1856[Scoi2010]字符串——组合数学+容斥的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 突然感觉好失落的说说,心情莫名失落的句子
- 下一篇: 坚持做一件事情真的很不容易,首先必须要克