CCF NOI1097 数列

问题链接：CCF NOI1097 数列。

---

时间限制: 1000 ms  空间限制: 262144 KB

题目描述 

  璋璋的数学又不及格了。他的数学老师很生气，罚他放学留下来做计算题。守老师给了璋璋一个递推数列：f[1]=1,f[2]=1,f[i]=3*f[i-1]+2*f[i-2]，让他求f[n]。
  本来这道题目是不难的，可是守老师为了惩罚璋璋，给了他一个非常大的n(n<=500)，这下子可难倒了璋璋。璋璋只有向身为计算机高手的你求助了，你能帮助他完成这个问题吗？
  最后请输出f[n]。

输入

  第1行为n。

输出

  仅一行，f[n]。

样例输入

4
样例输出

17

数据范围限制

  n<=500

提示

 

---

问题分析

  这是计算菲波拉契数列第n项的问题，是一个大数计算问题。

  大数计算一种是构建一个大数类进行计算，另外一种是直接计算。

  人们使用阿拉伯数字，据说是印度人发明的。需要注意的一点是，阿拉伯数字的高位在左边，阅读上是从左到右，而计算上人们则是从低位算起。

  大数可以放在数组中，为了计算上的方便，应该把低位放在下标小的地方，高位放在下标大的地方。

  计算时下标从小到大，输出时要倒过来。

程序说明

  程序中采用直接模拟的方法。

要点详解

• 使用宏定义可以使得代码可阅读性增强。
• 加法计算需要考虑进位，实际上每一位的加法是三个数相加。
  

---

参考链接：（略）。

---

100分通过的C语言程序：

#include <stdio.h>#include <string.h>#define N 500#define BASE 10char f1[N], f2[N], temp[N];int len;int main(void){    int n, sum, c, i, j;    scanf("%d", &n);    if(n <= 2) {        printf("1\n");        return 0;    }    memset(f1, 0, sizeof(f1));    memset(f2, 0, sizeof(f2));    memset(temp, 0, sizeof(temp));    len = 1;    f1[0] = 1;    f2[0] = 1;    for(i=3; i<=n; i++) {        c = 0;        for(j=0; j<len; j++) {            sum = c + 2 * f1[j] + 3 * f2[j];            temp[j] = sum % BASE;            c = sum / BASE;        }        while(c > 0) {            temp[len++] = c % BASE;            c /= BASE;        }        memcpy(f1, f2, len);        memcpy(f2, temp, len);    }    for(i=len-1; i>=0; i--)        printf("%d", f2[i]);    printf("\n");    return 0;}