题目描述
我们称一个长度为2n的数列是有趣的，当且仅当该数列满足以下三个条件：
(1)它是从1到2n共2n个整数的一个排列{ai}；
(2)所有的奇数项满足a1<a3<...<a2n-1，所有的偶数项满足a2<a4<...<a2n；
(3)任意相邻的两项a2i-1与a2i(1<=i<=n)满足奇数项小于偶数项，即：a2i-1<a2i。
现在的任务是：对于给定的n，请求出有多少个不同的长度为2n的有趣的数列。因为最后的答案可能很大，所以只要求输出答案 mod P的值。
输入输出格式
输入格式：
输入文件只包含用空格隔开的两个整数n和P。输入数据保证，50%的数据满足n<=1000，100%的数据满足n<=1000000且P<=1000000000。
输出格式：
仅含一个整数，表示不同的长度为2n的有趣的数列个数mod P的值。
输入输出样例
3 10
5
对应的5个有趣的数列分别为（1,2,3,4,5,6），（1,2,3,5,4,6），（1,3,2,4,5,6），（1,3,2,5,4,6），（1,4,2,5,3,6）。
卡特兰数。
可以在网格上模拟（以奇偶为横纵坐标），发现是卡特兰数的推导。
有模数，不需带高精。
现考虑怎样推：
卡特兰数四个公式：
1.h(n)=h(0)*h(n-1)+h(1)*h(n-2)+.....h(n-1)*h(0).
2.h(n)=C(2n,n)-C(2n,n-1).
3.h(n)=C(2n,n)/(n+1)
4.h(n)=h(n-1)*(4n+2)/(n+2)
现考虑用哪个？？
1肯定不行，2会炸时间，4模意义下的p不一定是质数，不能用逆元，也不可。
3是正解。
一通推导。（第一次算错了，懵了半天）
有趣的来了，接下来通过质因数分解与快速幂来求解。
通过题解得到至此的思路，接下来的靠阿白了（%%%）。
我一开始想的是每个数都质因数分解，分别记录，最后在统计。这样既费时间，也费空间。
正解是（来自阿白的悉心指导）：
直接开一个数组表示指数，分母为减，分子为加（初状态）， 再通过由大向小递推，统计指数。
质因数分解需要用到线性筛，会计算到mindiv（最小质因子），再次加以利用。
妙啊~~~~~~
1 #include<bits/stdc++.h> 2 #define ll long long 3 using namespace std; 4 const int maxn=2000010; 5 int pri[maxn],mindiv[maxn],num; 6 int k; 7 int n,p; 8 void xxs(){ 9 for(int i=2;i<=2*n;i++){ 10 if(!mindiv[i]) mindiv[i]=pri[++num]=i; 11 for(int j=1;j<=num&&((k=i*pri[j])<=2*n)&&pri[j]<=mindiv[i];j++) 12 mindiv[k]=pri[j]; 13 } 14 } 15 int zs[maxn]; 16 ll ans=1; 17 ll fast_pow(int a,int b){ 18 ll s=1; 19 while(b){ 20 if(b&1) s=(s*a)%p; 21 a=(a*a)%p; 22 b>>=1; 23 } 24 return s%p; 25 } 26 int main(){ 27 scanf("%d%d",&n,&p); 28 xxs(); 29 for(int i=1;i<=n;i++) zs[i]=-1; 30 for(int i=n+2;i<=2*n;i++) zs[i]=1; 31 for(int i=2*n;i>=2;i--) if(mindiv[i]!=i)zs[mindiv[i]]+=zs[i],zs[i/mindiv[i]]+=zs[i];//划重点☆☆☆ 32 for(int i=2;i<=2*n;i++) if(mindiv[i]==i)ans=(ans*fast_pow(i,zs[i])%p)%p; 33 printf("%lld",ans); 34 // for(int i=2;i<=n;i++) cout<<i<<" "<<mindiv[i]<<endl; 35 }