Usaco 4.3.1 Buy Low, Buy Lower 详细解题报告

 

Usaco 4.3.1

Buy Low, Buy Lower

By 小兔齐齐

描述

“逢低吸纳”是炒股的一条成功秘诀。如果你想成为一个成功的投资者，就要遵守这条秘诀:

"逢低吸纳,越低越买"

这句话的意思是：每次你购买股票时的股价一定要比你上次购买时的股价低.按照这个规则购买股票的次数越多越好，看看你最多能按这个规则买几次。

给定连续的N天中每天的股价。你可以在任何一天购买一次股票，但是购买时的股价一定要比你上次购买时的股价低。写一个程序，求出最多能买几次股票。

以下面这个表为例, 某几天的股价是:

天数 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

格式

PROGRAM NAME: buylow

INPUT FORMAT:

(file buylow.in)

第1行: N (1 <= N <=5000), 表示能买股票的天数。

第2行以下: N个正整数 (可能分多行) ，第i个正整数表示第i天的股价. 这些正整数大小不会超过longint(pascal)/long(c++).

OUTPUT FORMAT:

(file buylow.out)

只有一行，输出两个整数：

能够买进股票的天数 长度达到这个值的股票购买方案数量

在计算方案的数量的时候，如果两个方案的股价序列相同，那么这样的两个方案被认为是相同的（只能算做一个方案）。因此，两个不同的天数序列可能产生同一个股价序列，这样只能计算一次。

SAMPLE INPUT

12

68 69 54 64 68 64 70 67

78 62 98 87

SAMPLE OUTPUT

4 2

 

分析：

1）  第一问，最大下降子序列，经典dp，这个就不多说了，不知道的直接去搜索 “最大下降（上升）子序列”

2）  第二问，问最大下降长度的序列的种类，且单词完全相同的不重复计算。这个有点麻烦。

1.       对于这个问题我们还可以开一个数组num[N] ,num[i]记录第i个位置之前，对应dp[i]长度(表示以第i位结尾的最长下降序列长度)的序列种类数。举例：

原始	1	16	17	18	20	10	22	22	8	17	26	14	3	24	8	1	2	21	2	17
dp	1	1	1	1	1	2	1	1	3	2	1	3	4	2	4	5	5	3	5	4
Num	1	1	1	1	1	4	1	1	4	3	1	3	7	1	3	10	10	1	10	1

                    

2.       num[i]如何更新的呢？应该是累加前面满足dp[j]==dp[i]-1 的所有j（j即合法序列的前驱的那一位）的num[j]之和。但是注意一个问题，可能序列是重复的，例如：

9 8 7 6 2 6 5

第一个出现的6和第二个6，对应的递减序列 都是 9 8 7 6 ，属于重复的。此时要记录当前满足dp[i]== dp[j]+1（j即合法序列的前驱的那一位）num[j]是否统计过。

3.现在还有一个问题。就是遇到 前驱位 两个相同的数，是随便取一位么？有讲究么？

当然！例如：

原始	9	7	5	8	5	1
Dp	1	2	3	2	3	4
Num	1	1	1	1	2	2

 

想必大家已经注意到了：对于第一个出现的5和第二个5，他们的dp一样，但是num却不一样。这里我们可以取最后一个5，原因如下：

对于非最后一个5的序列，最后一个5，一定可以取得。例如对于第二个5，第一个5的9 7 5序列，第二个5同样可以取得。而且后面的5可能会有更多的取法，例如上例中的第二个5，还可以获得9 8 5这个序列。所以我们这里，最后1对应的num应该是2。

 

3.       实现方式：鉴于这种情况，我们可以从后往前搜索，并记录visited[],访问过表示累加过了，前面出现相同的就忽略了。可以保证正确性

 

 

3）细节方面，

1.      可以在序列末尾+个0，方便统计如果有多个最大的长度的总情况数。

2.      此题要求高精度.数太大了，我的代码用的是以前写的高精度模板，可能比较臃肿.呵呵

 

 

