码迷,mamicode.com
首页 > 其他好文 > 详细

[bzoj1485] [HNOI2009]有趣的数列

时间:2016-02-12 22:05:29      阅读:267      评论:0      收藏:0      [点我收藏+]

标签:

  找规律题。。。

  首先观察样例解释可得。我们只要确定奇数位的情况就行了,并且只要i<=第i个奇数<=2*i-1就是合法的= =

  然后我就一直在找规律..................

  最后弃疗跑去看题解....才发现答案不就是卡特兰数吗。。为啥我看半天都没看出来= =

  由搜索引擎可得,h(n)=C(2*n,n)/(n+1)。。。求组合数的话我写了分解质因数...

  顺便学(chao)习(xi)了一下黄学长的分解质因数。。我以前写的都是什么鬼QAQ

  写了个快速幂竟然比直接乘慢orz

技术分享
 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cstring>
 4 #define ll long long
 5 using namespace std;
 6 const int maxn=2000023;
 7 bool cant[maxn];
 8 int pri[maxn],num,mn[maxn],sum[maxn];
 9 ll ans;
10 int i,j,k,n,m,sxt;
11 inline void prerun(int x){
12     register int i,j,k;
13     for(i=2;i<=x;i++){
14         if(!cant[i])pri[++num]=i,mn[i]=num;
15         for(k=pri[j=1]*i;k<=x&&j<=num;k=pri[++j]*i){
16             cant[k]=1,mn[k]=j;
17             if(!(i%pri[j]))break;
18         }
19     }
20 }
21 inline void upd(int x,int v){
22     while(x>1)
23         sum[mn[x]]+=v,x/=pri[mn[x]];
24 }
25 inline int power(ll a,int b){
26     ll c=1;
27     while(b){
28         if(b&1){c*=a;if(c>=sxt)c%=sxt;}
29         a*=a;if(a>=sxt)a%=sxt;
30         b>>=1;
31     }
32     return c;
33 }
34 int main(){
35     scanf("%d%d",&n,&sxt);
36     prerun(n<<1);
37     for(i=n<<1;i>n;i--)upd(i,1);
38     for(;i;i--)upd(i,-1);
39     upd(n+1,-1);
40     ans=1;
41     for(i=1;i<=num;i++){
42         ans*=power(pri[i],sum[i]);
43         if(ans>=sxt)ans%=sxt;
44     }
45     printf("%lld\n",ans);
46     return 0;
47 }
View Code

其实那个合法的条件也可以变一下,将数列里的每个数都减去(i-1)之后就是1<=第i个奇数<=i 。。。就是卡特兰数了>_<

[bzoj1485] [HNOI2009]有趣的数列

标签:

原文地址:http://www.cnblogs.com/czllgzmzl/p/5187299.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!