SPOJ NWERC11A Binomial coefficients

            组合数单调性+二分。这道题犯了一个很二逼的错误，调了很久很久才能AC。。。好无语啊，以后要细心些了。看到这题的第一反应是，不会做。。。但是后来仔细想了想，其实我们很容易发现，对于其中的满足的情况C（n， k），我们只需要求k <= n/2的情况就可以了，然后就可以发现，其实k是不超过30的，因为C（60， 30）已经超过了1e15。然后就想着枚举k，然后二分n，求答案了！

            但是要是这样写的话我们发现，即使k 小于30，在相乘的时候也会超long long。怎么办呢？仔细想想，我们可以用二分求出小于6的情况（单独这样求是因为k越小n越大）。然后当k大于6的时候，即使k取最小，7，n只要大于500的时候C（n， k）就已经大于1e15了。这样的话问题就迎刃而解了。

            1、k小于7时二分n。

            2、k大于7时直接暴力枚举。

#include<algorithm>
#include<iostream>
#include<cstring>
#include<fstream>
#include<sstream>
#include<cstdlib>
#include<vector>
#include<string>
#include<cstdio>
#include<bitset>
#include<queue>
#include<stack>
#include<cmath>
#include<map>
#include<set>
#define FF(i, a, b) for(int i=a; i<b; i++)
#define FD(i, a, b) for(int i=a; i>=b; i--)
#define REP(i, n) for(int i=0; i<n; i++)
#define CLR(a, b) memset(a, b, sizeof(a))
#define debug puts("**debug**")
#define LL long long
#define PB push_back
#define MP make_pair
using namespace std;

const int N = 500;
const LL INF = 1e18;

LL c[N][N];

struct Node
{
    LL a, b;
    void pt()
    {
        cout << "(" << a << "," << b << ")";
    }
    bool operator < (const Node& rhs) const
    {
        return a < rhs.a || (a == rhs.a && b < rhs.b);
    }
}F[N];

LL ask(LL a, int b, LL cp)
{
    LL ret = 1;
    for(int i = 0; i < b; i ++)
    {
        if(ret > (cp + a + i - 1) / (a + i)) return cp + 1;//就是这里犯的2b错误，我当时用的乘法然后判断大小，果断超LL了。。。
        ret *= (a + i); 
    }
    return ret;
}

int solve(LL q)
{
    int ret = 0;
    for(int i = 1; i <= 6; i ++)
    {
        q *= i;
        LL l = 1, r = INF, m;
        while(l <= r)
        {
            m = (l + r) >> 1;
            LL cmp = ask(m, i, q);
            if(cmp == q)
            {
                if(i >= m) break;
                F[ret].a = m + i - 1;
                F[ret ++].b = i;
                if(m - 1 != i)
                {
                    F[ret].a = m + i - 1;
                    F[ret ++].b = m - 1;
                }
                break;
            }
            else if(cmp > q) r = m - 1;
            else l = m + 1;
        }
    }
    return ret;
}

int main()
{
    //freopen("in.txt", "r", stdin);
    //freopen("out.txt", "w", stdout);
    int t, cnt; LL m;
    c[1][0] = 1;c[1][1] = 1;
    for(int i = 2; i < N; i ++)
    {
        c[i][0] = c[i][i] = 1;
        for(int j = 1; j < i; j ++)
        {
            c[i][j] = c[i-1][j] + c[i-1][j-1];
        }
    }
    scanf("%d", &t);
    while(t --)
    {
        cin >> m;
        cnt = solve(m);
        for(int i = 7; i < N; i ++)
        {
            for(int j = 7; j <= i / 2; j ++)
            {
                if(c[i][j] > m) break;
                if(c[i][j] == m)
                {
                    F[cnt].a = i;
                    F[cnt ++].b = j;
                    if(i - j != j)
                    {
                        F[cnt].a = i;
                        F[cnt ++].b = i - j;
                    }
                }
            }
        }
        sort(F, F + cnt);
        printf("%d\n", cnt);
        for(int i = 0; i < cnt; i ++)
        {
            if(i)putchar(' ');
            F[i].pt();
        }puts("");
    }
}