test 8 Problem B: [noip2016十连测第八场]降雷皇 （dp+线段树）

Problem B: [noip2016十连测第八场]降雷皇

Time Limit: 10 Sec   Memory Limit: 512 MB
Submit: 29   Solved: 13
[ Submit][ Status][ Web Board]

Description

http://www.lydsy.com/JudgeOnline/upload/201610/test8888.rar

降雷皇哈蒙很喜欢雷电，他想找到神奇的电光。哈蒙有n条导线排成一排，每条导线有一个电阻值，神奇的电光只
能从一根导线传到电阻比它大的上面，而且必须从左边向右传导，当然导线不必是连续的。哈蒙想知道电光最多能
通过多少条导线，还想知道这样的方案有多少。

Input

第一行两个整数n和type。type表示数据类型
第二行n个整数表示电阻。
n<=100000,电阻值不超过100000

Output

第一行一个整数表示电光最多能通过多少条导线。
如果type=1则需要输出第二行，表示方案数，对123456789取模。

Sample Input

5 1
1 3 2 5 4

Sample Output

3
4

HINT

[ Submit][ Status]

题解：dp+线段树。

dp 求最长上升子序列 f[i]=max{f[j]}+1 (j<i,h[j]<h[i]) 如果暴力dp的话是O(N^2)的，但是我们可以用线段树优化dp。现将权值离散化，建一颗以权值为下标的线段树，每次查询[1,h[i]-1]这段区间的最大值。

第二问的话，其实还是可以利用线段树。

g[i]表示以i结尾的最长公共子序列的方案数。

g[i]=sigma g[j]  f[i]=f[j]+1,j<i,h[j]<h[i]

我们可以以f[i]为第一关键字，位置i为第二关键字。

然后统计f[i]=x的时候，设两个指针分别指在f[i]==x-1的开头，和f[i]==x的开头，不断的往后扫，保证计算答案的时候位置在他之前的都已经加入权值线段树。

每次用完打lazy标记，清线段树即可。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
#define N 100003
#define p 123456789
using namespace std;
int a1[N],b[N],c[N],cnt,n,g[N],p1[N],f[N];
struct data
{
	int h,pos,x;
}a[N];
int tr[N*4],delta[N*4],sum[N*4];
int cmp(int x,int y)
{
	return a1[x]<a1[y];
}
int cmp1(data a,data b)
{
	return a.x<b.x||a.x==b.x&&a.pos<b.pos;
}
void update(int now)
{
	tr[now]=max(tr[now<<1],tr[now<<1|1]);
}
void change(int now,int l,int r,int x,int v)
{
	if (l==r) {
		tr[now]=v;
		return;
	}
	int mid=(l+r)/2;
	if (x<=mid) change(now<<1,l,mid,x,v);
	else change(now<<1|1,mid+1,r,x,v);
	update(now);
}
int qjmax(int now,int l,int r,int ll,int rr)
{
	if (ll>rr) return 0;
	if (ll<=l&&r<=rr) return tr[now];
	int mid=(l+r)/2; int ans=0;
	if (ll<=mid) ans=max(ans,qjmax(now<<1,l,mid,ll,rr));
	if (rr>mid) ans=max(ans,qjmax(now<<1|1,mid+1,r,ll,rr));
	return ans;
}
void update1(int now)
{
	sum[now]=(sum[now<<1]+sum[now<<1|1])%p;
}
void pushdown(int now)
{
	if (delta[now]) {
		delta[now]=0;
		delta[now<<1]=1; delta[now<<1|1]=1;
		sum[now<<1]=0; sum[now<<1|1]=0;
	}
}
void pointchange(int now,int l,int r,int x,int v)
{
	if (l==r) {
		sum[now]=(sum[now]+v)%p;
		return ;
	}
	int mid=(l+r)/2;
	pushdown(now);
	if (x<=mid) pointchange(now<<1,l,mid,x,v);
	else pointchange(now<<1|1,mid+1,r,x,v);
	update1(now);
}
int qjsum(int now,int l,int r,int ll,int rr)
{
	if (ll>rr) return 0;
	if (ll<=l&&r<=rr) return sum[now]%p;
	int mid=(l+r)/2;
	pushdown(now); int ans=0;
	if (ll<=mid) ans=(ans+qjsum(now<<1,l,mid,ll,rr))%p;
	if (rr>mid) ans=(ans+qjsum(now<<1|1,mid+1,r,ll,rr))%p;
	return ans%p;
}
int main()
{
	freopen("hamon.in","r",stdin);
	freopen("my.out","w",stdout);
	int t;
	scanf("%d%d",&n,&t);
	for (int i=1;i<=n;i++) scanf("%d",&a1[i]),b[i]=i;
	sort(b+1,b+n+1,cmp);
	for (int i=1;i<=n;i++)
	 if(a1[b[i]]!=a1[b[i-1]])
	  c[b[i]]=++cnt;
	 else c[b[i]]=cnt;
	change(1,1,n,c[1],1);
	f[1]=1;
	for (int i=2;i<=n;i++) {
		f[i]=qjmax(1,1,n,1,c[i]-1)+1;
		change(1,1,n,c[i],f[i]);
    }
	for (int i=1;i<=n;i++) a[i].x=f[i],a[i].h=c[i],a[i].pos=i;
	sort(a+1,a+n+1,cmp1);
	if (t!=1) {
		printf("%d\n",a[n].x);
		return 0;
	}
	int i=1; 
	int ans=0;
	while (a[i].x==1) 
	{
	  g[a[i].pos]=1;
	  if (a[i].x==a[n].x) ans+=g[a[i].pos];
	  i++;
    }
	cnt=0;
	for (int i=1;i<=n;i++)
	 if (a[i].x!=a[i-1].x)  p1[++cnt]=i;
	for (int i=2;i<=cnt;i++)
	 {
	 	delta[1]=1; sum[1]=0;
	 	int j=p1[i-1]; int k=p1[i];
	 	while (a[k].x==i) {
	 		while (a[j].pos<a[k].pos&&j<=p1[i]-1)  pointchange(1,1,n,a[j].h,g[a[j].pos]),j++;
	 		g[a[k].pos]=qjsum(1,1,n,1,a[k].h-1);
	 		if (a[k].x==a[n].x) ans=(ans+g[a[k].pos])%p;
	 		k++;
		 }
	 }
	printf("%d\n%d\n",a[n].x,ans);
}