题目描述
“低价购买”这条建议是在奶牛股票市场取得成功的一半规则。要想被认为是伟大的投资者,你必须遵循以下的问题建议:“低价购买;再低价购买”。每次你购买一支股票,你必须用低于你上次购买它的价格购买它。买的次数越多越好!你的目标是在遵循以上建议的前提下,求你最多能购买股票的次数。你将被给出一段时间内一支股票每天的出售价(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
*/