组合数单调性+二分。这道题犯了一个很二逼的错误,调了很久很久才能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("");
}
}