卡特兰数
卡特兰数又称卡塔兰数,英文名Catalan number,是组合数学中一个常出现在各种计数问题中出现的数列。以比利时的数学家欧仁·查理·卡塔兰 (1814–1894)的名字来命名,其前几项为(从第零项开始) : 1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 208012, 742900, 2674440, 9694845, 35357670, 129644790, 477638700, 1767263190, 6564120420, 24466267020, 91482563640, 343059613650, 1289904147324, 4861946401452, …
卡特兰数Cn满足以下递推关系 :
例题:P1722 矩阵 II
原题地址
题目描述
众所周知,在中国古代算筹中,红为正,黑为负……
给定一个1*(2n)的矩阵,现让你自由地放入红色算筹和黑色算筹,使矩阵平衡[即对于所有的i(1<=i<=2n),使第1~i格中红色算筹个数大于等于黑色算筹]
问有多少种方案满足矩阵平衡。
见样例解释。
输入输出格式
输入格式:
正整数 n
输出格式:
方案数t对100取模
输入输出样例
输入样例#1:
2
输出样例#1:
2
说明
样例解释: 红 黑 红 黑
红 红 黑 黑
1<=n<=100
代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<string>
using namespace std;
int ka[110];
void tnt()
{
int i,j;
ka[0]=1;
ka[1]=1;
ka[2]=2;
for(i=3;i<=100;i++)
{
for(j=0;j<=i-1;j++)
ka[i]+=ka[j]*ka[i-1-j];
ka[i]%=100;
}
}
int main()
{
int n;
memset(ka,0,sizeof(ka));
tnt();
scanf("%d",&n);
printf("%d\n",ka[n]);
return 0;
}
例题:P1976 鸡蛋饼
原题地址
题目描述
最近小 x 又发现了一个关于圆的有趣的问题:在圆上有2N 个不同的点,小 x 想用 N 条线段把这些点连接起来(每个点只能连一条线段), 使所有的线段都不相交,他想知道这样的连接方案有多少种?
输入输出格式
输入格式:
有且仅有一个正整数 N 。 (N≤2999)
输出格式:
要求的方案数(结果 mod100000007)。
输入输出样例
输入样例#1:
24
输出样例#1:
4057031
代码:
那么假如我们在圆上画了2n个点,顺时针编号为1,2,3,4……,你便会发现,如果一个奇数点和另一个奇数点相连,一定会造成将剩下所没有连线的点分在两边的都是奇数个,而后两边必定有一个点没线连或穿越其中的一条线,那么就不可能完成了。那么我们可以把奇数点看成左括号,偶数点看成右括号,然后把圆切开,就变成了一个括号匹配的方案数问题。为什么圆可以切开呢?因为A连B和B连A是同一种切法。
在这里可以用h(n)=C(2n,n)/(n+1),用扩展欧几里德求逆元。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
using namespace std;
long long int n;
long long int qq[3030];
void ka(long long int x)
{
int i;
for(i=0;i<x;i++)
{
qq[x]=qq[x]+qq[i]*qq[x-1-i];
qq[x]%=100000007;
}
return;
}
int main()
{
long long int i,j;
scanf("%lld",&n);
memset(qq,0,sizeof(qq));
qq[0]=qq[1]=1;
for(i=2;i<=n;i++)
{
ka(i);
}
printf("%lld\n",qq[n]);
return 0;
}
汉诺塔公式
例题
ans=2^n-1
快速幂
思路见例题
例题:P1965 转圈游戏
原题地址
题目描述
n个小伙伴(编号从 0到 n-1)围坐一圈玩游戏。按照顺时针方向给 n个位置编号,从00 到 n-1。最初,第 0号小伙伴在第0号位置,第 1号小伙伴在第 1号位置,……,依此类推。游戏规则如下:每一轮第 0号位置上的小伙伴顺时针走到第 mm 号位置,第 1号位置小伙伴走到第 m+1号位置,……,依此类推,第n − m号位置上的小伙伴走到第 0 号位置,第n∼m+1 号位置上的小伙伴走到第 1 号位置,……,第n−1 号位置上的小伙伴顺时针走到第m-1号位置。
现在,一共进行了 10^k轮,请问 x 号小伙伴最后走到了第几号位置。
输入输出格式
输入格式:
共 1行,包含 4个整数 n,m,k,x每两个整数之间用一个空格隔开。
输出格式:
1个整数,表示 10^k
轮后 x号小伙伴所在的位置编号。
输入输出样例
输入样例#1:
10 3 4 5
输出样例#1:
5
说明
对于 30%的数据,0 < k < 7;
对于 80%的数据,0 < k < 10^7;
对于 100%的数据,1 <n < 1,000,000,0 < m < n,1 ≤ x ≤ n,0 < k < 10^9
代码:
这道题不难可以推出(x+m*10^k)%n
根据模运算的的分配率可以得到 (x%n+m%n*10k%n)%n只需用快速幂求出10k就可以了。
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<cstring>
using namespace std;
int n,m,k,x;
int ksm(int a,int b)//快速幂部分
{
int re=1;
int t=a;
while(b)
{
if(b&1) re=re*t%n;
t=t*t%n;
b>>=1;
}
return re;
}
int main()
{
cin>>n>>m>>k>>x;
cout<<(x%n+m%n*ksm(10,k)%n)%n<<endl;
return 0;
}
斯特林数
详细讲解
第一类斯特林数 S1(n,m) 表示的是将 n 个不同元素构成 m 个圆排列的数目。
第二类斯特林数 S2(n,m) 表示的是把 n 个不同元素划分到 m 个集合的方案数。
例题:P1287 盒子与球
原题地址
题目描述
现有r个互不相同的盒子和n个互不相同的球,要将这n个球放入r个盒子中,且不允许有空盒子。问有多少种方法?
例如:有2个不同的盒子(分别编为1号和2号)和3个不同的球(分别编为1、2、3号),则有6种不同的方法:
输入输出格式
输入格式:
两个整数,n和r,中间用空格分隔。(0≤n, r≤10)
输出格式:
仅一行,一个整数(保证在长整型范围内)。表示n个球放入r个盒子的方法。
输入输出样例
输入样例#1:
3 2
输出样例#1:
6
代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<iostream>
#include<cmath>
using namespace std;
int qq[15][15];
int main()
{
int n,r,ans,i,j;
scanf("%d%d",&n,&r);
memset(qq,0,sizeof(qq));
for(i=1;i<=n;i++)
qq[i][1]=1;
for(i=2;i<=n;i++)
for(j=2;j<=r;j++)
qq[i][j]=qq[i-1][j-1]+j*qq[i-1][j];
ans=qq[n][r];
for(i=2;i<=r;i++)
ans*=i;
printf("%d\n",ans);
return 0;
}
乘法逆元
简单来说就是这样:在(mod p) 意义下( p 是素数),如果 a*a’=1 ,那么我们就说 a’ 是 a 的逆元。当然啦,反过来, a 也是 a’ 的逆元。