QDU 良辰最喜欢良辰最喜欢对那些自认能力出众的人出手(Lucas定理推广)

良辰最喜欢对那些自认能力出众的人出手(感谢cqupt..

发布时间: 2015年11月8日 22:31   最后更新: 2015年11月8日 22:35   时间限制: 1000ms   内存限制: 128M

描述

果果一直很谦虚，但是良辰还是对他出手了（汗。。

良辰拿出来一个 Pascal 三角形（也叫杨辉三角形。（请不要问是怎么拿出来的。。

我们用一个矩阵来表示 Pascal 三角形

0: C(0, 0)

1: C(1, 0) C(1, 1)

2: C(2, 0) C(2, 1) C(2, 2)

...

其中最左边的数字表示行号。C(x, y) 表示一个组合数，即 x 个元素中取 y 个的方法数。

良辰只有一个问题：Pascal 三角形中的第 n 行有多少个奇数？

这样的问题对于果果来说太水了，于是他把这个问题给了你。

输入

第一行为一个整数 T，表示数据组数。
每组数据只有一行，包含一个整数 n，表示 Pascal 三角形的第 n 行。
T <= 2000,0 <= n <= 2^31

输出

每组数据输出一行，包含一个整数，表示Pascal 三角形的第 n 行中的奇数个数。

样例输入1  复制

2
0
1

样例输出1

1
2

思路：可以将问题分析，求m从0到n的C(n, m)%2为1的个数，而对于一个特定的m，求C(n, m)%2可通过Lucas定理去求。

Lucas定理是将n，m分解成p进制，则C(n, m)%p = C(a[n], b[n])*C(a[n-1], b[n-1])*...*C(a[0], b[0])(mod p)，本题p即是2。

令n=100101(二进制)，则m就是从000000到100101的枚举.

我们知道C(0,1)=0, C(0, 0) = C(1, 0) = C(1, 1) = 1.

因此如果n=100101对应的0的位置在m中对应的值为1，则此时根据Lucas定理C(n, m)=0.所以从0~n的m只有在0的位置填0，1的位置填0或1，才能保证C(n, m)%2 = 1;

所以求得n的二进制中1的个数cnt，然后2^cnt就是答案。

Code:

#include <iostream>
using namespace std;
int main()
{
	int t, n, cnt;
	cin >> t;
	for(;t--;)
	{
		cin >> n;
		cnt = 0;
		while(n)
		{
			if(n&1) ++cnt;
			n >>= 1;
		}
		cout << (1ll<<cnt) << endl;
	}
	return 0;
}

继续加油~