51Nod-1376
数组A包含N个整数(可能包含相同的值)。设S为A的子序列且S中的元素是递增的,则S为A的递增子序列。如果S的长度是所有递增子序列中最长的,则称S为A的最长递增子序列(LIS)。A的LIS可能有很多个。例如A为:{1 3 2 0 4},1 3 4,1 2 4均为A的LIS。给出数组A,求A的LIS有多少个。由于数量很大,输出Mod 1000000007的结果即可。相同的数字在不同的位置,算作不同的,例如 {1 1 2} 答案为2。
Input
第1行:1个数N,表示数组的长度。(1 <= N <= 50000) 第2 - N + 1行:每行1个数Ai,表示数组的元素(0 <= Ai <= 10^9)
Output
输出最长递增子序列的数量Mod 1000000007。
题解:本题如果我们按照一般的思维可以用dp的思想,先定义两个数组dp[maxn],num[maxn],dp[i]表示以第i个数结尾的最长上升子序列的长度,num[i]表示以第i个数为结尾的最长上升子序列的个数,显然当我们每次求dp[i]和num[i]的时候都很费时间,计算下来时间复杂度将会是O2,但是我们想到每次我们只需要1到i-1区间的两个最值,我们就可以想到用线段树或树状数组来维护这两个值,每次我们查询的时候要查1到i-1的区间来找最大值,更新的时候我们要更新i到n这段区间的dp[i]与num[i](dp的思想)
#include<bits/stdc++.h>
using namespace std;
const int maxn=50005;
typedef long long LL;
const LL mod=1000000007;
int n;
int a[maxn],b[maxn],c[maxn];
LL sum[maxn];
int lowbits(int x)
{
return x&(-x);
}
void query(int m,int &x,LL &y)
{
x=y=0;
for(int i=m;i>=1;i-=lowbits(i)){
if(c[i]>x){
x=c[i];
y=sum[i];
}
else
if(c[i]==x){
y+=sum[i];
if(y>mod){
y%=mod;
}
}
}
}
void update(int m,int x,LL y)
{
for(int i=m;i<=n;i+=lowbits(i)){
if(c[i]<x){
c[i]=x;
sum[i]=y;
}
else
if(c[i]==x){
sum[i]+=y;
if(sum[i]>mod){
sum[i]%=mod;
}
}
}
}
int main() {
scanf("%d",&n);
for(int i=1;i<=n;i++){
scanf("%d",&a[i]);
b[i]=a[i];
}
sort(b+1,b+1+n);
for(int i=1;i<=n;i++){
a[i]=lower_bound(b+1,b+1+n,a[i])-(b+1);
a[i]++;
}
LL y;
int x;
for(int i=1;i<=n;i++){
query(a[i]-1,x,y);
x++;
//printf("%d %d\n",x,y);
update(a[i],x,x==1?1:y);
}
LL ans;
query(n,x,ans);
printf("%lld\n",ans);
return 0;
}