MyException - 我的异常网
当前位置:我的异常网» 编程 » hdu 4893 Wow! Such Sequence

hdu 4893 Wow! Such Sequence

www.MyException.Cn  网友分享于:2014-08-05  浏览:0次
hdu 4893 Wow! Such Sequence!

http://acm.hdu.edu.cn/showproblem.php?pid=4893


三种操作:

1 k d - "add"
2 l r - "query sum"
3 l r - "change to nearest Fibonacci"


节点附件三个值:

s1:由lazy控制的区间的正确的和。

s2:区间内与所有数相近的fib数之和,随着单点更新而更新。

col:lazy,标记区间是否全部取fib数,是取1,否则取0。

询问区间的和时,找到相应区间直接返回s1,若有col为1的区间要先向下推送,表示要取该区间的fib数的和。

Mark.....

#include <stdio.h>
#include <iostream>
#include <map>
#include <set>
#include <stack>
#include <vector>
#include <math.h>
#include <string.h>
#include <queue>
#include <string>
#include <stdlib.h>
#include <algorithm>
#define LL long long
#define _LL __int64
#define eps 1e-12
#define PI acos(-1.0)
#define C 240
#define S 20
using namespace std;

const int maxn = 100010;
const int INF = 1 << 30;

 _LL fib[86];
void init()
{
    fib[0] = fib[1] = 1;
    for(int i = 2; i <= 85; i++)
        fib[i] = fib[i-1] + fib[i-2];
}

_LL Binsearch(_LL key)
{
    int low = 0,high = 85;
    while(low <= high)
    {
        int mid = (low + high) >> 1;
        if(fib[mid] > key)
            high = mid-1;
        else low = mid+1;
    }
    if(fabs(fib[low]-key) < fabs(fib[high]-key))
        return fib[low];
    else return fib[high];
    /*_LL Min = abs(fib[0]-key);
    int pos = 0;
    for(int i = 1; i <= 85; i++)
    {
        if(Min > abs(fib[i]-key))
        {
            Min = abs(fib[i]-key);
            pos = i;
        }
    }
    return fib[pos];*/
}

struct node
{
    int l, r;
    LL s1,s2;
    int col;
}tree[maxn*4];

void build(int v, int l, int r)
{
    tree[v].l = l;
    tree[v].r = r;
    tree[v].col = 0;
    tree[v].s1 = 0;
    tree[v].s2 = r-l+1;
    if(l == r)
    {
        return;
    }
    int mid = (l+r)>>1;
    build(v*2,l,mid);
    build(v*2+1,mid+1,r);
}

void update_fib(int v, int l, int r)
{
    if(tree[v].l == l && tree[v].r == r)
    {
        tree[v].col = 1;
        tree[v].s1 = tree[v].s2;
        return;
    }

    if(tree[v].col == 1)
        return;
    int mid = (tree[v].l + tree[v].r) >> 1;
    if(r <= mid)
        update_fib(v*2,l,r);
    else if(l > mid)
        update_fib(v*2+1,l,r);
    else
    {
        update_fib(v*2,l,mid);
        update_fib(v*2+1,mid+1,r);
    }
    tree[v].s1 = tree[v*2].s1 + tree[v*2+1].s1;
    tree[v].s2 = tree[v*2].s2 + tree[v*2+1].s2;
}

void update(int v, int pos, _LL d)
{
      if(tree[v].l == tree[v].r)
      {
          tree[v].s1 += d;
          tree[v].s2 = Binsearch(tree[v].s1);
          return;
      }

      if(tree[v].col == 1)
      {
          int mid = (tree[v].l + tree[v].r) >> 1;
          update_fib(v*2,tree[v].l,mid);
          update_fib(v*2+1,mid+1,tree[v].r);
          tree[v].col = 0;
      }

      int mid = (tree[v].l + tree[v].r) >> 1 ;
      if(pos <= mid)
        update(v*2,pos,d);
      else update(v*2+1,pos,d);
      tree[v].s1 = tree[v*2].s1 + tree[v*2+1].s1;
      tree[v].s2 = tree[v*2].s2 + tree[v*2+1].s2;
}

_LL query(int v, int l, int r)
{
    if(tree[v].l == l && tree[v].r == r)
    {
        return tree[v].s1;
    }

    if(tree[v].col == 1)
    {
        int mid = (tree[v].l+tree[v].r)>>1;
        update_fib(v*2,tree[v].l,mid);
        update_fib(v*2+1,mid+1,tree[v].r);
        tree[v].col = 0;
    }
    int mid = (tree[v].l + tree[v].r) >> 1;
    if(r <= mid)
        return query(v*2,l,r);
    else if(l > mid)
        return query(v*2+1,l,r);
    else return query(v*2,l,mid) + query(v*2+1,mid+1,r);
}

