有兴趣可以用矩阵加速做一下。

简单代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int n,f[1000005];

int main(){
    scanf("%d",&n);
    f[1]=1,f[2]=1;
    for(int i=3;i<=n;i++){
        f[i]=(f[i-1]+f[i-2])%10007;
    }
    cout<<f[n]<<endl;
    return 0;
}