大数

Hat's Fibonacci http://acm.hdu.edu.cn/showproblem.php?pid=1250

java

import java.math.BigInteger;
import java.util.Scanner;

public class Main {
    public static void main(String args[]){
        Scanner cin=new Scanner(System.in);
        BigInteger dp[]=new BigInteger[10010];
        dp[1]=dp[2]=dp[3]=dp[4]=BigInteger.ONE;
        for(int i=5;i<=10000;i++){
            dp[i]=dp[i-1];
            for(int j=2;j<=4;j++){
                dp[i]=dp[i].add(dp[i-j]);
            }
        }
        while(cin.hasNext()){
            int n=cin.nextInt();
            System.out.println(dp[n]);
        }
    }
}

 高精度，存不了，只能递推算

#include<cstdio>
#include<cstring>
#include<algorithm>
#define mt(a,b) memset(a,b,sizeof(a))
using namespace std;
const int M=2020;
class Hp { //高精度类
    int len,s[M];
public:
    void init() {
        len=1;
        mt(s,0);
    }
    void init(char ch[]) {
        len=1;
        mt(s,0);
        int i=0;
        while(ch[i]=='0'&&ch[i]!=0) i++;
        if(ch[i]!=0) {
            len=strlen(ch)-i;
            for(i=0; i<len; i++) {
                s[i]=ch[len-i-1]-48;
            }
        }
    }
    void init(int x) {
        len=1;
        mt(s,0);
        while(x) {
            s[len-1]=x%10;
            x/=10;
            if(x) len++;
        }
    }
    void init(const Hp& a) {
        len=a.len;
        for(int i=0; i<M; i++) {
            s[i]=a.s[i];
        }
    }
    void print() { //输出
        int i=len-1;
        while(s[i]==0&&i>0) i--;
        for(; i>=0; i--) {
            printf("%d",s[i]);
        }
    }
    int cmp(Hp a) { //比较，类似于char数组比较
        if(len<a.len) return -1;
        else if(len>a.len) return 1;
        else {
            int temp=len-1;
            while(temp>=0&&s[temp]==a.s[temp]) temp--;
            if(temp==-1) return 0;
            else if(s[temp]<a.s[temp]) return -1;
            else return 1;
        }
    }
    void add(Hp a) { //高精度与高精度相加
        int templen=max(len,a.len);
        for(int i=0; i<templen; i++) {
            s[i]+=a.s[i];
            if(s[i]>9) {
                s[i]-=10;
                s[i+1]++;
            }
        }
        len=templen;
        if(s[len]!=0) len++;
    }
    void add(int x) { //高精度加单精度
        int temp=0;
        s[0]+=x;
        while(s[temp]>9) {
            s[temp]-=10;
            temp++;
            s[temp]++;
        }
        if(s[len]!=0) len++;
    }
    void subtract(Hp a) { //高精度减高精度
        for(int i=0; i<len; i++) {
            s[i]-=a.s[i];
            if(s[i]<0) {
                s[i]+=10;
                s[i+1]--;
            }
        }
        while(len>1&&s[len-1]==0) len--;
    }
    void multiply10() { //高精度乘10
        for(int i=len-1; i>=0; i--) {
            s[i+1]=s[i];
        }
        len++;
        if(s[len-1]==0) len--;
    }
    void multiply(int b) { //高精度乘单精度
        int temp=0;
        for(int i=0; i<len; i++) {
            temp+=s[i]*b;
            s[i]=temp%10;
            temp/=10;
        }
        s[len++]=temp;
        while(s[len-1]>10) {
            s[len]+=s[len-1]/10;
            s[len-1]%=10;
            len++;
        }
        while(len>1&&s[len-1]==0) len--;
    }
    void multiply(Hp b) { //高精度乘高精度
        Hp c;
        mt(c.s,0);
        for(int i=0; i<len; i++) {
            for(int j=0; j<b.len; j++) {
                c.s[i+j]+=s[i]*b.s[j];
                c.s[i+j+1]+=c.s[i+j]/10;
                c.s[i+j]%=10;
            }
        }
        len=len+b.len;
        for(int i=0; i<len; i++) {
            s[i]=c.s[i];
        }
        while(len>1&&s[len-1]==0) len--;
    }
    int divide(int b) { //高精度除单精度
        int ret=0;
        for(int i=len-1; i>0; i--) {
            s[i-1]+=(s[i]%b)*10;
            s[i]/=b;
        }
        ret=s[0]%b;
        s[0]/=b;
        while(len>1&&s[len-1]==0) len--;
        return ret;
    }
    Hp divide(Hp b) { //高精度除高精度
        Hp c,d;
        mt(c.s,0);
        mt(d.s,0);
        d.len=1;
        for(int i=len-1; i>=0; i--) {
            d.multiply10();
            d.s[0]=s[i];
            while(d.cmp(b)>=0) {
                d.subtract(b);
                c.s[i]++;
            }
        }
        for(int i=0; i<len; i++) {
            s[i]=c.s[i];
        }
        while(len>1&&s[len-1]==0) len--;
        while(d.len>1&&d.s[d.len-1]==0) d.len--;
        return d;
    }
} dp[8];
int main() {
    int n;
    while(~scanf("%d",&n)){
        if(n<5){
            puts("1");
            continue;
        }
        dp[1].init(1);
        dp[2].init(1);
        dp[3].init(1);
        dp[4].init(1);
        for(int i=5;i<=n;i++){
            dp[5].init(dp[4]);
            for(int j=1;j<=3;j++){
                dp[5].add(dp[j]);
            }
            for(int j=1;j<=4;j++){
                dp[j].init(dp[j+1]);
            }
        }
        dp[5].print();
        puts("");
    }
    return 0;
}

end