BZOJ 1856: [Scoi2010]字符串 [Catalan数]
生活随笔
收集整理的這篇文章主要介紹了
BZOJ 1856: [Scoi2010]字符串 [Catalan数]
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
1856: [Scoi2010]字符串
Time Limit:?5 Sec??Memory Limit:?64 MBSubmit:?1418??Solved:?790
[Submit][Status][Discuss]
Description
lxhgww最近接到了一個生成字符串的任務,任務需要他把n個1和m個0組成字符串,但是任務還要求在組成的字符串中,在任意的前k個字符中,1的個數不能少于0的個數。現在lxhgww想要知道滿足要求的字符串共有多少個,聰明的程序員們,你們能幫助他嗎?Input
輸入數據是一行,包括2個數字n和mOutput
輸出數據是一行,包括1個數字,表示滿足要求的字符串數目,這個數可能會很大,只需輸出這個數除以20100403的余數?
?
對于100%的數據,保證1<=m<=n<=1000000
?
又是裸題,Catalan數變式...
結果我這個傻逼沒預處理逆元階乘T了一次
?
#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <cmath> using namespace std; typedef long long ll; const int N=2e6+5,MOD=20100403; inline int read(){char c=getchar();int x=0,f=1;while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}return x*f; } int n,m; int inv[N],fac[N],facInv[N]; void ini(int n){inv[1]=1;fac[0]=facInv[0]=1;for(int i=1;i<=n;i++){if(i!=1) inv[i]=(ll)-MOD/i*inv[MOD%i]%MOD;inv[i]+=inv[i]<0?MOD:0;fac[i]=(ll)fac[i-1]*i%MOD;facInv[i]=(ll)facInv[i-1]*inv[i]%MOD;} } inline int C(int n,int m){return (ll)fac[n]*facInv[m]%MOD*facInv[n-m]%MOD; } int main(){freopen("in","r",stdin);n=read();m=read();ini(n+m);printf("%d",(C(n+m,n)-C(n+m,n+1)+MOD)%MOD); }?
?
?
?
?
轉載于:https://www.cnblogs.com/candy99/p/6406710.html
總結
以上是生活随笔為你收集整理的BZOJ 1856: [Scoi2010]字符串 [Catalan数]的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: RenderTextureFormat.
- 下一篇: 迟来的2017年计划