51nod 1376 最长递增子序列的数量(dp+cdq分治)

51nod 1376 最长递增子序列的数量(dp+cdq分治)

数组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个数A[i],表示数组的元素(0 <= A[i] <= 10^9)

Output

输出最长递增子序列的数量Mod 1000000007。
Input示例
5
1
3
2
0
4
Output示例
2

题解:

这题的做法有很多,我就拿来练练cdq分治。
分治的套路:
  第一维排序,第二维分治
  1. 先递归求解区间左半边
  2. 用左边的值更新右边,相当于处理跨越区间中点的情况
  3. 递归求解区间右半边

代码:

#include
using namespace std;
const int maxn=5e5+5;
const int mod=1e9+7; 
int a[maxn];
int idx[maxn];
typedef pair<int,int> P;
P f[maxn];
bool cmp(int x,int y)
{
	return a[x]==a[y]?x>y:a[x]<a[y];
}
void Max(P& x, P y)
{
	if(x.first==y.first)
	{
		x.second=(x.second+y.second)%mod;
	}
	else if(x.first<y.first) x=y;
}
void cdq(int l,int r)
{
	if(l==r) return;
	int m=(l+r)>>1;
	cdq(l,m);
	for(int i=l;i<=r;i++) idx[i]=i;
	sort(idx+l,idx+r+1,cmp);
	P _max(0,0);
	for(int i=l;i<=r;i++)
	{
		int index=idx[i];
//		cout<
		if(index<=m)
		{
			Max(_max,f[index]) ;
		}
		else
		{
			P temp=_max;
			temp.first+=1; 
			Max(f[index],temp);
		}
	 } 
	 cdq(m+1 ,r); 
 } 
int main()
{
	int T,n,m,i,j,k;
	scanf("%d",&n);
	for(i=1;i<=n;i++) scanf("%d",&a[i]);
	for(i=1;i<=n;i++)f[i]=P(1,1) ;
	cdq(1,n);
	P ans(0,0);
	for(int i=1;i<=n;i++) Max(ans,f[i]);
	 printf("%d\n",ans.second);
}

你可能感兴趣的:(解题报告,DP,cdq分治)