[bzoj3287] Mato的刷屏计划

　　第一眼以为是傻逼斜率优化>_<

　　f[i]表示按i次最多可输出字符数。。f[i]=max{ f[i-1]+1,(i-j-1)*f[j] }，j<i-2

　　结果n在100+的时候就喜闻乐见地爆了longlong

　　根据网上题解可得（T_T）。。这题大概是要FFT优化？（跟着ccz大爷刷题果然高风险TAT

　　然而题解表示，这题求出f的前几项后，令x=n%5+15,最终答案f[n]=f[x]*4^((n-x)/5)..........

　　我死活没看出来这是为啥。。。

　　我显然是连FFT都不会的傻逼。。只好写了个压位高精度快速幂。。结果跑了#3。。233

　　突然觉得正解显然不是FFT = =。。另外status里面一堆奇怪语言也是感人>_<

 #include<cstdio>
 #include<iostream>
 #include<cstring>
 #define ll long long
 #define d double
 using namespace std;
 const int maxn=;
 const int modd=;
 ll f[];
 int dl[maxn],l,r;
 ll a[],c[],b[];
 int i,j,k,n,m,lena,lenc;

 inline d xl(int k,int j){//k<j<i-2
     return ( f[j]*(j+)-f[k]*(k+) )/ (d)(f[j]-f[k]);
 }
 int ra,fh;char rx;
 inline int read(){
     rx=getchar(),ra=,fh=;
     while((rx<''||rx>'')&&rx!='-')rx=getchar();
     if(rx=='-')fh=-,rx=getchar();
     while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra*fh;
 }
 inline void multoa(){
     register int i,j,k;
     memset(b,,(lena*+)<<);
     for(i=;i<=lena;i++)for(j=,k=i;j<=lena;j++,k++){
         b[k]+=a[i]*a[j];
         if(b[k]>=modd)b[k+]+=b[k]/modd,b[k]%=modd;
     }
     for(k=lena<<;!b[k];k--);
     memcpy(a,b,(k+)<<);lena=k;
 }
 inline void multoc(){
     register int i,j,k;
     memset(b,,(lena+lenc+)<<);
     for(i=;i<=lena;i++)for(j=,k=i;j<=lenc;j++,k++){
         b[k]+=a[i]*c[j];
         if(b[k]>=modd)b[k+]+=b[k]/modd,b[k]%=modd;
     }
     for(k=lena+lenc;!b[k];k--);
     memcpy(c,b,(k+)<<);lenc=k;
 }
 inline void poi(int b){//计算4^b
     c[lenc=]=;a[lena=]=;
     while(b){
         if(b&)multoc();
         b>>=;
         if(b)multoa();
     }
 }

 int main(){
     n=read();
     f[]=,f[]=,f[]=;
     l=,r=;
     for(i=;i<=;i++){
         while(l<r&&xl(dl[r-],dl[r])>=xl(dl[r],i-))r--;
         dl[++r]=i-;
         while(l<r&&i>xl(dl[l],dl[l+]))l++;
         f[i]=f[dl[l]]*(i-dl[l]-);
         if(f[i-]+>f[i])f[i]=f[i-]+;
 //      printf(" i:%d   %lld\n",i,f[i]);
     }
     if(n<=)printf("%lld\n",f[n]);else{
         int x=n%+;
         poi((n-x)/);
         a[lena=]=f[x];
         multoc();
         for(printf("%lld",c[lenc]),i=lenc-;i>;i--){
             for(j=modd/;j;j/=)if(c[i]<j)putchar('');
             printf("%lld",c[i]);
         }puts("");//printf("  %d\n",lenc);
     }
     return ;
 }