HDU 4991 Ordered Subsequence DP(套路,BIT优化)

HDU 4991
题意:给出长度为n的序列a,问a有多少个单调递增的子序列并且其长度为m
n<=1e4,m<=100,a[i]<=1e9.问总的方案数?


设dp[i][j]为递增子序列以a[i]结尾,长度为j的个数.
dp[i][j]=dp[i][j]+dp[k][j-1] (a[k]<a[i]) O(m*n^2) 能不能优化?(这种累加和的形式,感觉都出了好几次了)
a[i]<=1e9 先将a[i]离散化一下 
外层枚举j,则算dp[i][j]时 将dp[i-1][j-1]压入到对应的a[i-1]位置上.

则此时BIT中 存的是dp[1~i-1][j-1]的dp值,BIT.query求出a[i]-1位置前的dp值之和即可

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e4+5,M=110,mod=123456789;
int dp[N][M],c[N],pos[N];
int n,m,a[N];
vector<int> b;
int lowbit(int x){return x&-x;}
void update(int x,int val)
{
	for(int i=x;i<N;i+=lowbit(i))
		c[i]=(c[i]+val)%mod;
}
int query(int x)
{
	int res=0;
	for(int i=x;i>0;i-=lowbit(i))
		res=(res+c[i])%mod;
	return res;
}
int main()
{
	while(~scanf("%d%d",&n,&m))
	{
		b.clear();
		memset(dp,0,sizeof(dp));
		for(int i=1;i<=n;i++)
			scanf("%d",&a[i]),b.push_back(a[i]),dp[i][1]=1;
		sort(b.begin(),b.end());
		b.erase(unique(b.begin(),b.end()),b.end());
		for(int i=1;i<=n;i++)
			pos[i]=lower_bound(b.begin(),b.end(),a[i])-b.begin()+1;
		ll ans=0;
		for(int j=2;j<=m;j++)
		{
			memset(c,0,sizeof(c));
			for(int i=1;i<=n;i++)
			{
				if(i>1)
					update(pos[i-1],dp[i-1][j-1]);
				dp[i][j]=query(pos[i]-1);
				if(j==m)
					ans=(ans+dp[i][j])%mod;
			}
		}
		if(m==1)
			ans=n;
		printf("%lld\n",ans);
	}
	return 0;
}