P1363 火车进出栈问题
时间: 1000ms / 空间: 131072KiB / Java类名: Main
描述
一列火车n节车厢,依次编号为1,2,3,…,n。每节车厢有两种运动方式,进栈与出栈,问n节车厢出栈的可能排列方式有多少种。
输入格式
一个数,n(n<=30000)
输出格式
一个数s表示n节车厢出栈的可能排列方式
测试样例1
输入
3
输出
5
思路:首先知道catalan数通项公式:h(n)=c(2n,n)/(n+1);
说白了就是求c(2n,n)(这个公式再不知道我也是没辙),
将每个数分解质因数,n+1 to 2*n 的每个质因数个数加上(因为要除(n+1),所以从n+2开始),1 to n-m的每个质因数个数减去,数组记录,最后高精乘以单精,轻松AC Σ(っ °Д °;)っ
program df;
var i,j,n,m,x,y,z,k,t:longint;
a,d:array[0..100000] of longint;
procedure deal(n,m:longint);
var i,j,x,y:longint;
begin
for i:=m+2 to n do
begin
x:=i;
for j:=2 to trunc(sqrt(x)) do
while x mod j=0 do
begin
inc(d[j]);
x:=x div j;
end;
if x>1 then inc(d[x]);
end;
for i:=1 to n-m do
begin
x:=i;
for j:=2 to trunc(sqrt(x)) do
while x mod j=0 do
begin
dec(d[j]);
x:=x div j;
end;
if x>1 then dec(d[x]);
end;
end;
begin
readln(n);
deal(2*n,n);
a[1]:=1; t:=1;
for i:=2 to 3*n do
while d[i]>0 do
begin
for j:=1 to t do
a[j]:=a[j]*i;
for x:=1 to t do
begin
a[x+1]:=a[x+1]+a[x] div 10;
a[x]:=a[x] mod 10;
end;
j:=t+1;
while a[j]<>0 do
begin
a[j+1]:=a[j+1]+a[j] div 10;
a[j]:=a[j] mod 10;
inc(j);
end;
t:=j-1;
dec(d[i]);
end;
for i:=t downto 1 do
write(a[i]);
end.