SJTU OnlineJudge 1012 使用DP的一道题

SJTU OnlineJudge 1012 使用DP的一道题

1. 问题描述
   1.1 Description
   有一个数列，它是由自然数组成的，并且严格单调上升。最小的数不小于S，最大的不超过T。现在知道这个数列有一个性质：后一个数相对于前一个数的增长率总是百分比下的整数（如5相对于4的增长率是25%，25为整数；而9对7就不行了）。现在问：这个数列最长可以有多长？满足最长要求的数列有多少个？
   1.2 Input Format
   输入仅有一行，包含S和T两个数（ $0<S<T\leq 200000$）
   30%的数据，$0<S<T\leq 100$
   100%的数据，$0<S<T\leq 200000$
   1.3 Output Format
   输出有2行。第一行包含一个数表示长度，第二行包含一个数表示个数。
   1.4 Sample Input
   2
   10
   1.5 Sample Output
   5
   2
   1.6 样例解释
   2 4 5 6 9以及2 4 5 8 10

2. 思路
   看到林学长的思路，思考从S到T，考虑如何将问题变为从S到T-1的问题

   这里考虑的状态为以i（$S\leq i\leq T$）结尾的满足条件的最长序列当考虑其下一个数字的时候，如果原有的以下一个数字结尾的序列要比当前序列加上这个数字要长（这种情况在从S推到T的时候并不存在，因为下一个数字被我们考虑的时候总是晚于当前的i的，而我们初始化以每个数字结尾的最长序列都是出现一次的单元素序列（即数字本身）），或者两者相等，那么以下一个数字结尾的最长序列的个数要加上以当前i结尾的最长序列的个数，又或者下一个数字原有的最长序列要比当前的以i结尾的序列要短，那么以下一个数字结尾的最长序列的长度就会变成以当前i结尾的最长序列的长度+1，而出现次数则会变成当前以i结尾的最长序列的出现次数。

   那么我们如何寻找满足符合条件的下一个数字呢？根据题意，这个数字必须不能超过T，而且必须符合对当前i的增长率的百分数为整数，那么设增长率为j% ，当前数为i，下一个数字为tmp，我们有
   

   $$\begin{equation} \dfrac{tmp-i}{i}=\dfrac{j}{100} \end{equation}$$

   整理得到

$$\begin{equation} tmp = i+\dfrac{i*j}{100} \end{equation}$$
那么只要从令i从S遍历到T，在每个loop中把tmp加到序列上面去。
设
d[i]为以i结尾的最长序列的长度
number[i]为以i结尾的最长序列的出现次数
cnt[k]为长度为k的满足题设条件的序列的个数
ans为从S到T中最长的满足条件的序列的长度可以得出状态转移方程

$$\begin{cases} d[i]+1 & d[i]+1 \geq d[tmp]\\ d[tmp] & d[i]+1 = d[tmp] \end{cases}$$

$$\begin{cases} number[i] & d[i]+1 \geq d[tmp]\\ number[i]+number[tmp] & d[i]+1 = d[tmp] \end{cases}$$

$$\begin{equation} ans = max ( ans,d[tmp] ) \end{equation}$$
每一次循环
$$\begin{equation} cnt[d[i]+1]+=number[i] \end{equation}$$
3. 实现

#include<iostream>
#include<iostream>
#include<cstring>
#include<algorithm>
#define MAXSIZE 200010
using namespace std;

int d[MAXSIZE];
long long cnt[MAXSIZE],number[MAXSIZE];

int main()
{
    int S,T;
    cin>>S>>T;

    memset(cnt,0,sizeof(cnt));
    int i,j,tmp,ans=1;
    cnt[1] = T-S+1;

    for(i=S;i<=T;i++)
        d[i]=number[i]=1;

    for(int i=S;i<=T;i++)
        for(j=1;j<=100;j++)
            if((i*j)%100==0)
            {
                tmp = i+i*j/100;
                if(tmp<=T)
                {
                    if(d[i]+1>d[tmp])
                    {
                        d[tmp] = d[i]+1;
                        number[tmp] = number[i];
                    }
                    else if(d[i]+1==d[tmp])
                    {
                        number[tmp]+=number[i];
                    }
                    ans = max(ans,d[tmp]);
                    cnt[d[i]+1]+=number[i];
                }
            }
    cout<<ans<<endl<<cnt[ans];

    return 0;
}