数的计数(Noip2001)

【问题描述】
我们要求找出具有下列性质数的个数（包括输入的自然数n）。先输入一个自然数n（n≤1000），然后对此自然数按照如下方法进行处理：
不作任何处理；
在它的左边加上一个自然数，但该自然数不能超过原数的一半；
加上数后，继续按此规则进行处理，直到不能再加自然数为止。
输入：自然数n（n≤1000）
输出：满足条件的数
【输入样例】
6 满足条件的数为 6 (此部分不必输出)
16
26
126
36
136
【输出样例】
6
【方法一】
用递归，f(n)=1+f(1)+f(2)+…+f(div/2)，当n较大时会超时，时间应该为指数级。

【参考程序】
#include<iostream>
using namespace std;
int ans;
void dfs(int m)                           //统计m所扩展出的数据个数
{
    int i;
    ans++;                                  //每出现一个原数，累加器加1；
    for (i = 1; i <= m/2; i++)         //左边添加不超过原数一半的自然数，作为新原数
      dfs(i); 
}
int main()
{
    int n;
    cin >> n;
    dfs(n);
    cout << ans;
    return 0;
} 
【方法二】：用记忆化搜索，实际上是对方法一的改进。设h[i]表示自然数i满足题意三个条件的数的个数。如果用递归求解，会重复来求一些子问题。例如在求h[4]时，需要再求h[1]和h[2]的值。现在我们用h数组记录在记忆求解过程中得出的所有子问题的解，当遇到重叠子问题时，直接使用前面记忆的结果。
【参考程序】
#include<iostream>
using namespace std;
int h[1001];
void dfs(int m)
{
    int i;
    if (h[m] != -1) return;         //说明前面已经求得h[m]的值，直接引用即可，不需要再递归
    h[m] = 1;                           //将h[m]置为1，表示m本身为一种情况
    for (i = 1; i <= m/2; i++)
    {
        dfs(i);
        h[m] += h[i];
    }   
}
int main()
{
    int n;
    cin >> n;
    for (int i = 1; i <= n; i++)
      h[i] = -1;                         //h数组初始化为-1
    dfs(n);                              //由顶到下记忆化递归求解
    cout << h[n];   
    return 0;
}
方法三】
           用递推，用h(n)表示自然数n所能扩展的数据个数，则h(1)=1, h(2)=2, h(3)=2, h(4)=4, h(5)=4, h(6)=6, h(7)=6, h(8)=10, h(9)=10.分析以上数据，可得递推公式：h(i)=1+h(1)+h(2)+…+h(i/2)。此算法的时间度为O(n*n)。
设h[i]-i按照规则扩展出的自然数个数(1≤i≤n)。下表列出了h[i]值及其方案：
【参考程序】
#include<iostream>
using namespace std;
int h[10001];
int main()
{
    int n;
    cin >> n;
    for (int i = 1; i <= n; i++)   //按照递增顺序计算扩展出的自然数的个数
    {
        h[i] = 1;                        //扩展出的自然数包括i本身
        for (int j = 1; j <= i/2; j++)         
                                 //i左边分别加上1…自然数 按规则扩展出的自然数
          h[i] += h[j]; 
    }
    cout << h[n];
    return 0;
}
【方法四】
是对方法三的改进，我们定义数组s，s(x)=h(1)+h(2)+…+h(x),h(x)=s(x)-s(x-1)，此算法的时间复杂度可降到O(n)。
【参考程序】
#include<iostream>
using namespace std;
int h[1001],s[1001];
int main()
{
    int n;
    cin >> n;
    for (int i = 1; i <= n; i++)
    {
        h[i] = 1 + s[i/2];
        s[i] = s[i-1] + h[i];                  //s是h的前缀累加和
    }
    cout << h[n];
    return 0;
}
【方法五】
     还是用递推，只要作仔细分析，其实我们还可以得到以下的递推公式: (1)当i为奇数时，h(i)=h(i-1);
     (2)当i为偶数时,h(i)=h(i-1)+h(i/2).
【参考程序】
#include<iostream>
using namespace std;
int h[1001];
int main()
{
    int n;
    cin >> n;
    h[1] = 1;
    for (int i = 2; i <= n; i++)
    {
        h[i] = h[i-1];
        if (i % 2 == 0) h[i] = h[i-1] + h[i/2];
    }
    cout << h[n];
    return 0;
}