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

版权声明：本文博客原创文章。博客，未经同意，不得转载。