代码：

//核心代码在main函数里面，前面的大段是高精度模板，请无视

/* ID:dingyag1 LANG:C++ TASK:buylow */ #include <stdio.h> #include <iostream> #include <string> #include <vector> using namespace std; void func_add(vector<int> ∑,vector<int> a,int n) { vector<int>::size_type i; for(i=0;i<a.size();i++) { if(i+n+1>sum.size()) sum.insert(sum.end(),i+n+1-sum.size(),0); sum[i+n]+=a[i]; } // cout<<"输出add函数中sum:"<<endl; // print(sum,4); for(i=0;i<sum.size();i++) { if(sum[i]>=10000) { if(i==sum.size()-1) sum.push_back(0); sum[i+1]+=sum[i]/10000; sum[i]=sum[i]%10000; } } } vector<int> & operator +=(vector<int> &a,const vector<int> &b ) { func_add(a,b,0); return a; } void mitoa(int a,char *s,int c) { int n,i=0; char s2[10]; do { s2[i++]=a%c+'0'; a=a/c; }while(a>0); for(n=0;n<i;n++) s[n]=s2[i-n-1]; s[n]='/0'; } void print(const vector<int> ∑,int nlog) { vector<int>::size_type i; string s; int tmp; bool f=0; char ss[10]={0}; for (i=0;i<sum.size();++i) { if (sum[i]!=0) { f=1; break; } } if (!f) { cout<<0; } else { //取头 mitoa(sum[sum.size()-1],ss,10); s+=ss; if(sum.size()>1){ //取中间 for(i=sum.size()-2;i>0;i--) { if(sum[i]<10) { s+="000"; } else if(sum[i]<100) { s+="00"; } else if(sum[i]<1000) { s+="0"; } mitoa(sum[i],ss,10); s+=ss; } //取尾 if(sum[0]<10) { s+="000"; } else if(sum[0]<100) { s+="00"; } else if(sum[0]<1000) { s+="0"; } mitoa(sum[0],ss,10); s+=ss; } //处理小数点 tmp=s.size(); if(nlog==0) { } else if(nlog>0) { s.append(nlog,'0'); } else if(nlog>-tmp) { s.insert(tmp+nlog,1,'.'); i=s.find_last_not_of('0'); s.erase(i+1,s.size()-i-1); if(s[s.size()-1]=='.') s.erase(s.size()-1,1); } else { s.insert((string::size_type)0,-(tmp+nlog),'0'); s.insert((string::size_type)0,"0."); i=s.find_last_not_of('0'); s.erase(i+1,s.size()-i-1); if(s[s.size()-1]=='.') s.erase(s.size()-1,1); } //删除前面的0 i=s.find_first_not_of('0'); s.erase(0,i); cout<<""<<s; // cin>>s; } } int main() { freopen("buylow.in","r",stdin); freopen("buylow.out","w",stdout); int N,S[5003],dp[5003]={0}; int i,j,max=1,maxp; vector<int> num[5003]; cin>>N; for (i=1;i<=N;++i) { cin>>S[i]; dp[i]=1; //初始化dp为1 num[i].push_back(1); } ///末尾+0 方便统计 N++; S[N]=0;num[N].push_back(1);dp[N]=1; /// if (N==1) { maxp=1; } for (i=2;i<=N;++i) { bool visited[50001]={0}; for (j=i-1;j>0;--j) { if (S[j]>S[i] && dp[j]+1>=dp[i]) { if (dp[j]+1>dp[i]) { dp[i]=dp[j]+1; num[i]=num[j]; if (dp[i]>max) { max=dp[i]; maxp=i; } //record }else if (dp[i]==dp[j]+1) { if (!visited[S[j]]) { num[i]+=num[j]; }else { // num[i]+=num[j]; } } visited[S[j]]=1; } } } if (max==1) { cout<<1<<" "<<N<<endl; return 0; } cout<<max-1<<" "; print(num[maxp],0); cout<<endl; return 0; }