西天取经之杨辉三角形

思路源自yxc

经历了99八十一难终于花了一晚上理解了杨辉三角形
时间复杂度O(16logn) y总牛逼
思路是枚举每一斜列二分查找n
左边界是2k 因为 C(2k,k)是每一斜列的第一个
有边界是max(n, l)与l取最大值只是特判
因为C(n, k)至少比 n 大有可能会取到 n**
            1  ---> C(0, 0)
          1 
        1   2  ---> C(2, 1)
      1   3                             ---> C(2n, n)
    1   4   6  ---> C(4, 2)
  1   5   10
1   6   15  20 ---> C(6, 3)
 

#include<iostream>
using namespace std;

int n;

long long C(int a, int b){
    long long res = 1;
    for(int i = a, j = 1; j <= b; j++, i--){
        res = res * i / j;
        if(res > n) return res;
    }
    return res;
}

bool check(int k){
    int l = k * 2, r = max(l, n);
    while(l < r){
        int mid = l + r >> 1;
        if(C(mid, k) >= n) r = mid;
        else l = mid + 1;
    }
    if(C(l, k) != n) return false;
    cout << (long long)l * (l + 1) / 2 + k + 1; 
    //C(l, k)这一个是满足条件的 前面有l列在加上k + 1就是本身位置
    return true;
}

int main(){
    cin >> n;
    for(int i = 16; ;i--){
        if(check(i)) break;
    }
    return 0;
}