POJ1430 Binary Stirling Numbers

@(POJ)[Stirling數, 排列組合, 數形結合]

Description

The Stirling number of the second kind S(n, m) stands for the number of ways to partition a set of n things into m nonempty subsets. For example, there are seven ways to split a four-element set into two parts:
{1, 2, 3} U {4}, {1, 2, 4} U {3}, {1, 3, 4} U {2}, {2, 3, 4} U {1}
{1, 2} U {3, 4}, {1, 3} U {2, 4}, {1, 4} U {2, 3}.
There is a recurrence which allows to compute S(n, m) for all m and n.
S(0, 0) = 1; S(n, 0) = 0 for n > 0; S(0, m) = 0 for m > 0;
S(n, m) = m S(n - 1, m) + S(n - 1, m - 1), for n, m > 0.
Your task is much "easier". Given integers n and m satisfying 1 <= m <= n, compute the parity of S(n, m), i.e. S(n, m) mod 2.
Example:
S(4, 2) mod 2 = 1.
Task
Write a program which for each data set:
reads two positive integers n and m,
computes S(n, m) mod 2,
writes the result.

Input

The first line of the input contains exactly one positive integer d equal to the number of data sets, 1 <= d <= 200. The data sets follow.
Line i + 1 contains the i-th data set - exactly two integers ni and mi separated by a single space, 1 <= mi <= ni <= 10^9.

Output

The output should consist of exactly d lines, one line for each data set. Line i, 1 <= i <= d, should contain 0 or 1, the value of S(ni, mi) mod 2.

Sample Input

1
4 2

Sample Output

1

Solution

題意:
求斯特林數$$ \left{ \begin{array}{} n \ k \end{array}{} \right} % 2$$$$n, m \in [1, 10^9]$$
這題直接求解肯定是會T的, 因此考慮優化.

轉載自sdchr博客
侵刪

代碼附上:

#include<cstdio>
#include<cctype>
using namespace std;

inline int read()
{
	int x = 0, flag = 1;
	char c;
	while(! isdigit(c = getchar()))
		if(c == '-')
			flag *= - 1;
	while(isdigit(c))
		x = x * 10 + c - '0', c = getchar();
	return x * flag;
}

void println(int x)
{
	if(x < 0)
		putchar('-'), x *= - 1;
	if(x == 0)
		putchar('0');
	int ans[1 << 5], top = 0;
	while(x)
		ans[top ++] = x % 10, x /= 10;
	for(; top; top --)
		putchar(ans[top - 1] + '0');
	putchar('\n');
}

long long getQuantity(int x)
{
	long long ret = 0;
	
	for(int i = 2; i <= x; i <<= 1)
		ret += x / i;
		
	return ret;
}

int calculate(int x, int y)
{
	return getQuantity(x) - getQuantity(y) - getQuantity(x - y) == 0;
}

int main()
{
	int T = read();
	
	while(T --)
	{
		int n = read(), m = read();
		int d = n - m, oddQua = (m + 1) / 2;
		println(calculate(d + oddQua - 1, oddQua - 1));
	}
}