vijos P1375 大整数（高精不熟的一定要做！）

/*
我尼玛这题不想说啥了
亏了高精写的熟.....
加减乘除max都写了
高精二分
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#define maxn 1010
#define memcpy(a,b); for(int i=0;i<=1000;i++)a[i]=b[i];
using namespace std;
int n[maxn],len,l[maxn],r[maxn],mid[maxn],ans[maxn];
char s[maxn];
void Plus(int a[maxn],int b[maxn]){
    int L=max(a[],b[]);
    int c[maxn];memset(c,,sizeof(c));
    for(int i=;i<=L;i++)
        c[i]=a[i]+b[i];
    for(int i=;i<=L;i++)
        if(c[i]>){c[i+]++;c[i]%=;}
    if(c[L+])L++;c[]=L;
    memcpy(a,c);
}
void Sub(int a[maxn]){
    a[]--;
    int p=;
    while(a[p]<){
        a[p]=;p++;a[p]--;
    }
    if(a[a[]]==)a[]--;
}
void Mul(int a[maxn],int b[maxn]){
    int c[maxn];memset(c,,sizeof(c));
    int l1=a[],l2=b[],l3=l1+l2;
    for(int i=;i<=l1;i++){
          int x=;
          for(int j=;j<=l2;j++){
              c[j+i-]+=a[i]*b[j]+x;
            x=c[i+j-]/;
            c[i+j-]=c[i+j-]%;
        }
        c[i+l2]=x;
    }
    int k=;
    for(int i=l3;i>=;i--)
        if(c[i]){k=i;break;}
    c[]=k;memcpy(a,c);
}
void Div(int a[maxn],int x){
    int b[maxn],c=;memset(b,,sizeof(b));
    for(int i=a[];i>=;i--){
        b[i]=(c*+a[i])/x;c=a[i]%x;
    }
    int k=;
    for(int i=a[];i>=;i--)
        if(b[i]){k=i;break;}
    b[]=k;memcpy(a,b);
}
bool Cmp(int a[maxn],int b[maxn]){
    if(a[]<b[])return ;
    else if(a[]>b[])return ;
    for(int i=a[];i>=;i--){
        if(a[i]<b[i])return ;
        if(a[i]>b[i])return ;
    }
    return ;
}
void Cal(int x[maxn]){
    int a[maxn],b[maxn],c[maxn],d[maxn],e[maxn];
    e[]=;e[]=;memset(a,,sizeof(a));
    memcpy(b,x);memcpy(c,x);memcpy(d,x);
    Mul(b,x);Mul(b,x);Mul(c,x);Mul(d,e);
    Plus(a,b);Plus(a,c);Plus(a,d);
    memcpy(x,a);
}
bool Judge(int a[maxn]){
    int b[maxn];memcpy(b,a);
    Cal(b);
    return Cmp(b,n);
}
void cal(int a[maxn]){
    int b[maxn],c[maxn];
    memcpy(b,l);memcpy(c,r);
    Plus(b,c);Div(b,);
    memcpy(a,b);
}
int main()
{
    scanf("%s",s);
    len=strlen(s);
    for(int i=;i<=len;i++)
        n[i]=s[len-i]-'';
    n[]=len;
    l[]=;r[]=;
    for(int i=;i<=;i++)r[i]=;
    while(Cmp(l,r)){
        cal(mid);int a[maxn];
        memset(a,,sizeof(a));
        a[]=;a[]=;
        if(Judge(mid)){
            memcpy(l,mid);
            Plus(l,a);
            memcpy(ans,mid);
        }
        else{
            memcpy(r,mid);
            Sub(r);
        }
    }
    for(int i=max(ans[],);i>=;i--)
        printf("%d",ans[i]);
    return ;
}