UVa12333 Revenge of Fibonacci

高精度 trie

暴力预处理出前100000个fibonacci数，将每个数的前40位数字串插入到trie中，记录每个结点最早可以由哪个数字串到达。

然后依次回答询问即可。

存fibonacci数的数组当然要滚动起来。

时限是10秒。本机试着卡了常数后跑了18秒（这个高精度写法好像本来就很慢），干脆交一波，UVA上7s AC

　　　　　　　　↑for循环i+=4 和直接一个register跑出来一样快，果然WC2017那道题只能在特定机子上做吗233

　　　　　　　　↑本机CPU是Intel I5-4590

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#include<cstdlib>
#include<time.h>
using namespace std;
const int mxn=200010;
struct Num{
    int x[24501];
    int len;
    friend Num operator + (Num a,Num b){
        a.len=max(a.len,b.len);
//        printf("len:%d
",a.len);
        int ed=a.len/4*4;
        for(int i=1;i<=ed;i+=4){
            a.x[i]+=b.x[i];
            a.x[i+1]+=b.x[i+1];
            a.x[i+2]+=b.x[i+2];
            a.x[i+3]+=b.x[i+3];
        }
        for(int i=ed+1;i<=a.len;i++)
            a.x[i]+=b.x[i];

//        for(register int i=1;i<=a.len;i++) a.x[i]+=b.x[i];
        
        for(int i=1;i<=a.len;i++){
            if(a.x[i]>=10){
                a.x[i+1]+=a.x[i]/10;
                a.x[i]%=10;
            }
        }
        if(a.x[a.len+1]) ++a.len;
        return a;
    }
}a[2];
struct Trie{
    int t[4850007][10];
    int end[4850007];
    int cnt,rt;
    void init(){
        memset(end,0x3f,sizeof end);
        cnt=rt=1;
    }
    void insert(int id,Num a){
        int ed=a.len;
        int st=max(1,a.len-40+1);//取高位40位 
        int now=rt;
        for(int i=ed;i>=st;i--){
            if(!t[now][a.x[i]])t[now][a.x[i]]=++cnt;
            now=t[now][a.x[i]];
            end[now]=min(end[now],id);
        }
    }
    int query(char *s){
        int now=rt,len=strlen(s),tmp;
        for(int i=0;i<len;i++){
            tmp=s[i]-'0';
            if(!t[now][tmp])return -1;
            now=t[now][tmp];
            if(end[now]>100000)return -1;
        }
        return end[now];
    }
}tr;
void init(){
    a[0].len=1;a[0].x[1]=1;
    a[1].len=1;a[1].x[1]=1;
    tr.insert(0,a[0]);
    tr.insert(1,a[1]);
    for(int i=2;i<100000;i++){
//        printf("init:%d 
",i);
        int tmp=i&1;
        a[tmp]=a[tmp]+a[tmp^1];
//        for(int j=1;j<=a[tmp].len;j++)printf("%d ",a[tmp].x[j]);
//        printf("%d cnt:%d
",i,tr.cnt);
        tr.insert(i,a[tmp]);
    }
    return;
}
char s[mxn];
int n;
int main(){
    int i,j;
//    clock_t a=clock();
//    printf("%d
",a);
    tr.init();
    init();
//    clock_t b=clock();
//    printf("%d
",b);
//    printf("fin: %d 
",b-a);
    scanf("%d",&n);
    for(i=1;i<=n;i++){
        scanf("%s",s);
        printf("Case #%d: %d
",i,tr.query(s));
    }
    return 0;
}