A Simple Problem with Integers POJ 3468

线段树

抄的博客

#include<cstdio>
#include<algorithm>
#include<iostream>

using namespace std;
const int MAXN=100000;
int num[MAXN];
struct Node
{
    int l,r;//区间的左右端点
long long nSum;//区间上的和
long long Inc;//区间增量的累加 
}segTree[MAXN*3];
void Build(int i,int l,int r)
{
    segTree[i].l=l;
    segTree[i].r=r;
    segTree[i].Inc=0;
    if(l==r)
    {
        segTree[i].nSum=num[l];
        return;
    }    
    int mid=(l+r)>>1;
    Build(i<<1,l,mid);
    Build(i<<1|1,mid+1,r);
    segTree[i].nSum=segTree[i<<1].nSum+segTree[i<<1|1].nSum;
} 
void Add(int i,int a,int b,long long c)//在结点i的区间（a,b）上增加c
{
    if(segTree[i].l==a&&segTree[i].r==b)
    {
        segTree[i].Inc+=c;
        return;
    }
    segTree[i].nSum+=c*(b-a+1);
    int mid=(segTree[i].l+segTree[i].r)>>1;
    if(b<=mid)  Add(i<<1,a,b,c);
    else if(a>mid)  Add(i<<1|1,a,b,c);
    else
    {
        Add(i<<1,a,mid,c);
        Add(i<<1|1,mid+1,b,c);
    }        
}  
long long Query(int i,int a,int b)//查询a-b的总和
{
    if(segTree[i].l==a&&segTree[i].r==b)
    {
        return segTree[i].nSum+(b-a+1)*segTree[i].Inc;
    }  
    segTree[i].nSum+=(segTree[i].r-segTree[i].l+1)*segTree[i].Inc;
    int mid=(segTree[i].l+segTree[i].r)>>1;
    Add(i<<1,segTree[i].l,mid,segTree[i].Inc);
    Add(i<<1|1,mid+1,segTree[i].r,segTree[i].Inc);
    segTree[i].Inc=0;
    if(b<=mid)  return Query(i<<1,a,b);
    else if(a>mid)  return Query(i<<1|1,a,b);
    else return Query(i<<1,a,mid)+Query(i<<1|1,mid+1,b);
}  
int main()
{
    int n,q;
    int i;
    int a,b,c;
    char ch;
    while(scanf("%d%d",&n,&q)!=EOF)
    {
        for(i=1;i<=n;i++)  scanf("%d",&num[i]);
        Build(1,1,n);
        for(i=1;i<=q;i++)
        {
            cin>>ch;
            if(ch=='C')
            {
                scanf("%d%d%d",&a,&b,&c);
                Add(1,a,b,c);
            }    
            else
            {
                scanf("%d%d",&a,&b);
                printf("%I64d\n",Query(1,a,b));
            }    
        }    
    } 
    return 0;   
}

树状数组

原博客

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>

using namespace std;

const int MAXN=101000;

long long a[MAXN<<1];
long long b[MAXN<<1];
long long pos[MAXN<<1];
int n;

int lowbit(int x)
{
	return x&(-x);
}

void update(long long *arry,int x,int d)
{
	while(x<=n)
	{
		arry[x]+=d;
		x+=lowbit(x);
	}
}

long long add(long long *arry,int x)
{
	long long sum=0;
	while(x>0)
	{
		sum+=arry[x];
		x-=lowbit(x);
	}
	return sum;
}

int main()
{
	int k;
	int x,y,z;
	scanf("%d%d",&n,&k);
	memset(a,0,sizeof(a));
	memset(b,0,sizeof(b));
	memset(pos,0,sizeof(pos));
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&x);
		pos[i]+=pos[i-1]+x;
	}
	char c;
	while(k--)
	{
		cin>>c;
		if(c=='C')
		{
			scanf("%d%d%d",&x,&y,&z);
			update(a,x,z);
			update(a,y+1,-z);
			update(b,x,z*x);
			update(b,y+1,-z*(y+1));
		}
		else
		{
			scanf("%d%d",&x,&y);
			long long sum=-pos[x-1]-x*add(a,x-1)+add(b,x-1);
			sum+=pos[y]+(y+1)*add(a,y)-add(b,y);
			printf("%I64d\n",sum);
		}
	}
	return 0;
}