【2014 Multi-University Training Contest 3 1007】/【HDU 4893】 Wow! Such Sequence!

第三场了，还是只有一道题能在赛后补出来。

题目大意：

给出n个数，初始值为0.有三种操作：

1.将某一个数加d。

2.询问区间[l,r]之间的数总和。

3，将区间[l,r]之间的数变成与其距离最近的斐波那契数，若两边距离相同则选择小的。

解题思路：

线段树

对于第一种操作使用点更新，对于第三种操作使用区间更新。

居然上首页了o(∩_∩)o ！！

下面是代码：

#include <stdio.h>#include <string.h>#include <iostream>#include <algorithm>#include <vector>#include <queue>#include <set>#include <map>#include <string>#include <math.h>#include <stdlib.h>#define clear(A, X, SIZE) memset(A, X, sizeof(A[0]) * (SIZE))#define clearall(A, X) memset(A, X, sizeof(A))#define max( x, y )  ( ((x) > (y)) ? (x) : (y) )#define min( x, y )  ( ((x) < (y)) ? (x) : (y) )using namespace std;struct node1{    long long nsum,fsum;    int sta,l,r,m; //all none -1 all f 1} node[100000<<2];long long f[100];void build(int l,int r,int tr){    node[tr].nsum=0;    node[tr].fsum=1;    node[tr].sta=-1;    node[tr].l=l;    node[tr].r=r;    node[tr].m=(l+r)>>1;    if(l==r)return;    int m=(l+r)>>1;    build(l,m,tr<<1);    build(m+1,r,tr<<1|1);    node[tr].fsum=node[tr<<1].fsum+node[tr<<1|1].fsum;}void update1(int p,long long num,int tr){    if(node[tr].l==node[tr].r&&node[tr].l==p)    {        if(node[tr].sta==1)node[tr].nsum=node[tr].fsum;        node[tr].nsum+=num;        int l=2,r=90,i;        while(l<r)        {            i=(l+r)>>1;            if(f[i]>node[tr].nsum)r=i;            else l=i+1;        }        i=r;        if(f[i]-node[tr].nsum<node[tr].nsum-f[i-1])node[tr].fsum=f[i];        else node[tr].fsum=f[i-1];        node[tr].sta=-1;        return;    }    if(node[tr].sta==1)    {        node[tr<<1].sta=1;        node[tr<<1|1].sta=1;    }    if(p>=node[tr].l&&p<=node[tr].m)update1(p,num,tr<<1);    else if(node[tr].m<p&&p<=node[tr].r)update1(p,num,tr<<1|1);    node[tr].nsum=node[tr<<1].nsum+node[tr<<1|1].nsum;    node[tr].fsum=node[tr<<1].fsum+node[tr<<1|1].fsum;    if(node[tr<<1].sta+node[tr<<1|1].sta==-2)node[tr].sta=-1;    else if(node[tr<<1].sta+node[tr<<1|1].sta==2)node[tr].sta=1;    else node[tr].sta=0;}void update(int L,int R,int tr){    if(node[tr].sta==1)    {        return ;    }    if(L<=node[tr].l&&node[tr].r<=R)    {        node[tr].sta=1;        return ;    }    if(L<=node[tr].m)update(L,R,tr<<1);    if(node[tr].m<R)update(L,R,tr<<1|1);    if(node[tr<<1].sta+node[tr<<1|1].sta==-2)node[tr].sta=-1;    else if(node[tr<<1].sta+node[tr<<1|1].sta==2)node[tr].sta=1;    else node[tr].sta=0;}long long query(int L,int R,int tr){    if(L<=node[tr].l&&node[tr].r<=R&&node[tr].sta==-1)return node[tr].nsum;    else if(L<=node[tr].l&&node[tr].r<=R&&node[tr].sta==1)return node[tr].fsum;    if(node[tr].sta==1)    {        node[tr<<1].sta=1;        node[tr<<1|1].sta=1;    }    long long ans=0;    if(L<=node[tr].m)ans+=query(L,R,tr<<1);    if(node[tr].m<R)ans+=query(L,R,tr<<1|1);    if(node[tr<<1].sta+node[tr<<1|1].sta==-2)node[tr].sta=-1;    else if(node[tr<<1].sta+node[tr<<1|1].sta==2)node[tr].sta=1;    else node[tr].sta=0;    return ans;}int main(){    int n,m;    f[1]=1;    f[2]=1;    for(int i=3; i<93; i++)    {        f[i]=f[i-1]+f[i-2];    }    int x,l,r;    long long tt;    while(scanf("%d%d",&n,&m)!=EOF)    {        build(1,n,1);        for(int i=0; i<m; i++)        {            scanf("%d",&x);            if(x==1)            {                scanf("%d%I64d",&l,&tt);                update1(l,tt,1);            }            else if(x==2)            {                scanf("%d%d",&l,&r);                printf("%I64d\n",query(l,r,1));            }            else            {                scanf("%d%d",&l,&r);                update(l,r,1);            }        }    }    return 0;}