Drolca

Apologize To Drolca
随笔 - 28, 文章 - 1, 评论 - 6, 引用 - 0
数据加载中……

2012年8月22日

Topcoder SRM550 div1 500


Sierpinski三角形
      

杨辉三角,二项式系数奇偶性的判定 C(k,n)| k&n==k , 或者比较n!与{k!,(n-k)!}中2的个数

posted @ 2012-08-22 10:05 Drolca 阅读(64) | 评论 (0)编辑 收藏

2012年4月13日

poj 3691 AC自动机+DP

     摘要: Code highlighting produced by Actipro CodeHighlighter (freeware)http://www.CodeHighlighter.com/-->#include <iostream>#include <math.h>using namespace std;const int...  阅读全文

posted @ 2012-04-13 19:28 Drolca 阅读(121) | 评论 (0)编辑 收藏

hdu 2222 AC自动机

#include <iostream>
using namespace std;

const int MAXN=50*10005;
const int MAXL=1005;
const int K=26;

struct Node
{
    Node 
*next[K], *fail;
    
int flag, id;
    
void Init(int index)
    
{
        id
=index;
        flag
=0;
        fail
=NULL;
        
for(int i=0; i<K; i++)next[i]=NULL;
    }

}
* Q[MAXN/2], *root, T[MAXN];// Q for queue, root&T for tree;

int index=0;
Node 
* newNode()
{
    T[index].Init(index);
    
return &T[index++];
}


int tokind(char k){return k-'a';}

void insert(char *str)
{    
    
if(root==NULL)
        root
=newNode();

    Node 
*now=root;
    
for(int i=0; str[i]; i++)
    
{
        
int kind=tokind(str[i]);
        
if(now->next[kind]==NULL)
            now
->next[kind]=newNode();
        now
=now->next[kind];
    }

    now
->flag++;
}


void buildAC()
{
    
int head=0, tail=0;
    root
->fail=NULL;

    Q[tail
++]=root;
    
while(head<tail)
    
{
        Node 
*now=Q[head++];
        
for(int i=0; i<K; i++)
        
{
            
if(now->next[i]!=NULL)
            
{
                
if(now==root)now->next[i]->fail=root;
                
else
                
{
                    Node 
*p=now->fail;
                    
while(p->next[i]==NULL&&p!=root)p=p->fail;
                    p
=p->next[i];
                    now
->next[i]->fail=(p==NULL)?root:p;
                }

                Q[tail
++]=now->next[i];
            }

        }

    }

}


int query(char *str)
{
    
int res=0;
    Node 
*p=root;
    
for(int i=0; str[i]; i++)
    
{
        
int kind=tokind(str[i]);
        
while(p->next[kind]==NULL&&p!=root)p=p->fail;
        p
=p->next[kind];
        p
=(p==NULL)?root:p;
        Node 
*now=p;
        
while(now!=root&&now->flag!=-1)
        
{
            res
+=now->flag;
            now
->flag=-1;
            now
=now->fail;
        }

    }

    
return res;
}


int main()
{
    
//freopen ("in.txt", "r", stdin);
    int T;
    scanf(
"%d"&T);
    
while(T--)
    
{
        
int n;
        root
=NULL;
        scanf(
"%d",&n);getchar();        
        
while(n--)
        
{
            
char str[55];
            scanf(
"%s",str);
            insert(str);
        }

        buildAC();
        
char ch[1000005];
        scanf(
"%s",ch);
        printf(
"%d\n",query(ch));
    }

    
return 0;
}

posted @ 2012-04-13 15:17 Drolca 阅读(114) | 评论 (0)编辑 收藏

2012年4月9日

Topcoder SRM 538 1050P

 1#include <iostream>
 2#include <vector>
 3#include <string>
 4using namespace std;
 5
 6#define madd(a,b) a=(a+b)%MOD
 7
 8const int MOD=1000000007;
 9const int MAXC=55, MAXH=75, MAXW=10, MAXB=15;
