首页 > 代码库 > BZOJ4403 序列统计
BZOJ4403 序列统计
本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作。
本文作者:ljh2000
作者博客:http://www.cnblogs.com/ljh2000-jump/
转载请注明出处,侵权必究,保留最终解释权!
题目链接:BZOJ4403
正解:$Lucas$定理+组合数学
解题报告:
考虑这类不降的问题,我们通常把第$i$个数$+i$变成求升序序列问题。
那么权值范围变成了$[l+1,r+n]$,一共$r-l+n$个数,我们要选$i$个数的话就是$C_{r-l+i}^{i}$。
显然答案就是对它求和,简化之后式子变成了$C_{n+R-L+1}^{R-L+1}-1$。
//It is made by ljh2000 //有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。 #include <iostream> #include <cstdlib> #include <cstring> #include <cstdio> #include <cmath> #include <algorithm> #include <ctime> #include <vector> #include <queue> #include <map> #include <set> #include <string> #include <complex> #include <bitset> using namespace std; typedef long long LL; typedef long double LB; const double pi = acos(-1); const int mod = 1000003; int n,len,L,R,jie[mod+12],nj[mod+12]; inline int fast_pow(int x,int y){ int r=1; while(y>0) { if(y&1) r=1LL*r*x%mod; x=1LL*x*x%mod; y>>=1; } return r; } inline int C(int n,int m){ if(n<m) return 0; return 1LL*jie[n]*nj[m]%mod*nj[n-m]%mod; } inline int getint(){ int w=0,q=0; char c=getchar(); while((c<‘0‘||c>‘9‘) && c!=‘-‘) c=getchar(); if(c==‘-‘) q=1,c=getchar(); while (c>=‘0‘&&c<=‘9‘) w=w*10+c-‘0‘,c=getchar(); return q?-w:w; } inline void Init(){ jie[0]=nj[0]=1; for(int i=1;i<mod;i++) jie[i]=1LL*jie[i-1]*i%mod; nj[mod-1]=fast_pow(jie[mod-1],mod-2); for(int i=mod-2;i>=1;i--) nj[i]=1LL*nj[i+1]*(i+1)%mod; } inline int Lucas(int n,int m){ if(n<m) return 0; if(m==0) return 1; if(n<mod && m<mod) return C(n,m); return 1LL*Lucas(n/mod,m/mod)*C(n%mod,m%mod)%mod; } inline void work(){ Init(); int T=getint(); while(T--) { n=getint(); L=getint(); R=getint(); len=R-L+1; printf("%d\n",(Lucas(len+n,n)-1+mod)%mod); } } int main() { #ifndef ONLINE_JUDGE freopen("4403.in","r",stdin); freopen("4403.out","w",stdout); #endif work(); return 0; } //有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。
BZOJ4403 序列统计
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。