poj 3378 Crazy Thairs 树状数组+高精度+dp

题意:

给一个长为n的序列,求里面长度为5的上升子序列有多少个。

分析:

树状数组c[i][j]表示以i结尾长度为j的序列数量,要用高精度。

代码:

//poj 3378
//sep9
#include <iostream>
#include <algorithm>  
using namespace std;  
const int maxN=50012;    
pair<int,int> x[maxN];
int v[maxN];
int n; 
struct H
{
	int a[10],len;
	H(){memset(a,0,sizeof(a));len=1;};
	void clear(){
		memset(a,0,sizeof(a));len=1;
	}
	void print(){
		printf("%d",a[len]);
		for(int i=len-1;i>=1;--i)
			printf("%04d",a[i]);				
		printf("\n");
	}		
};
H add(H x,H y)
{
	H z;
	int i,r=0,len=max(x.len,y.len);
	z.len=len;
	for(i=1;i<=len;++i)	
	{
		z.a[i]=(x.a[i]+y.a[i]+r)%10000;
		r=(x.a[i]+y.a[i]+r)/10000;
	}
	if(r!=0)
		z.a[++z.len]=r;
	return z;
}

struct BIT
{
	H c[maxN][6];
	int n;
	BIT(){}
	void clear(){
		for(int i=1;i<=n;++i)
			for(int j=1;j<5;++j)
				c[i][j].clear();
	}
	int lowbit(int x)
	{
		return x&(x^(x-1));
	}
	void modify(int i,int j,H d)
	{
		while(i<=n){
			c[i][j]=add(c[i][j],d);
			i+=lowbit(i);
		}
	}
	H q(int i,int j)
	{
		H sum;
		for(;i>0;i-=lowbit(i))
			sum=add(sum,c[i][j]);
		return sum;
	}
}bit;

void solve()  
{  
    int i,j;
    for(i=1;i<=n;++i){
        scanf("%d",&x[i].first);  
    	x[i].second=i;
	}
	sort(x+1,x+1+n);
	int t=0;
	for(i=1;i<=n;++i)
		if(i==1||x[i].first!=x[i-1].first)
			v[x[i].second]=++t;
		else
			v[x[i].second]=t;
	bit.n=n;
	bit.clear();	
	H ans;
	for(i=1;i<=n;++i){
		ans=add(ans,bit.q(v[i]-1,4));	
		for(j=4;j>=2;--j){
			H tmp;
			tmp=bit.q(v[i]-1,j-1);
			bit.modify(v[i],j,tmp); 
		}
		H one;
		one.a[1]=1;
		one.len=1;
		bit.modify(v[i],1,one);
	}
	ans.print();
}   

int main()
{

	while(scanf("%d",&n)==1)
		solve();
	return 0;	
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值