[洛谷luogu]P1108 低价购买

题目描述

“低价购买”这条建议是在奶牛股票市场取得成功的一半规则。要想被认为是伟大的投资者，你必须遵循以下的问题建议:“低价购买；再低价购买”。每次你购买一支股票,你必须用低于你上次购买它的价格购买它。买的次数越多越好!你的目标是在遵循以上建议的前提下，求你最多能购买股票的次数。你将被给出一段时间内一支股票每天的出售价(2^16范围内的正整数)，你可以选择在哪些天购买这支股票。每次购买都必须遵循“低价购买；再低价购买”的原则。写一个程序计算最大购买次数。

这里是某支股票的价格清单：

日期 1 2 3 4 5 6 7 8 9 10 11 12

价格 68 69 54 64 68 64 70 67 78 62 98 87

最优秀的投资者可以购买最多4次股票，可行方案中的一种是：

日期 2 5 6 10

价格 69 68 64 62

输入输出格式

输入格式：
第1行: N (1 <= N <= 5000)，股票发行天数

第2行: N个数，是每天的股票价格。

输出格式：
输出文件仅一行包含两个数:最大购买次数和拥有最大购买次数的方案数(<=2^31)当二种方案“看起来一样”时（就是说它们构成的价格队列一样的时候）,这2种方案被认为是相同的。

输入输出样例

输入样例#1：
12
68 69 54 64 68 64 70 67 78 62 98 87
输出样例#1：
4 2

---

f[i]表示以i结尾最长的长度，g[i]表示以i结尾最长长度的方案数,初值都为1

第一问

• 一个最长下降子序列问题
• 朴素方法都会哈：f[i]=max(f[i],f[j]+1), a[j]>a[i]
• 还可以用nlogn二分查找算法

第二问

• 每次更新完f[i]后，此时f[i]已经为最优解
• 遍历1~i-1 若a[j]>a[i]&&f[j]+1==f[i],g[i]+=g[j]
• 若a[j]==a[i],g[j]变为0，两者有相同方案，且后面的比前面的更优
• 最后输出前 遍历一遍 若f[i]为最长的 第二问答案加上g[i]（有重复的已经为0了 不影响）

第一问朴素方法的代码:

#include <iostream>
#include <cstring>
using namespace std;
long long n,f[5555],g[5555],ans1,ans2,a[5555];
int bo[77177];
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)
        cin>>a[i];
    f[1]=g[1]=1;
    for(int i=1;i<=n;i++)
    {
        f[i]=1;
        for(int j=1;j<i;j++)    if(a[j]>a[i]) f[i]=max(f[j]+1,f[i]);
        if(f[i]==1) g[i]=1;
        for(int j=i-1;j>=1;j--)
        {
            if(a[j]>a[i]&&f[j]+1==f[i])
                    g[i]+=g[j];
            if(a[j]==a[i])    g[j]=0;
        }
        if(f[i]>ans1)    ans1=f[i];
    }
    for(int i=1;i<=n;i++)
        if(f[i]==ans1)
        {
            ans2+=g[i];
            bo[a[i]]=1;
        }
    cout<<ans1<<' '<<ans2<<endl;
}

 /*
12
68 69 54 64 68 64 70 67 78 62 98 87
*/

升级之后 →.→

二分详解：

• 设一个d[i]为长度为i的序列的最后一个数（若有多个长度相同的，d[i]为这些序列中的最大的末项元素），len为最长序列的长度
• 则在查找中，若d[len]>a[i], 则令len2=++len;否则，就说明在长度为0~len的d[x]中一定有一个长度，使得其结尾的元素大于a[i]，并把此长度+1后赋值给len2.

• len2即为f[i]的值。

• 二分：

• 设l=0，r=len ，now为a[i] ，mid=（l+r)/2
• 循环条件先设为l<=r
• 首先，d数组是一个递减的数组。

• 思想：我们要在d数组中找一个刚好大于a[i]的元素并返回其下标（即d[x]>now且x尽量的大（上下界，在x插入此元素后，不破坏原序列）），这样我们就能用此元素更新d[x+1],使得d[i]中的元素都是长度为i的子序列中末项最大的元素

• 考虑一下几种情况：

• 1.d[mid]< now 则说明区间[mid,l]中一定没有，答案在[l,mid-1]中。
• 2.d[mid]==now 同上 也说明答案在[l,mid-1]中。
• 3.d[mid]>now 此时，mid是一个符合条件的下标，然而在[mid+1,r]中会不会存在更优的呢（下标x更大，且d[x]也大于now），因此我们在这里进行一下判断，若d[mid+1]<=now, return now ,否则l=mid+1（[mid+1，r]中有更有解）。

最后return 的应该是l还是r 可以测试一下 在此程序中都可以。

代码：

#include <iostream>
#include <cstring>
using namespace std;
long long n,f[5555],g[5555],ans1,ans2,a[5555];
int bo[77177];
#define INF 2147483647
int d[5555];
int search(int now,int l,int r)
{
    int mid=(l+r)>>1;
    while(l<=r)
    {
        if(d[mid]<=now) r=mid-1;
        else
            if(d[mid]>now)  
            {
                if(d[mid+1]<=now) return mid;
                else l=mid+1;
            }
        mid=(l+r)>>1;
    //  cout<<l<<' '<<mid<<' '<<r<<endl; 
    }
    return l;
}
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)
        cin>>a[i];
    d[0]=INF;
    f[1]=g[1]=1;
    d[1]=a[1];
    ans1=1;
    int len=1;
    for(int i=2;i<=n;i++)
    {
        f[i]=1;
        int len2=1;
        if(d[len]>a[i]) len2=++len;
        else
            len2=search(a[i],0,len)+1;
        d[len2]=a[i];
        f[i]=len2;
        if(f[i]==1) g[i]=1;
        for(int j=i-1;j>=1;j--)
        {
            if(a[j]>a[i]&&f[j]+1==f[i])
                    g[i]+=g[j];
            if(a[j]==a[i])  g[j]=0;
        }
        if(f[i]>ans1)   ans1=f[i];
    //  cout<<f[i]<<' '<<i<<' '<<search(a[i],1,len)<<endl;
    }
    for(int i=1;i<=n;i++)
        if(f[i]==ans1)
        {
            ans2+=g[i];
            bo[a[i]]=1;
        }
    cout<<ans1<<' '<<ans2<<endl;
}   

 /*
12
68 69 54 64 68 64 70 67 78 62 98 87
*/