hdu 1316

Description

我们定义斐波那契数列如下： 

f1=1 

f2=2 

   f(n)=f(n-1)+f(n-2)(n>=3) 

     

现在，给定两个数a和b，计算有多少个斐波那契数列中的数在a和b之间（包含边界）。 

Input

输入包含多组测试数据，每组测试数据都是两个非负整数a和b，当a和b都等于0时，程序结束。0<=a<=b<=10^100,注意a和b都不会出现前导0的情况（即首位为0）。

Output

输出在a和b之间（满足a<=f(i)<=b)数的数量。

Sample Input

   10 100 

   1234567890 9876543210 

   0 0  

Sample Output

   5 

    4 

单纯的高精度加法，先打表，再找到第一个大于或等于下界的斐波拉契数，再找到最后一个小于等于上界的斐波拉契数，相减再加一。

代码

#include <iostream>
#include <string.h>
using namespace std;
int a[1100][110],aa[1100][110];
int head[1100];
void dabiao()
{
  memset(a,0,sizeof(a));
  a[1][105]=1;
  a[2][105]=2;
  int c;
  for(int i=3;i<=1050;i++)
  {
   for(int j=105;j>=1;j--)
   {
    a[i][j]=a[i-1][j]+a[i-2][j];
   }
   for(int j=105;j>=1;j--)
   {
    c=a[i][j]/10;
   if(c!=0)
   {
    a[i][j]=a[i][j]%10;
    a[i][j-1]=a[i][j-1]+c;
   }
   }
  }
  for(int i=1;i<=1050;i++)
  {
   for(int j=0;j<=105;j++)
   {
    if(a[i][j]!=0)
    {
     head[i]=j;
     break;
    }
   }
  }
}
int compare1(char *string,int len,int p[],int n)
{
  if(n>len)  return 1;
  if(n<len)  return 0;
  for(int i=0;i<=n-1;i++)
  {
   if(string[i]-'0'<p[i]) return 1;
   if(string[i]-'0'>p[i]) return 0;
  }
  return 1;
}
int compare2(char *string,int len,int p[],int n)
{
  if(n<len)  return 1;
  if(n>len)  return 0;
  for(int i=0;i<=n-1;i++)
  {
   if(string[i]-'0'>p[i]) return 1;
   if(string[i]-'0'<p[i]) return 0;
  }
  return 1;
}
int main()
{
    char str1[110],str2[110];
    int len1,len2,res1,res2;
    dabiao();
    for(int i=1;i<=1050;i++)
    {
     for(int j=head[i];j<=105;j++)
     {
      aa[i][j-head[i]]=a[i][j];
     }
    }
    while(cin>>str1>>str2)
    {
     len1=strlen(str1);
     len2=strlen(str2);
     if(len1==1&&len2==1&&str1[0]=='0'&&str2[0]=='0')
      break;
     for(int i=1;i<=1050;i++)
     {
      if(compare1(str1,len1,aa[i],105-head[i]+1))
      {
        res1=i;
        break;
        }
     }
     for(int i=1050;i>=1;i--)
     {
      if(compare2(str2,len2,aa[i],105-head[i]+1))
      {
       res2=i;
       break;
      }
     }
    cout <<res2-res1+1<< endl;
    }
    return 0;
}