LightOJ1170 - Counting Perfect BST（卡特兰数）

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<set>
 4 #include<algorithm>
 5 using namespace std;
 6 long long pow(long long x,int y){
 7     long long res=1;
 8     for(int i=0; i<y; ++i) res*=x;
 9     return res;
10 }
11 long long pownum[111111];
12 int pn;
13 int getCnt(long long a,long long b){
14     return (upper_bound(pownum,pownum+pn,b)-pownum)-(lower_bound(pownum,pownum+pn,a)-pownum);
15 }
16 long long d[1111];
17 int main(){
18     for(long long i=2; i*i<=10000000000L; ++i){
19         for(int j=2; ; ++j){
20             if(pow(i,j)>10000000000L) break;
21             pownum[pn++]=pow(i,j);
22         }
23     }
24     sort(pownum,pownum+pn);
25     pn=unique(pownum,pownum+pn)-pownum;
26     d[0]=d[1]=1;
27     for(int i=2; i<1111; ++i){
28         for(int j=1; j<=i; ++j){
29             d[i]+=d[j-1]*d[i-j];
30             d[i]%=100000007;
31         }
32     }
33     d[0]=0;
34     long long a,b;
35     int t;
36     scanf("%d",&t);
37     for(int cse=1; cse<=t; ++cse){
38         scanf("%lld%lld",&a,&b);
39         printf("Case %d: %lld\n",cse,d[getCnt(a,b)]); 
40     }
41     return 0;
42 }