int main()
{
    init();
    int n,m;
    int num,l,r,pos;
    _LL d;

    while(~scanf("%d %d",&n,&m))
    {
        build(1,1,n);
        while(m--)
        {
            scanf("%d",&num);
            if(num == 1)
            {
                scanf("%d %I64d",&pos,&d);
                update(1,pos,d);
               // for(int i = 1; i <= 3; i++)
                 //  printf("l : %2d r : %2d col : %2d s1 : %I64d s2 : %I64d\n",tree[i].l,tree[i].r,tree[i].col,tree[i].s1,tree[i].s2);
            }
            else if(num == 2)
            {
                scanf("%d %d",&l,&r);
                printf("%I64d\n",query(1,l,r));
            }
            else
            {
                scanf("%d %d",&l,&r);
                update_fib(1,l,r);
                // for(int i = 1; i <= 3; i++)
                   // printf("l : %2d r : %2d col : %2d s1 : %I64d s2 : %I64d\n",tree[i].l,tree[i].r,tree[i].col,tree[i].s1,tree[i].s2);
            }
        }
    }
    return 0;
}


文章评论

5款最佳正则表达式编辑调试器
5款最佳正则表达式编辑调试器
程序猿的崛起——Growth Hacker
程序猿的崛起——Growth Hacker
程序员都该阅读的书
程序员都该阅读的书
旅行,写作,编程
旅行,写作,编程
程序员最害怕的5件事 你中招了吗?
程序员最害怕的5件事 你中招了吗?
当下全球最炙手可热的八位少年创业者
当下全球最炙手可热的八位少年创业者
 程序员的样子
程序员的样子
要嫁就嫁程序猿—钱多话少死的早
要嫁就嫁程序猿—钱多话少死的早
科技史上最臭名昭著的13大罪犯
科技史上最臭名昭著的13大罪犯
老程序员的下场
老程序员的下场
程序员眼里IE浏览器是什么样的
程序员眼里IE浏览器是什么样的
“肮脏的”IT工作排行榜
“肮脏的”IT工作排行榜
总结2014中国互联网十大段子
总结2014中国互联网十大段子
为什么程序员都是夜猫子
为什么程序员都是夜猫子
我是如何打败拖延症的
我是如何打败拖延症的
程序员的鄙视链
程序员的鄙视链
程序员应该关注的一些事儿
程序员应该关注的一些事儿
十大编程算法助程序员走上高手之路
十大编程算法助程序员走上高手之路
漫画:程序员的工作
漫画:程序员的工作
团队中“技术大拿”并非越多越好
团队中“技术大拿”并非越多越好
那些争议最大的编程观点
那些争议最大的编程观点
我的丈夫是个程序员
我的丈夫是个程序员
聊聊HTTPS和SSL/TLS协议
聊聊HTTPS和SSL/TLS协议
程序员周末都喜欢做什么?
程序员周末都喜欢做什么?
如何区分一个程序员是“老手“还是“新手“?
如何区分一个程序员是“老手“还是“新手“?
程序员必看的十大电影
程序员必看的十大电影
Web开发者需具备的8个好习惯
Web开发者需具备的8个好习惯
代码女神横空出世
代码女神横空出世
亲爱的项目经理,我恨你
亲爱的项目经理,我恨你
程序员的一天:一寸光阴一寸金
程序员的一天:一寸光阴一寸金
中美印日四国程序员比较
中美印日四国程序员比较
老美怎么看待阿里赴美上市
老美怎么看待阿里赴美上市
什么才是优秀的用户界面设计
什么才是优秀的用户界面设计
初级 vs 高级开发者 哪个性价比更高?
初级 vs 高级开发者 哪个性价比更高?
Java 与 .NET 的平台发展之争
Java 与 .NET 的平台发展之争
如何成为一名黑客
如何成为一名黑客
60个开发者不容错过的免费资源库
60个开发者不容错过的免费资源库
写给自己也写给你 自己到底该何去何从
写给自己也写给你 自己到底该何去何从
鲜为人知的编程真相
鲜为人知的编程真相
做程序猿的老婆应该注意的一些事情
做程序猿的老婆应该注意的一些事情
不懂技术不要对懂技术的人说这很容易实现
不懂技术不要对懂技术的人说这很容易实现
2013年中国软件开发者薪资调查报告
2013年中国软件开发者薪资调查报告
看13位CEO、创始人和高管如何提高工作效率
看13位CEO、创始人和高管如何提高工作效率
Web开发人员为什么越来越懒了?
Web开发人员为什么越来越懒了?
10个帮程序员减压放松的网站
10个帮程序员减压放松的网站
一个程序员的时间管理
一个程序员的时间管理
我跳槽是因为他们的显示器更大
我跳槽是因为他们的显示器更大
“懒”出效率是程序员的美德
“懒”出效率是程序员的美德
编程语言是女人
编程语言是女人
软件开发程序错误异常ExceptionCopyright © 2009-2015 MyException 版权所有