51Nod - 1376（求最长上升子序列的数量+树状数组维护（线段树也可以））

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;
}