蒜头君爬楼梯（1） - 计蒜客

目录

题目 

题解

---

题目 

•  1000ms
•  32767K

 

蒜头君很喜欢爬楼梯，但是蒜头君腿不够长，每次蒜头君最多只能一步跨越两个阶梯。比如他初始在楼底，跨越一个阶梯到达 1 号阶梯，或者跨越两个阶梯到达 2 号阶梯。如下图

为了选出一种最轻松的爬楼梯的方式，蒜头君想把所有不同的到达楼顶的方式都尝试一遍。对于一共有 n 个阶梯的楼梯，蒜头君一共有多少总方法从楼底到达楼顶。

由于最后答案可能很大，输出最后的答案对 100007 取模的结果。

输入格式

第一行输入一个整数 n(1≤n≤1000)。

输出格式

输出最后答案对于 100007 取模的结果。

输出时每行末尾的多余空格，不影响答案正确性

要求使用「文件输入输出」的方式解题，输入文件为 stairs.in，输出文件为 stairs.out

样例输入

4

样例输出

5

 

 

题解：

知识点：递推的应用

分析：其实本题主要是找规律，其次如果学过小学奥数的可以用一个法子来做。这道题比较容易得出是求斐波拉契数列，所以只需要递推即可。

代码：

#include<iostream>
#include<cstdio>
using namespace std;
const int mod=100007;
int main(){
    freopen("stairs.in","r",stdin);
    freopen("stairs.out","w",stdout);
    int n;
    cin>>n;
    long long a=1,b=1,c=0;//用long long防止数据溢出
    if (n==1){//先搞个判断，以免影响接下来的处理
        cout<<1<<endl;//输出当只有一级台阶时的答案
        return 0;
    }
    for (int i=2;i<=n;i++){//用斐波拉契数列来做
        c=a+b;
        a=b;
        b=c;
        a%=mod;//不要忘了取余
        b%=mod;//不要忘了取余
        c%=mod;//不要忘了取余
    }
    cout<<c<<endl;
    return 0;
}

拓展（若能力有限不必掌握，这里本人不用dfs求解）

方法二：

知识点：记忆化搜索的应用

分析：其实本题主要是找规律，其次如果学过小学奥数的可以用一个法子来做。这道题比较容易得出是求斐波拉契数列，所以只需要记忆化搜索即可。

什么是记忆化搜索？记忆化搜索 是level4的内容，是一种搜索的形式，对搜索的结果用数组或其它数据结构记录下来。若当前状态搜索过，则返回已存储的答案。这样，每个状态至多计算一次，与DP（动态规划）差不多。

代码：

#include<iostream>
#include<cstdio>
using namespace std;
const int mod=100007;
int f[10000];
int fib(int x) {//记忆化搜索
    if (f[x]){
        return f[x];//返回已记录的状态
    }
    if (x <= 2) {
        return f[x]=1;//记下状态
    }
    else {
        return f[x]=(fib(x - 1) + fib(x - 2))%mod;//同上，莫忘取余
    }
}

int main() {
    freopen("stairs.in","r",stdin);
    freopen("stairs.out","w",stdout);
    int n;
    cin >> n;
    cout << fib(n+1) << endl;
    return 0;
}

方法三：

知识点：矩阵快速幂的应用

分析：其实本题主要是找规律，其次如果学过小学奥数的可以用一个法子来做。这道题比较容易得出是求斐波拉契数列，所以只需要矩阵快速幂求解即可。

什么是记忆化搜索？矩阵快速幂 是level4的内容，就是把结果初始化为 1 的部分改为了初始化为单位矩阵，乘法部分换成了矩阵乘法，参数和返回值相应换成了矩阵而已。

代码：

/*
因为Fib(n)至于最近的俩个序列有关(及Fib(n-1)和Fib(n-2))，所以我们保存最近的那俩个就行了。
设f(n)表示一个1*2的矩阵，f(n)=[Fib(n),Fib(n+1)]，可以看成【a,b】–>【a+b,b】;
所以可以变成f(n)=f(n-1)*A; (A表示一个二维矩阵) A[2][2]={{0,1},{1,1}};
然后就可以得到最终的表达式 f(n)={0,1}*A^n(表示矩阵A的N次方)；
难点也就在于如何求矩阵A的n次方。这就需要结合快速幂的方法求矩阵的幂次。
为了方便起见，这里我用c语言风格来书写
*/
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int MOD=10000;
void mul(int f[2] ,int a[2][2]){
	int c[2];
	memset(c,0,sizeof(c));
	for(int i=0;i<2;i++)
		for(int j=0;j<2;j++)
			c[i]=(c[i]+(long long)f[j]*a[j][i])%MOD;
	memcpy(f,c,sizeof(c));
}
void mulself(int a[2][2]){
	int c[2][2];
	memset(c,0,sizeof(c));
	for(int i=0;i<2;i++)
		for(int j=0;j<2;j++)
			for(int k=0;k<2;k++)
				c[i][j]=(c[i][j]+(long long)a[i][k]*a[k][j])%MOD;
	memcpy(a,c,sizeof(c));
}
int main(){
	int n;
	while(~scanf("%d",&n)){
		if(n==-1) break;
		int f[2]={0,1};
		int a[2][2]={{0,1},{1,1}};
		while(n){
			if(n&1) mul(f,a);
			mulself(a);
			n>>=1;
		}
		printf("%d\n",f[0]);
	}
	return 0;
}