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

bzoj1485

时间:2015-05-09 14:43:43      阅读:129      评论:0      收藏:0      [点我收藏+]

标签:

首先考虑dp,设f[i,j]表示1~i用过了,期中j个放在偶数位然后转移大家都会

这显然TLE,我们观察这个dp,任意前i个数,无论怎么放,放在奇数位的数的个数一定要大于等于放在偶数位的个数

于是很明显这是经典的卡特兰数模型

注意这里涉及到了除法取模,而模数不一定是质数

很显然的想法是分解质因数然后约分

但有更漂亮的做法,http://blog.csdn.net/jasonzhu8/article/details/5949622

技术分享
 1 var v,b,p:array[0..2000010] of longint;
 2     n,i,j,mo,t:longint;
 3     ans:int64;
 4 
 5 procedure calc(x,w:longint);
 6   begin
 7     while v[x]<>0 do
 8     begin
 9       inc(b[v[x]],w);
10       x:=x div v[x];
11     end;
12     inc(b[x],w);
13   end;
14 
15 function quick(x:int64;y:longint):int64;
16   begin
17     quick:=1;
18     while y>0 do
19     begin
20       if y mod 2=1 then quick:=quick*x mod mo;
21       y:=y div 2;
22       x:=x*x mod mo;
23     end;
24   end;
25 
26 begin
27   readln(n,mo);
28   for i:=2 to 2*n do
29   begin
30     if v[i]=0 then
31     begin
32       inc(t);
33       p[t]:=i;
34     end;
35     for j:=1 to t do
36     begin
37       if int64(i)*int64(p[j])>2*n then break;
38       v[i*p[j]]:=p[j];
39       if i mod p[j]=0 then break;
40     end;
41   end;
42   for i:=n+2 to 2*n do
43     calc(i,1);
44   for i:=2 to n do
45     calc(i,-1);
46   ans:=1;
47   for i:=1 to t do
48     ans:=ans*quick(p[i],b[p[i]]) mod mo;
49   writeln(ans);
50 end.
View Code

 

bzoj1485

标签:

原文地址:http://www.cnblogs.com/phile/p/4489857.html

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