hdu 4099 Revenge of Fibonacci（字典树+大数加法）

最新推荐文章于 2020-08-26 17:18:11 发布

風吹过

最新推荐文章于 2020-08-26 17:18:11 发布

阅读量580

点赞数

CC 4.0 BY-SA版权

分类专栏：数论字符串

本文链接：https://blog.youkuaiyun.com/hlmfjkqaz/article/details/13631925

数论同时被 2 个专栏收录

42 篇文章

订阅专栏

字符串

5 篇文章

订阅专栏

本文介绍了一种利用字典树与大数优化的方法来解决特定问题：即找到一个整数首次出现在斐波那契数列中的位置。通过预处理斐波那契数列并构建字典树，实现了高效的查询机制。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意：给你一个数字字符串，表示一个整数的前几位（<40），求这个整数出现斐波那契数列中的最早的位置n

解析：赤裸裸的字典树，不过要用大数优化，坑爹的是插入时要严格控制小于100000，不能等于；就因为这个wa了好久

#include<iostream>
#include<cstdio>
#include<string.h>
#include<string>
using namespace std;
#define N 100001
int maxx=45;
struct node
{
    int minn;
    node  *next[11];
    node()
    {
        minn=N;
        memset(next,0,sizeof(next));
    }
}*root;

int search(node *root,char *s)
{
    node *p=root;
    for(int i=0;s[i];i++)
    {
        int x=s[i]-'0';
        if(p->next[x]==NULL) return -1;

         p=p->next[x];
         //cout<<x<<"   "<<p->minn;
    }
    //cout<<endl;
    return p->minn;
}
void insert(node *root,int *s,int m,int len)
{
    node *p=root;
    int i;
    for(i=len;i>=0&&i>=len-40;i--)
    {
        int x=s[i];
        //cout<<x;
        if(p->next[x]==NULL)
          p->next[x]=new node;
        p=p->next[x];
        p->minn=min(p->minn,m);
       // if(x<30)
      //  cout<<p->minn<<"********"<<endl;
    }
   // cout<<"   "<<p->minn<<endl;
}

void inint()
{
   int a[70],b[70],c[70];
   memset(a,0,sizeof(a));
   memset(b,0,sizeof(b));
   memset(c,0,sizeof(c));
   int i,j,k,ma,mb,mc,w;
   a[0]=b[0]=c[0]=1;
   ma=mb=mc=0;
   insert(root,c,0,0);
   for(i=2;i<N-1;i++)
   {
       memcpy(a,b,sizeof(a));
       memcpy(b,c,sizeof(b));
       memset(c,0,sizeof(c));
       ma=mb;
       mb=mc;
       if(ma>=60)
       {
          for(k=0;k<70;k++)
          {
              if(k<=ma)
              a[k]=a[k+2];
              else a[k]=0;
          }
           for(k=0;k<70;k++)
          {
              if(k<=mb)
              b[k]=b[k+2];
              else b[k]=0;
          }
          ma-=2;
          mb-=2;
       }
       mc=-1;
       for(j=0;j<=ma;j++)
       {
          c[j]+=a[j]+b[j];
       }
       for(j;j<=mb;j++)
         c[j]+=b[j];
       for(j=0;j<70;j++)
       {
          w=c[j];
          c[j]=w%10;
          c[j+1]+=w/10;
          if(c[j])
          mc=max(mc,j);
       }
       insert(root,c,i,mc);
       /*
       for(j=mc;j>=0;j--)
        cout<<c[j];
        cout<<endl;
        */
   }
}
int main()
{
    root=new node;
    inint();
    int i,t=1,T,f;
    char s[100];
    scanf("%d",&T);
    while(T--)
    {
       scanf("%s",&s);
       f=search(root,s);
       printf("Case #%d: %d\n",t++,f);
    }
    return 0;
}