VIjos 大整数(超级恶心的高精)

背景

很久很久以前,有个整数很大很囧

描述

一个k(1<=k<=80)位的十进制正整数N,就是所谓的大整数.
请你设计程序,对于给出的某一个大整数N,找到满足p^3+p^2+3p<=n的p的最大值.

格式

输入格式

输入数据只有一行，是一个K位的大整数N，行首行未无多余空格

输出格式

输出第一行为你所找到的P最大值，行首行末别乱加东西

样例1

样例输入1

1000000000000001000000000000003000000000000001

Copy

样例输出1

1000000000000000

  1 /*
  2     一道不错的高精练手题
  3     先确定p的位数
  4     然后从大到小找合适的p
  5     k最多80位
  6     p那最多 28位 
  7     假设每一位都用0-9试一遍
  8     最多试不超过300遍 
  9     即使加上高精复杂度也不会太高
 10 */
 11 #include<cstdio>
 12 #include<cstring> 
 13 #include<iostream>
 14 #define MAXN 1010
 15 
 16 using namespace std;
 17 
 18 char s[MAXN];
 19 
 20 int a[MAXN],b[MAXN],c[MAXN],e[MAXN],n[MAXN],_w;
 21 
 22 inline bool judge() {
 23     if(e[0]<n[0]) return true;
 24     if(e[0]>n[0]) return false;
 25     for(int i=n[0];i>=1;i--) {
 26         if(e[i]>n[i]) return false;
 27         if(e[i]<n[i]) return true;
 28     }
 29     return true;
 30 }
 31 
 32 inline void check() {
 33     for(int i=1;i<=a[0];i++) {
 34         if(a[i]<0) {
 35             a[i+1]--;
 36             a[i]+=10;
 37         }
 38     }
 39     while(a[0]>1&&!a[a[0]]) a[0]--;
 40     return;
 41 }
 42 
 43 inline void _pluss() {
 44     a[_w]++;
 45     a[0]++;
 46     for(int i=1;i<=a[0];i++) {
 47         if(a[i]>=10) {
 48             a[i+1]++;
 49             a[i]-=10;
 50         } 
 51     }
 52     while(a[0]>0&&!a[a[0]]) a[0]--;
 53     return;
 54 }
 55 
 56 inline void __mul() {
 57     for(int i=0;i<=a[0];i++) b[i]=a[i];
 58     c[0]=a[0]+b[0]+1;
 59     for(int i=1;i<=a[0];i++)
 60       for(int j=1;j<=b[0];j++)
 61           c[i+j-1]+=a[i]*b[j];
 62     for(int i=1;i<=c[0];i++) 
 63       if(c[i]>=10) {
 64           c[i+1]+=c[i]/10;
 65           c[i]%=10;
 66       }
 67     while(c[0]&&!c[c[0]]) c[0]--;
 68     for(int i=c[0];i>=0;i--) b[i]=c[i],c[i]=0;
 69     c[0]=a[0]+b[0]+1;
 70     for(int i=1;i<=a[0];i++)
 71       for(int j=1;j<=b[0];j++)
 72         c[i+j-1]+=a[i]*b[j];
 73     for(int i=1;i<=c[0];i++) 
 74       if(c[i]>=10) {
 75         c[i+1]+=c[i]/10;
 76         c[i]%=10;
 77     }
 78     while(c[0]&&!c[c[0]]) c[0]--;
 79     e[0]=max(e[0],c[0])+1;
 80     for(int i=1;i<=e[0];i++) {
 81         e[i]+=c[i];
 82         if(e[i]>=10) {
 83             e[i+1]+=e[i]/10;
 84             e[i]%=10;
 85         }
 86     }
 87     while(e[0]&&!e[e[0]]) e[0]--;
 88     memset(c,0,sizeof c);
 89     memset(b,0,sizeof b);
 90 }
 91 
 92 inline void _mul() {
 93     for(int i=0;i<=a[0];i++) b[i]=a[i];
 94     c[0]=a[0]*2+1;
 95     for(int i=1;i<=a[0];i++)
 96       for(int j=1;j<=b[0];j++)
 97           c[i+j-1]+=a[i]*b[j];
 98     for(int i=1;i<=c[0];i++) 
 99       if(c[i]>=10) {
100           c[i+1]+=c[i]/10;
101           c[i]%=10;
102       }
103     while(c[0]&&!c[c[0]]) c[0]--;
104     e[0]=max(e[0],c[0])+1;
105     for(int i=1;i<=e[0];i++) {
106         e[i]+=c[i];
107         if(e[i]>=10) {
108             e[i+1]+=e[i]/10;
109             e[i]%=10;
110         }
111     }
112     while(e[0]&&!e[e[0]]) e[0]--;
113     memset(c,0,sizeof c);
114     memset(b,0,sizeof b);
115 }
116 
117 inline void pluss() {
118     for(int i=1;i<=a[0];i++) b[i]=0;
119     b[0]=a[0]+1;
120     for(int i=1;i<=b[0];i++) {
121         b[i]+=3*a[i];
122         if(b[i]>=10) {
123             b[i+1]+=b[i]/10;
124             b[i]%=10;
125         }
126     }
127     while(b[0]&&!b[b[0]]) b[0]--;
128     e[0]++;
129     for(int i=1;i<=b[0];i++) {
130         e[i]+=b[i];
131         if(e[i]>=10) {
132             e[i+1]+=e[i]/10;
133             e[i]%=10;
134         }
135     }
136     while(e[0]&&!e[e[0]]) e[0]--;
137 }
138 
139 int main() {
140     scanf("%s",s);
141     int len=strlen(s);
142     n[0]=len;
143     for(int i=1;i<=len;i++) n[i]=s[len-i]-48;
144     int _len=len/3+1;//因为有p^3 
145     a[0]=_len;
146     a[_len]=0;
147     _w=_len;
148     while(true) {
149         if(!judge()) {
150             a[_w]--;
151             check();
152             _w--;
153             if(_w==0) break;
154             a[_w]=1;
155         }
156         else _pluss();
157         memset(e,0,sizeof e);
158         __mul();
159         _mul();
160         pluss();
161     }
162     for(int i=a[0];i>=1;i--) printf("%d",a[i]);
163     return 0;
164 } 

代码