10long g[MAXC], Comb[MAXC][MAXC];
11int dp[MAXW][MAXB][MAXC][MAXH][MAXW][2];
12
13class SkewedPerspective
14{
15
16public:
17    int countThem(vector <int> cubes, int B, int w)
18    {
19        int n=cubes.size(), total=0;
20        int i, j, k;
21        for(i=0; i<n; i++) total+=cubes[i];
22        
23        for(i=0; i<=total; i++)
24            for(j=0; j<=i; j++)
25                Comb[i][j]=(j?(Comb[i-1][j]+Comb[i-1][j-1])%MOD:1);
26        g[0]=1;
27        for(i=0; i<n; i++)
28            for(j=total; j; j--)
29                for(k=1; k<=&& k<=cubes[i]; k++)
30                    g[j]=(g[j]+g[j-k]*Comb[j][k])%MOD;
31
32        long ans=0;
33        dp[1][0][0][0][0][0]=1;
34        for(int tower=1; tower<=w; tower++)for(int black=0; black<=B; black++)
35        for(int color=0; color<=total; color++)for(int need=0; need<=total+black*2; need++)
36        for(int needOdd=0; needOdd<=tower; needOdd++)for(int lastBlack=0; lastBlack<2; lastBlack++)
37        {
38            int x=dp[tower][black][color][need][needOdd][lastBlack];
39            if(!x)continue;
40            //get result
41            if(black+color>0 && (B-black)*2+total-color>=need && total-color>=needOdd)
42                ans=(ans+g[color]*x)%MOD;
43            //put colored
44            madd(dp[tower][black][color+1][need][needOdd][0], x);
45            if(lastBlack) continue;
46            //put black
47            int layer=black*2+color-(tower-1);
48            for(int blackSize=1; blackSize+black*2<=B*2; blackSize++)
49                if(blackSize%2==0)
50                    madd(dp[tower][black+blackSize/2][color][need][needOdd][1], x);
51                else
52                {
53                    if(!layer && blackSize==1continue;  //"b1bbbb"
54                    int needNow=(layer?layer-1:layer+1);
55                    if(need+needNow<=total+B*2)
56                        madd(dp[tower+1][black+(blackSize+1)/2][color][need+needNow][needOdd+needNow%2][1], x);
57                }

58
59        }

60        return int(ans);
61    }

62}
;
63
64int main()
65{
66    SkewedPerspective a;
67    int cubes[]={110};
68    vector<int> t(cubes, cubes+3);
69
70    cout<<a.countThem(t, 12)<<endl;
71    return 0;
72}

73

posted @ 2012-04-09 17:30 Drolca 阅读(100) | 评论 (0)编辑 收藏

2011年5月21日

分享一篇好文章《主题:说说字符集和编码》

转载自:http://www.iteye.com/topic/398782
很久很久以前,有一群人,他们决定用8个可以开合的晶体管来组合成不同的状态,以表示世界上的万物。他们看到8个开关状态是好的,于是他们把这称为"字节"。
 

再后来,他们又做了一些可以处理这些字节的机器,机器开动了,可以用字节来组合出很多状态,状态开始变来变去。他们看到这样是好的,于是它们就这机器称为"计算机"。 



开始计算机只在美国用。八位的字节一共可以组合出256(2的8次方)种不同的状态。 

他们把其中的编号从0开始的32种状态分别规定了特殊的用途,一但终端、打印机遇上约定好的这些字节被传过来时,就要做一些约定的动作。遇上00x10, 终端就换行,遇上0x07, 终端就向人们嘟嘟叫,例好遇上0x1b, 打印机就打印反白的字,或者终端就用彩色显示字母。他们看到这样很好,于是就把这些0x20以下的字节状态称为"控制码"。  

他们又把所有的空格、标点符号、数字、大小写字母分别用连续的字节状态表示,一直编到了第127号,这样计算机就可以用不同字节来存储英语的文字了。大家看到这样,都感觉很好,于是大家都把这个方案叫做 ANSI 的"Ascii"编码(American Standard Code for Information Interchange,美国信息互换标准代码)。当时世界上所有的计算机都用同样的ASCII方案来保存英文文字。 

后来,就像建造巴比伦塔一样,世界各地的都开始使用计算机,但是很多国家用的不是英文,他们的字母里有许多是ASCII里没有的,为了可以在计算机保存他们的文字,他们决定采用127号之后的空位来表示这些新的字母、符号,还加入了很多画表格时需要用下到的横线、竖线、交叉等形状,一直把序号编到了最后一个状态255。从128到255这一页的字符集被称"扩展字符集"。从此之后,贪婪的人类再没有新的状态可以用了,美帝国主义可能没有想到还有第三世界国家的人们也希望可以用到计算机吧!  

等中国人们得到计算机时,已经没有可以利用的字节状态来表示汉字,况且有6000多个常用汉字需要保存呢。但是这难不倒智慧的中国人民,我们不客气地把那些127号之后的奇异符号们直接取消掉, 规定:一个小于127的字符的意义与原来相同,但两个大于127的字符连在一起时,就表示一个汉字,前面的一个字节(他称之为高字节)从0xA1用到0xF7,后面一个字节(低字节)从0xA1到0xFE,这样我们就可以组合出大约7000多个简体汉字了。在这些编码里,我们还把数学符号、罗马希腊的字母、日文的假名们都编进去了,连在 ASCII 里本来就有的数字、标点、字母都统统重新编了两个字节长的编码,这就是常说的"全角"字符,而原来在127号以下的那些就叫"半角"字符了。  

中国人民看到这样很不错,于是就把这种汉字方案叫做 "GB2312"。GB2312 是对 ASCII 的中文扩展。 

但是中国的汉字太多了,我们很快就就发现有许多人的人名没有办法在这里打出来,特别是某些很会麻烦别人的国家领导人。于是我们不得不继续把 GB2312 没有用到的码位找出来老实不客气地用上。 

后来还是不够用,于是干脆不再要求低字节一定是127号之后的内码,只要第一个字节是大于127就固定表示这是一个汉字的开始,不管后面跟的是不是扩展字符集里的内容。结果扩展之后的编码方案被称为 GBK 标准,GBK 包括了 GB2312 的所有内容,同时又增加了近20000个新的汉字(包括繁体字)和符号。  

后来少数民族也要用电脑了,于是我们再扩展,又加了几千个新的少数民族的字,GBK 扩成了 GB18030。从此之后,中华民族的文化就可以在计算机时代中传承了。 

中国的程序员们看到这一系列汉字编码的标准是好的,于是通称他们叫做 "DBCS"(Double Byte Charecter Set 双字节字符集)。在DBCS系列标准里,最大的特点是两字节长的汉字字符和一字节长的英文字符并存于同一套编码方案里,因此他们写的程序为了支持中文处理,必须要注意字串里的每一个字节的值,如果这个值是大于127的,那么就认为一个双字节字符集里的字符出现了。那时候凡是受过加持,会编程的计算机僧侣们都要每天念下面这个咒语数百遍:  

"一个汉字算两个英文字符!一个汉字算两个英文字符......" 



因为当时各个国家都像中国这样搞出一套自己的编码标准,结果互相之间谁也不懂谁的编码,谁也不支持别人的编码,连大陆和台湾这样只相隔了150海里,使用着同一种语言的兄弟地区,也分别采用了不同的 DBCS 编码方案。当时的中国人想让电脑显示汉字,就必须装上一个"汉字系统",专门用来处理汉字的显示、输入的问题,但是那个台湾的愚昧封建人士写的算命程序就必须加装另一套支持 BIG5 编码的什么"倚天汉字系统"才可以用,装错了字符系统,显示就会乱了套!这怎么办?而且世界民族之林中还有那些一时用不上电脑的穷苦人民,他们的文字又怎么办? 

真是计算机的巴比伦塔命题啊! 

正在这时,大天使加百列及时出现了:一个叫 ISO (国际标谁化组织)的国际组织决定着手解决这个问题。他们采用的方法很简单:废了所有的地区性编码方案,重新搞一个包括了地球上所有文化、所有字母和符号的编码!他们打算叫它"Universal Multiple-Octet Coded Character Set",简称 UCS, 俗称 "UNICODE"。 

UNICODE 开始制订时,计算机的存储器容量极大地发展了,空间再也不成为问题了。于是 ISO 就直接规定必须用两个字节,也就是16位来统一表示所有的字符,对于ascii里的那些"半角"字符,UNICODE 包持其原编码不变,只是将其长度由原来的8位扩展为16位,而其他文化和语言的字符则全部重新统一编码。由于"半角"英文符号只需要用到低8位,所以其高8位永远是0,因此这种大气的方案在保存英文文本时会多浪费一倍的空间。  

这时候,从旧社会里走过来的程序员开始发现一个奇怪的现象:他们的strlen函数靠不住了,一个汉字不再是相当于两个字符了,而是一个!是的,从 UNICODE 开始,无论是半角的英文字母,还是全角的汉字,它们都是统一的"一个字符"!同时,也都是统一的"两个字节",请注意"字符"和"字节"两个术语的不同,"字节"是一个8位的物理存贮单元,而"字符"则是一个文化相关的符号。在UNICODE 中,一个字符就是两个字节。一个汉字算两个英文字符的时代已经快过去了。 

从前多种字符集存在时,那些做多语言软件的公司遇上过很大麻烦,他们为了在不同的国家销售同一套软件,就不得不在区域化软件时也加持那个双字节字符集咒语,不仅要处处小心不要搞错,还要把软件中的文字在不同的字符集中转来转去。UNICODE 对于他们来说是一个很好的一揽子解决方案,于是从 Windows NT 开始,MS 趁机把它们的操作系统改了一遍,把所有的核心代码都改成了用 UNICODE 方式工作的版本,从这时开始,WINDOWS 系统终于无需要加装各种本土语言系统,就可以显示全世界上所有文化的字符了。  

但是,UNICODE 在制订时没有考虑与任何一种现有的编码方案保持兼容,这使得 GBK 与UNICODE 在汉字的内码编排上完全是不一样的,没有一种简单的算术方法可以把文本内容从UNICODE编码和另一种编码进行转换,这种转换必须通过查表来进行。 

如前所述,UNICODE 是用两个字节来表示为一个字符,他总共可以组合出65535不同的字符,这大概已经可以覆盖世界上所有文化的符号。如果还不够也没有关系,ISO已经准备了UCS-4方案,说简单了就是四个字节来表示一个字符,这样我们就可以组合出21亿个不同的字符出来(最高位有其他用途),这大概可以用到银河联邦成立那一天吧!  



UNICODE 来到时,一起到来的还有计算机网络的兴起,UNICODE 如何在网络上传输也是一个必须考虑的问题,于是面向传输的众多 UTF(UCS Transfer Format)标准出现了,顾名思义,UTF8就是每次8个位传输数据,而UTF16就是每次16个位,只不过为了传输时的可靠性,从UNICODE到UTF时并不是直接的对应,而是要过一些算法和规则来转换。 

受到过网络编程加持的计算机僧侣们都知道,在网络里传递信息时有一个很重要的问题,就是对于数据高低位的解读方式,一些计算机是采用低位先发送的方法,例如我们PC机采用的 INTEL 架构,而另一些是采用高位先发送的方式,在网络中交换数据时,为了核对双方对于高低位的认识是否是一致的,采用了一种很简便的方法,就是在文本流的开始时向对方发送一个标志符。如果之后的文本是高位在位,那就发送"FEFF",反之,则发送"FFFE"。不信你可以用二进制方式打开一个UTF-X格式的文件,看看开头两个字节是不是这两个字节?  



讲到这里,我们再顺便说说一个很著名的奇怪现象:当你在 windows 的记事本里新建一个文件,输入"联通"两个字之后,保存,关闭,然后再次打开,你会发现这两个字已经消失了,代之的是几个乱码!呵呵,有人说这就是联通之所以拼不过移动的原因。 

其实这是因为GB2312编码与UTF8编码产生了编码冲撞的原因。 

从网上引来一段从UNICODE到UTF8的转换规则: 

Unicode 

UTF-8  
0000 - 007F 

0xxxxxxx 



0080 - 07FF 

110xxxxx 10xxxxxx 



0800 - FFFF 

1110xxxx 10xxxxxx 10xxxxxx 



例如"汉"字的Unicode编码是6C49。6C49在0800-FFFF之间,所以要用3字节模板:1110xxxx 10xxxxxx 10xxxxxx。将6C49写成二进制是:0110 1100 0100 1001,将这个比特流按三字节模板的分段方法分为0110 110001 001001,依次代替模板中的x,得到:1110-0110 10-110001 10-001001,即E6 B1 89,这就是其UTF8的编码。  

而当你新建一个文本文件时,记事本的编码默认是ANSI, 如果你在ANSI的编码输入汉字,那么他实际就是GB系列的编码方式,在这种编码下,"联通"的内码是: 

c1 1100 0001 

aa 1010 1010 

cd 1100 1101 

a8 1010 1000 

注意到了吗?第一二个字节、第三四个字节的起始部分的都是"110"和"10",正好与UTF8规则里的两字节模板是一致的,于是再次打开记事本时,记事本就误认为这是一个UTF8编码的文件,让我们把第一个字节的110和第二个字节的10去掉,我们就得到了"00001 101010",再把各位对齐,补上前导的0,就得到了"0000 0000 0110 1010",不好意思,这是UNICODE的006A,也就是小写的字母"j",而之后的两字节用UTF8解码之后是0368,这个字符什么也不是。这就是只有"联通"两个字的文件没有办法在记事本里正常显示的原因。  

而如果你在"联通"之后多输入几个字,其他的字的编码不见得又恰好是110和10开始的字节,这样再次打开时,记事本就不会坚持这是一个utf8编码的文件,而会用ANSI的方式解读之,这时乱码又不出现了。 

posted @ 2011-05-21 21:36 Drolca 阅读(84) | 评论 (1)编辑 收藏

2010年1月1日

hdu 2102

#include <iostream>
using namespace std;
const int M=10;
char map[2][M][M];
bool vis[2][M][M];
int dx[]={1,0,-1,0};
int dy[]={0,1,0,-1};

struct point
{
    
int layer;
    
int x,y;
    
int time;
}
Q[200];

bool BFS(int m,int n,int t)
{
    point now,next;
    now.layer
=now.x=now.y=now.time=0;
    
int Front=0;
    
int Near=1;
    Q[Front]
=now;
    vis[now.layer][now.x][now.y]
=true;
    
while(Front<Near)
    
{
        now
=Q[Front++];
        
if(map[now.layer][now.x][now.y]=='P')
        
{
            
if(now.time<=t)
                
return true;
            
return false;
        }

        
if(map[now.layer][now.x][now.y]=='#')
            now.layer
=!now.layer;

        
if(map[now.layer][now.x][now.y]=='P')
        
{
            
if(now.time<=t)
                
return true;
            
return false;
        }

        
if(map[now.layer][now.x][now.y]=='*'||map[now.layer][now.x][now.y]=='#')
            
continue;
        
int k;
        
for(k=0;k<4;k++)
        
{
            next.layer
=now.layer;
            next.time
=now.time+1;
            next.x
=now.x+dx[k];
            next.y
=now.y+dy[k];
            
if(!vis[next.layer][next.x][next.y]&&next.x>=0&&next.x<m&&next.y>=0&&next.y<n)
            
{
                Q[Near
++]=next;
                vis[next.layer][next.x][next.y]
=true;
            }

        }

        
    }

    
return false;
}

int main()
{
    
int T;
    scanf(
"%d",&T);
    
while(T--)
    
{
        
int m,n,t;
        scanf(
"%d%d%d",&m,&n,&t);
        
int i,j,k;
        
for(k=0;k<2;k++)
        
{
            
for(i=0;i<m;i++)
            
{
                
char s[M];
                scanf(
"%s",&s);
                
for(j=0;j<n;j++)
                
{
                    map[k][i][j]
=s[j];
                    vis[k][i][j]
=false;
                }

            }

        }

        
if(BFS(m,n,t))
            printf(
"YES\n");
        
else 
            printf(
"NO\n");
    }

    system(
"pause");
    
return 0;
}

posted @ 2010-01-01 22:21 Drolca 阅读(212) | 评论 (0)编辑 收藏

9*9数独游戏

#include <iostream>
using namespace std;
const int M=10;
bool userow[M][M],usecol[M][M],useblock[M][M];
int map[M][M];

struct node{
    
int x,y;
    
int num;
}
sudu[M*M];

int find(int x,int y)
{
    
int row=x/3;
    
int col=y/3;
    
return 3*row+col;
}


bool dfs(int n,int cnt)
{
    
if(n==cnt)return 1;
    
int i;
    
for(i=1;i<M;i++)
    
{
        
if(!userow[sudu[n].x][i]&&!usecol[sudu[n].y][i]&&!useblock[find(sudu[n].x,sudu[n].y)][i])
        
{
            userow[sudu[n].x][i]
=true;
            usecol[sudu[n].y][i]
=true;
            useblock[find(sudu[n].x,sudu[n].y)][i]
=true;
            sudu[n].num
=i;
            
if(dfs(n+1,cnt))
                
return 1;
            userow[sudu[n].x][i]
=false;
            usecol[sudu[n].y][i]
=false;
            useblock[find(sudu[n].x,sudu[n].y)][i]
=false;
            sudu[n].num
=0;

        }

    }

    
return 0;

}


int main()
{
    
int T;
    scanf(
"%d",&T);
    
while(T--)
    
{
        memset(userow,
false,sizeof(userow));
        memset(usecol,
false,sizeof(usecol));
        memset(useblock,
false,sizeof(useblock));

        
int i,j;
        
int cnt=0;
        
for(i=0;i<M-1;i++)
        
{
            
char mess[M];
            scanf(
"%s",&mess);
            
for(j=0;j<M-1;j++)
            
{
                map[i][j]
=mess[j]-'0';
                
if(map[i][j])
                
{
                    userow[i][map[i][j]]
=true;
                    usecol[j][map[i][j]]
=true;
                    useblock[find(i,j)][map[i][j]]
=true;
                }

                
else
                
{
                    sudu[cnt].x
=i;
                    sudu[cnt].y
=j;
                    sudu[cnt].num
=0;
                    cnt
++;
                }

                
            }


        }


        dfs(
0,cnt);

        
for(i=0;i<cnt;i++)
            map[sudu[i].x][sudu[i].y]
=sudu[i].num;
        
        
for(i=0;i<M-1;i++){
            
for(j=0;j<M-1;j++)
                printf(
"%d",map[i][j]);
            printf(
"\n");
        }

    }


    
return 0;
}

posted @ 2010-01-01 21:19 Drolca 阅读(286) | 评论 (1)编辑 收藏

2009年11月21日

有一个悲剧...最小点割

     摘要:   #include <iostream>using namespace std;const int maxn=200;const int INF=1000000;int g[maxn][maxn];int f[maxn][maxn];int r[maxn][maxn];in...  阅读全文

posted @ 2009-11-21 21:08 Drolca 阅读(66) | 评论 (0)编辑 收藏

2009年10月5日

pku 3744 Scout YYF I

#include <iostream>
#include 
<algorithm>
using namespace std;

const int maxn=12;
int n;
double p;
struct matrix
{
    
double m[2][2];
}
;
int mine[maxn];
matrix 
operator*(const matrix&a,const matrix&b)
{
    matrix tmp;
    
int i,j,k;
    
for(i=0;i<2;i++)
        
for(j=0;j<2;j++)
        
{
            tmp.m[i][j]
=0;
            
for(k=0;k<2;k++)
                tmp.m[i][j]
+=a.m[i][k]*b.m[k][j];
        }

    
return tmp;
}


matrix power(
int k)
{
    matrix tmp,res;
    matrix A;
    A.m[
0][0]=p,A.m[0][1]=1-p,A.m[1][0]=1,A.m[1][1]=0;
    matrix B;
    B.m[
0][0]=1,B.m[0][1]=0,B.m[1][0]=0,B.m[1][1]=1;
    
if(k==0)
        
return B;
    
if(k==1)
        
return A;
    
else 
    
{
        tmp
=power(k/2);
        res
=tmp*tmp;
        
if(k%2==1)
            res
=res*A;
        
return res;
    }

}


void slove(int n,double p)
{
    
double a=1,b=0;
    
int i;
    
for(i=0;i<n;i++)
        scanf(
"%d",&mine[i]);
    sort(mine,mine
+n);

    
double f2=1.0,f1=0.0;
    
int now=1;
    
for(i=0;i<n;i++)
    
{
        
if((mine[i]-1)-now>=0)
        
{
            matrix tmp
=power(mine[i]-1-now);
            f2
=(tmp.m[0][0]*f2+tmp.m[0][1]*f1)*(1-p);
            f1
=0;
            now
=mine[i]+1;
        }

        
else
        
{
            printf(
"%.7lf\n",0.0);
            
return;
        }

    }

    printf(
"%.7lf\n",f2);
}


int main()
{
    
while(scanf("%d %lf",&n,&p)!=EOF)
        slove(n,p);
    
return 0;
}


posted @ 2009-10-05 18:48 Drolca 阅读(107) | 评论 (0)编辑 收藏

hdu 2292 Minimum Heap

#include <iostream>
using namespace std;

const int maxn=1005;
__int64 n,m;
__int64 F[maxn];
int c[maxn][maxn];

int cal(int n)
{
    
int t = 1;
    
while (t <= n) t = t * 2 + 1;
    t 
= (t - 1/ 2;
    t 
= (t - 1/ 2;
    
int sum = n - 1 - t;
    
if (sum > 2 * t + 1{
        sum 
= 2 * t + 1;
    }

    
return sum;
}


void calc_c() {
    
for (int i = 0; i < maxn; i++
    
{
        c[i][
0= c[i][i] = 1;
        
for (int j = 1; j < i; j++
        
{
            c[i][j] 
= (c[i - 1][j - 1+ c[i - 1][j]) % m;
        }

    }

}


__int64 slove(
int n)
{
    
if(F[n])
        
return F[n];
    
if(n==0||n==1)
        
return 1;
    
int left=cal(n);
    
int right=(n-1)-left;
    
return F[n]=( (slove(left)*slove(right) )%m )*(__int64)c[n-1][left]%m;
}

int main()
{
    
int T;
    scanf(
"%d",&T);
    
while(T--)
    
{
        memset(F,
0,sizeof(F));
        scanf(
"%I64d%I64d",&n,&m);
        calc_c();
        __int64 ans
=slove(n);
        printf(
"%I64d\n",ans);
    }

    
return 0;
}

posted @ 2009-10-05 11:19 Drolca 阅读(163) | 评论 (1)编辑 收藏