【二分法】 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);
}
}