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));
    }
}