POJ 2413 How many Fibs?

大数区间Fibonacci数的个数


const int maxn = 9999 ; // 每位储存4位
const int dig = 4 ;

class BigNum{
 public  :
       int  a[60] ;
       int  len ;
       BigNum(){len = 1 ; memset(a , 0 , sizeof(a)) ;}
       BigNum(const int) ;
       BigNum(const char*) ;
       BigNum operator +(const BigNum &) const ;
       bool operator >= (const BigNum &) const ;
       friend ostream& operator <<(ostream &out , const BigNum&T) ;
};

BigNum::BigNum(const int T){
      int t = T ;
      len = 0 ;  memset(a , 0 , sizeof(a)) ;
      while(t > maxn){
            a[len++] = t - t/(maxn+1)*(maxn+1) ;
            t = t/(maxn+1) ;
      }
      a[len++] = t ;
}

BigNum::BigNum(const char* s){
      memset(a , 0 , sizeof(a)) ;
      int t , i , j , k , id = 0 , l = strlen(s) ;
      len = l/dig ;
      if(l%dig) len++ ;
      for(i = l-1 ; i >= 0 ; i -= dig){
          t = 0 ;
          k = max(0  , i - dig + 1);
          for(j = k ; j <= i ; j++)  t = t*10 + s[j] - '0' ;
          a[id++] = t ;
      }
}

ostream& operator <<(ostream &out , const BigNum&T){
      printf("%d" , T.a[T.len-1]) ;
      for(int i = T.len-2 ; i >= 0 ; i--)  printf("%04d" ,T.a[i]) ;
      return out ;
}

BigNum BigNum::operator + (const BigNum &T) const{
       BigNum s ;
       s.len = max(len , T.len) ;
       for(int i = 0 ; i < s.len ; i++){
           s.a[i] += a[i] + T.a[i] ;
           if(s.a[i] > maxn){
                s.a[i+1]++ ;
                s.a[i] %= (maxn+1) ;
           }
       }
       if(s.a[s.len]!=0) s.len++ ;
       return s ;
}

bool BigNum::operator>= (const BigNum &T) const{
     if(len > T.len)  return 1 ;
     else if(len == T.len){
           int i = len-1 ;
           while(i >= 0 && a[i] == T.a[i]) i-- ;
           if(i >= 0 && a[i] < T.a[i]) return 0 ;
           else return 1 ;
     }
     else  return 0 ;
}

vector <BigNum> lis ;

int  main(){
     lis.clear()  ;
     lis.push_back(BigNum("1")) ;
     lis.push_back(BigNum("2")) ;
     int i  , s ;
     for(i = 2 ;  ; i++){
          BigNum b = lis[i-1] + lis[i-2] ;
          if(b.len > 100)  break  ;
          lis.push_back(b) ;
     }
     BigNum l , r ;
     char a[108] , b[108] ;
     while(scanf("%s%s" ,a , b)){
          if(strcmp(a , "0")== 0 &&  strcmp(b , "0") == 0)  break ;
          l = BigNum(a) , r = BigNum(b) ;
          for(i = s = 0  ; i < lis.size() ; i++){
              if(lis[i] >= l && r >= lis[i]) s++ ;
          }
          printf("%d\n" , s) ;
     }
     return 0 ;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值