[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里面一堆奇怪语言也是感人>_<

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cstring>
 4 #define ll long long
 5 #define d double
 6 using namespace std;
 7 const int maxn=1002333;
 8 const int modd=1000000000;
 9 ll f[233];
10 int dl[maxn],l,r;
11 ll a[32333],c[32333],b[32333];
12 int i,j,k,n,m,lena,lenc;
13
14 inline d xl(int k,int j){//k<j<i-2
15     return ( f[j]*(j+1)-f[k]*(k+1) )/ (d)(f[j]-f[k]);
16 }
17 int ra,fh;char rx;
18 inline int read(){
19     rx=getchar(),ra=0,fh=1;
20     while((rx<‘0‘||rx>‘9‘)&&rx!=‘-‘)rx=getchar();
21     if(rx==‘-‘)fh=-1,rx=getchar();
22     while(rx>=‘0‘&&rx<=‘9‘)ra*=10,ra+=rx-48,rx=getchar();return ra*fh;
23 }
24 inline void multoa(){
25     register int i,j,k;
26     memset(b,0,(lena*2+1)<<3);
27     for(i=1;i<=lena;i++)for(j=1,k=i;j<=lena;j++,k++){
28         b[k]+=a[i]*a[j];
29         if(b[k]>=modd)b[k+1]+=b[k]/modd,b[k]%=modd;
30     }
31     for(k=lena<<1;!b[k];k--);
32     memcpy(a,b,(k+1)<<3);lena=k;
33 }
34 inline void multoc(){
35     register int i,j,k;
36     memset(b,0,(lena+lenc+1)<<3);
37     for(i=1;i<=lena;i++)for(j=1,k=i;j<=lenc;j++,k++){
38         b[k]+=a[i]*c[j];
39         if(b[k]>=modd)b[k+1]+=b[k]/modd,b[k]%=modd;
40     }
41     for(k=lena+lenc;!b[k];k--);
42     memcpy(c,b,(k+1)<<3);lenc=k;
43 }
44 inline void poi(int b){//计算4^b
45     c[lenc=1]=1;a[lena=1]=4;
46     while(b){
47         if(b&1)multoc();
48         b>>=1;
49         if(b)multoa();
50     }
51 }
52
53 int main(){
54     n=read();
55     f[1]=1,f[2]=2,f[3]=3;
56     l=1,r=0;
57     for(i=4;i<=25;i++){
58         while(l<r&&xl(dl[r-1],dl[r])>=xl(dl[r],i-3))r--;
59         dl[++r]=i-3;
60         while(l<r&&i>xl(dl[l],dl[l+1]))l++;
61         f[i]=f[dl[l]]*(i-dl[l]-1);
62         if(f[i-1]+1>f[i])f[i]=f[i-1]+1;
63 //      printf(" i:%d   %lld\n",i,f[i]);
64     }
65     if(n<=25)printf("%lld\n",f[n]);else{
66         int x=n%5+15;
67         poi((n-x)/5);
68         a[lena=1]=f[x];
69         multoc();
70         for(printf("%lld",c[lenc]),i=lenc-1;i>0;i--){
71             for(j=modd/10;j;j/=10)if(c[i]<j)putchar(‘0‘);
72             printf("%lld",c[i]);
73         }puts("");//printf("  %d\n",lenc);
74     }
75     return 0;
76 }