hdoj 1316 How Many Fibs?(大数)

题目大意：给出a ，b(a <= b <= 10^100)，求区间[a,b]内含fib数的个数。

思路：初始化fib数列至100位（第480个）。从头开始比较 小于a的的个数，从尾开始比较小于等于b的个数，相减即为所求。

 

需先把a，b化为与数列格式相同的整型。比较的时候先比较长度，长度相同，从最高为开始比较。

 

#include <iostream>using namespace std;char a[101],b[101];int A[15],B[15];int f[481][15]={0};int len[481];void init(){ int i,j,l; f[0][0]=1; f[1][0]=2; len[0]=len[1]=1; for(i=2;i<=480;i++) { l=len[i-1]>len[i-2]? len[i-1]:len[i-2]; for(j=0;j<l;j++) { f[i][j]+=f[i-1][j]+f[i-2][j]; if(f[i][j]>=100000000) { f[i][j+1]++; f[i][j]%=100000000; } } if(f[i][j]) len[i]=l+1; else len[i]=l; } return;}int bigger(int l){ int i,j; for(i=0;i<=480;i++) { if(len[i]>l) return i-1; else if(len[i]==l) { for(j=l-1;j>=0;j--) if(f[i][j]>A[j]) return i-1; else if(f[i][j]<A[j]) break; if(j==-1 && f[i][0]==A[0]) return i-1; } }}int smaller(int l){ int i,j; for(i=480;i>=0;i--) { if(len[i]<l) return i; else if(len[i]==l) { for(j=l-1;j>=0;j--) if(f[i][j]<B[j]) return i; else if(f[i][j]>B[j]) break; if(j==-1 && f[i][0]==B[0]) return i; } }} int main(){ int i,j,t,ans; int la,lb,lA, lB; init(); while(scanf("%s%s",&a,&b)!=EOF) { la=strlen(a); lb=strlen(b); if(la==1 && a[0]=='0' && lb==1 && b[0]=='0') break; memset(A,0,sizeof(A)); memset(B,0,sizeof(B)); for(i=la-1,lA=0;i>=0;lA++,i-=8) for(j=0,t=1;j<8 && i-j>=0;j++,t*=10) A[lA]+=(a[i-j]-'0')*t; for(i=lb-1,lB=0;i>=0;lB++,i-=8) for(j=0,t=1;j<8 && i-j>=0;j++,t*=10) B[lB]+=(b[i-j]-'0')*t; ans=smaller(lB)-bigger(lA); printf("%d/n",ans); } return 0;}