【二分法】 HDU 2446 Shell Pyramid-优快云博客

意甲冠军：非常多，形成一个金字塔球

文章x层 x*（x+1）/ 2 球

给你个S 金字塔的一层代表第一数字向下S球

它是其中这层中的第几行第几列

公式 1 ： x*（x+1）*（x+2）/ 6

公式 2 ：x*（x+1）/ 2

公式1为公式2 的前n项和

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <string>
#include <iostream>
#include <algorithm>
#include <sstream>
#include <cmath>
using namespace std;
#include <queue>
#include <stack>
#include <vector>
#include <deque>
#include <set>
#include <map>
#include <time.h>;
#define cler(arr, val)    memset(arr, val, sizeof(arr))
#define FOR(i,a,b)  for(int i=a;i<=b;i++)
#define IN   freopen ("in.txt" , "r" , stdin);
#define OUT  freopen ("out.txt" , "w" , stdout);
typedef long long  LL;
const int MAXN = 10323;
const int MAXM = 201;
const int INF = 0x3f3f3f3f;
const int mod = 1000000007;
LL getnum1(LL num)
{
    LL x=num;
    LL y=num+1;
    LL z=num+2;
    if(x%2==0) x/=2;
    else if(y%2==0) y/=2;
    else if(z%2==0) z/=2;

    if(x%3==0) x/=3;
    else if(y%3==0) y/=3;
    else if(z%3==0) z/=3;
    return x*y*z;
}
LL getnum2(LL num)
{
    return num*(num+1)/2;
}
int main()
{
    #ifndef ONLINE_JUDGE
        freopen("in.txt", "r", stdin);
    #endif
    int t;
    LL n;
    scanf("%d",&t);
    while(t--)
    {
        cin>>n;
        LL l=0,r=3810776;//r表示极限了。

。 while(l<=r)//求第几层 { LL mid=(l+r)>>1;//getnum1(3810776) if(getnum1(mid)>=n) r=mid-1; else l=mid+1; } LL ans1=l; n-=getnum1(l-1); l=0,r=3810776; while(l<=r)//求第几行 { LL mid=(l+r)>>1; if(getnum2(mid)>=n) r=mid-1; else l=mid+1; } n-=getnum2(l-1); LL ans2=l; printf("%I64d %I64d %I64d\n",ans1,ans2,n); } }