hdu 4578 Transformation 2013ACM-ICPC杭州赛区全国邀请赛

最新推荐文章于 2021-02-16 16:46:11 发布
最新推荐文章于 2021-02-16 16:46:11 发布
阅读量400 收藏
点赞数
分类专栏： ACM
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/aksdrxzd/article/details/45602977
版权

TLE 代码（错在询问操作不必更新到叶节点）：

#include<cctype>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<string>
#include<queue>
#include<stack>
#include<set>
#include<map>
#define CLEAR(a) memset((a),0,sizeof((a)))

using namespace std;

typedef long long LL;
const double pi = acos(-1.0);
const int maxn=1e5+10;
const int inf=99999999;
const double eps=1e-3;
const int mod=10007;

struct Node
{
    int l,r,mid;
    LL v,add,mul;
} SegTree[maxn*4];

void build(int rt,int l,int r)
{
    SegTree[rt].l=l;SegTree[rt].r=r;
    SegTree[rt].mid=(l+r)>>1;
    SegTree[rt].add=SegTree[rt].v=0;
    SegTree[rt].mul=1;
    if (l<r)
    {
        build(rt<<1,l,SegTree[rt].mid);
        build(rt<<1|1,SegTree[rt].mid+1,r);
    }
}

void add(int rt,int x,int y,int num)
{
    int &l=SegTree[rt].l,&r=SegTree[rt].r,&mid=SegTree[rt].mid;
    if (x>y) return;
    //cout<<"+ "<<rt<<' '<<l<<' '<<r<<' '<<mid<<"  "<<x<<' '<<y<<' '<<num<<endl;
    if (x<=SegTree[rt].l&&SegTree[rt].r<=y)
    {
        SegTree[rt].add+=num;
        SegTree[rt].add%=mod;
        return;
    }
    int tmp=SegTree[rt].add;
    if (SegTree[rt].add)
    {

        SegTree[rt].add=0;
        add(rt,SegTree[rt].l,x-1,tmp);
        add(rt,y+1,SegTree[rt].r,tmp);
    }
    if (x<=mid) add(rt<<1,x,min(mid,y),tmp+num);
    if (y>mid) add(rt<<1|1,max(mid+1,x),y,tmp+num);
}

void mul(int rt,int x,int y,int num)
{
    int &l=SegTree[rt].l,&r=SegTree[rt].r,&mid=SegTree[rt].mid;
    if (x>y) return;
    if (x<=SegTree[rt].l&&SegTree[rt].r<=y)
    {
        SegTree[rt].mul*=num;SegTree[rt].mul%=mod;
        SegTree[rt].add*=num;SegTree[rt].add%=mod;
        return;
    }
    int tmp=SegTree[rt].mul;
    if (SegTree[rt].mul!=1)
    {
        SegTree[rt].mul=1;
        add(rt,SegTree[rt].l,x-1,SegTree[rt].add*tmp);
        mul(rt,SegTree[rt].l,x-1,tmp);
        add(rt,x,y,SegTree[rt].add*tmp*num);
        add(rt,y+1,SegTree[rt].r,SegTree[rt].add*tmp);
        mul(rt,y+1,SegTree[rt].r,tmp);
    }
    if (x<=mid) mul(rt<<1,x,min(mid,y),tmp*num);
    if (y>mid) mul(rt<<1|1,max(mid+1,x),y,tmp*num);
}

void cover(int rt,int x,int y,int num)
{
    int &l=SegTree[rt].l,&r=SegTree[rt].r,&mid=SegTree[rt].mid;
    if (x>y) return;
    //cout<<"c "<<rt<<' '<<l<<' '<<r<<' '<<mid<<"  "<<x<<' '<<y<<' '<<num<<endl;
    if (x<=SegTree[rt].l&&SegTree[rt].r<=y)
    {
        SegTree[rt].v=num;
        SegTree[rt].add=0;
        SegTree[rt].mul=1;
        return;
    }
    if (SegTree[rt].add)
    {
        int tmp=SegTree[rt].add;
        SegTree[rt].add=0;
        add(rt,SegTree[rt].l,x-1,tmp);
        add(rt,y+1,SegTree[rt].r,tmp);
    }
    if (SegTree[rt].mul!=1)
    {
        int tmp=SegTree[rt].mul;
        SegTree[rt].mul=1;
        mul(rt,SegTree[rt].l,x-1,tmp);
        mul(rt,y+1,SegTree[rt].r,tmp);
    }
    if (x<=mid) cover(rt<<1,x,min(mid,y),num);
    if (y>mid) cover(rt<<1|1,max(mid+1,x),y,num);
}

LL query(int rt,int x,int y,int p)
{
    int &l=SegTree[rt].l,&r=SegTree[rt].r,&mid=SegTree[rt].mid;
    //cout<<"q "<<rt<<' '<<l<<' '<<r<<' '<<mid<<"  "<<SegTree[rt].v<<' '<<SegTree[rt].add<<' '<<SegTree[rt].mul<<endl; //<<"  "<<x<<' '<<y<<' '<<endl;
    if (SegTree[rt].l==SegTree[rt].r)
    //if (SegTree[rt].v)
    {
        SegTree[rt].v*=SegTree[rt].mul;SegTree[rt].v%=mod;
        SegTree[rt].v+=SegTree[rt].add;SegTree[rt].v%=mod;
        SegTree[rt].add=0;
        SegTree[rt].mul=1;
        switch(p)
        {
            case 1:return SegTree[rt].v;
            case 2:return SegTree[rt].v*SegTree[rt].v%mod;
            case 3:return SegTree[rt].v*SegTree[rt].v*SegTree[rt].v%mod;
        }

    }
    if (SegTree[rt].v)
    {
        add(rt<<1,l,mid,SegTree[rt].v*SegTree[rt].mul+SegTree[rt].add);
        add(rt<<1|1,mid+1,r,SegTree[rt].v*SegTree[rt].mul+SegTree[rt].add);
        SegTree[rt].v=0;
    }
    if (SegTree[rt].add)
    {
        //puts("update");
        add(rt<<1,l,mid,SegTree[rt].add);
        add(rt<<1|1,mid+1,r,SegTree[rt].add);
        SegTree[rt].add=0;
    }
    if (SegTree[rt].mul!=1)
    {
        mul(rt<<1,l,mid,SegTree[rt].mul);
        mul(rt<<1|1,mid+1,r,SegTree[rt].mul);
        SegTree[rt].mul=1;
    }
    return (((x<=mid)?query(rt<<1,x,min(mid,y),p):0)+((y>mid)?query(rt<<1|1,max(mid+1,x),y,p):0))%mod;

}

LL fnd(int rt,int x)
{
    if(SegTree[rt].l==SegTree[rt].r) return SegTree[rt].v;
    else if (x<=SegTree[rt].mid) return fnd(rt<<1,x);
    else return fnd(rt<<1|1,x);
}

int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        if (!n&&!m) break;
        build(1,1,n);
        for(int i=1;i<=m;i++)
        {
            int op,x,y,c;
            scanf("%d%d%d%d",&op,&x,&y,&c);
            switch(op)
            {
                case 1:
                    add(1,x,y,c);//query(1,x,y,1);
                    break;
                case 2:
                    mul(1,x,y,c);//query(1,x,y,1);
                    break;
                case 3:
                    cover(1,x,y,c);//query(1,x,y,1);
                    break;
                case 4:
                    printf("%I64d\n",query(1,x,y,c));
                    break;
            }
            //for(int i=1;i<=n;i++) printf("%lld ",fnd(1,i));printf("\n");
        }
    }
    return 0;
}



AC 4882ms:

#include<cctype>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<string>
#include<queue>
#include<stack>
#include<set>
#include<map>
#define CLEAR(a) memset((a),0,sizeof((a)))

using namespace std;

typedef long long LL;
const double pi = acos(-1.0);
const int maxn = 1e5 + 10;
const int inf = 99999999;
const double eps = 1e-3;
const int mod = 10007;

struct Node
{
    int l, r, mid;
    LL v, add, mul;
    bool up;
} SegTree[maxn * 4];

void build(int rt, int l, int r)
{
    SegTree[rt].l = l;
    SegTree[rt].r = r;
    SegTree[rt].mid = (l + r) >> 1;
    SegTree[rt].add = SegTree[rt].v = 0;
    SegTree[rt].mul = 1;
    if (l < r)
    {
        SegTree[rt].up = 0;
        build(rt << 1, l, SegTree[rt].mid);
        build(rt << 1 | 1, SegTree[rt].mid + 1, r);
    }
    else
    {
        SegTree[rt].up = 1;
    }
}

void pushdown(int rt)
{
    int mid = SegTree[rt].mid;
    if (SegTree[rt].up)
    {
        SegTree[rt << 1].add = SegTree[rt << 1 | 1].add = 0;
        SegTree[rt << 1].mul = SegTree[rt << 1 | 1].mul = 1;
        SegTree[rt << 1].v = SegTree[rt << 1 | 1].v = SegTree[rt].v;
        SegTree[rt << 1].up = SegTree[rt << 1 | 1].up = 1;
        SegTree[rt].up = 0;
    }
    else
    {
        if (SegTree[rt].add)
        {
            if (SegTree[rt << 1].up)
            {
                SegTree[rt << 1].v += SegTree[rt].add;
                SegTree[rt << 1].v %= mod;
            }
            else
            {
                pushdown(rt << 1);
                SegTree[rt << 1].add += SegTree[rt].add;
                SegTree[rt << 1].add %= mod;
            }
            if (SegTree[rt << 1 | 1].up)
            {
                SegTree[rt << 1 | 1].v += SegTree[rt].add;
                SegTree[rt << 1 | 1].v %= mod;
            }
            else
            {
                pushdown(rt << 1 | 1);
                SegTree[rt << 1 | 1].add += SegTree[rt].add;
                SegTree[rt << 1 | 1].add %= mod;
            }
            SegTree[rt].add = 0;
        }
        if (SegTree[rt].mul != 1)
        {
            if (SegTree[rt << 1].up)
            {
                SegTree[rt << 1].v *= SegTree[rt].mul;
                SegTree[rt << 1].v %= mod;
            }
            else
            {
                pushdown(rt << 1);
                SegTree[rt << 1].mul *= SegTree[rt].mul;
                SegTree[rt << 1].mul %= mod;
            }
            if (SegTree[rt << 1 | 1].up)
            {
                SegTree[rt << 1 | 1].v *= SegTree[rt].mul;
                SegTree[rt << 1 | 1].v %= mod;
            }
            else
            {
                pushdown(rt << 1 | 1);
                SegTree[rt << 1 | 1].mul *= SegTree[rt].mul;
                SegTree[rt << 1 | 1].mul %= mod;
            }
            SegTree[rt].mul = 1;
        }
    }
}

void update(int rt, int x, int y, LL num,int op)
{
    int& l = SegTree[rt].l, &r = SegTree[rt].r, &mid = SegTree[rt].mid;
    //cout<<"+ "<<rt<<' '<<l<<' '<<r<<' '<<mid<<"  "<<x<<' '<<y<<' '<<num<<endl;
    if (x <= SegTree[rt].l && SegTree[rt].r <= y)
    {
        if (op==3)
        {
            SegTree[rt].add=0;
            SegTree[rt].mul=1;
            SegTree[rt].v=num;
            SegTree[rt].up=1;
        }
        else
        {
            if (SegTree[rt].up)
            {
                if (op==1)
                {
                    SegTree[rt].v+=num;SegTree[rt].v%=mod;
                }
                else
                {
                    SegTree[rt].v*=num;SegTree[rt].v%=mod;
                }
            }
            else
            {
                pushdown(rt);
                if (op==1)
                {
                    SegTree[rt].add+=num;SegTree[rt].add%=mod;
                }
                else
                {
                    SegTree[rt].mul*=num;SegTree[rt].mul%=mod;
                }
            }
        }
        return;
    }
    pushdown(rt);
    if (x <= mid)
        update(rt << 1, x, min(mid, y), num,op);
    if (y > mid)
        update(rt << 1 | 1, max(mid + 1, x), y, num,op);
}

LL query(int rt, int x, int y, int p)
{
    int& l = SegTree[rt].l, &r = SegTree[rt].r, &mid = SegTree[rt].mid;
    //cout<<"q "<<rt<<' '<<l<<' '<<r<<' '<<mid<<"  "<<SegTree[rt].v<<' '<<SegTree[rt].add<<' '<<SegTree[rt].mul<<endl; //<<"  "<<x<<' '<<y<<' '<<endl;
    if (x <= SegTree[rt].l && SegTree[rt].r <= y&&SegTree[rt].up)
    {
        LL tmp=LL(r-l+1)%mod;
        for(int i=1;i<=p;i++)
        {
            tmp*=SegTree[rt].v;tmp%=mod;
        }
        return tmp;
    }
    pushdown(rt);
    return (((x <= mid) ? query(rt << 1, x, min(mid, y), p) : 0LL) + ((y > mid) ? query(rt << 1 | 1, max(mid + 1, x), y, p) : 0LL)) % mod;
}

LL fnd(int rt, int x)
{
    if (SegTree[rt].l == SegTree[rt].r)
    {
        return SegTree[rt].v;
    }
    else if (x <= SegTree[rt].mid)
    {
        return fnd(rt << 1, x);
    }
    else
    {
        return fnd(rt << 1 | 1, x);
    }
}

int main()
{
    int n, m;
    while (~scanf("%d%d", &n, &m))
    {
        if (!n && !m)
        {
            break;
        }
        build(1, 1, n);
        for (int i = 1; i <= m; i++)
        {
            int op, x, y;
            LL c;
            scanf("%d%d%d%I64d", &op, &x, &y, &c);
            switch (op)
            {
                case 1:
                case 2:
                case 3:
                    update(1, x, y, c,op); //query(1,x,y,1);
                    break;
                case 4:
                    printf("%I64d\n", query(1, x, y, c));
                    break;
            }
            //for(int i=1;i<=n;i++) printf("%lld ",fnd(1,i));printf("\n");
        }
    }
    return 0;
}


aksdrxzd
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
[2018 ACM-ICPC 焦作赛区网络赛] G - Give Candies（找规律+快速幂）
04-18
思路：这个题和HDU - 5703类似，其实就是拆数问题，一个数的拆法其实就是2^(N-1)，具体证明过程可以直接搜刚才杭电那道题的题解，所以这道题其实就是让你算2^(N-1)，但是题目给的N特别大， 可以达到10^100000，我们...
ACM-ICPC-OJ-Code:流行的在线裁判系统的一些解决方案代码
06-19
ACM-ICPC OJ 解决方案代码 介绍 我的解决方案代码适用于流行的在线裁判系统，例如 POJ、HDOJ、SGU 和 ACM-ICPC Live Archive。 但是，我忘记了我用来解决这些问题的算法！ 我将通过竞赛来组织代码和解决方案。 一些...
2013ACM-ICPC杭州赛区全国邀请赛——Random Walk
Confidence的专栏
05-09 1550
n个点，按照题中给的公式可以求出任意两个点转移的概率。求从1到n的期望转移次数
2013ACM-ICPC杭州赛区全国邀请赛——X-Boxes
Confidence的专栏
05-06 1358
给定n个球（标号1 - n）和m个盒子（标号1 - m），现在要求：当一个球x放到了i号盒子中，如果有i + 1盒子，那么x * 2球必须在里边；如果有i - 1盒子，必然有x / 2球在里边。求能1号盒子中球数量最大值 ( 1≤n≤10 10000, 2≤k≤25 )
2013ACM-ICPC杭州赛区全国邀请赛（题目重现)
u010211389的专栏
08-10 2921
Robot Time Limit: 8000/4000 MS (Java/Others)    Memory Limit: 102400/102400 K (Java/Others) Total Submission(s): 0    Accepted Submission(s): 0 Problem Description Michael has a telecontrol
2013ACM-ICPC杭州赛区全国邀请赛 ——Robot
Confidence的专栏
05-06 1541
给定一个长度为n的环形（1 - n），初始在1点。m次操作，每次给一个整数t，表示从当前位置向左向右t的概率相等，求最后停在[l, r]的概率是多少
2013ACM-ICPC杭州赛区全国邀请赛（C待补）
qq_43461168的博客
02-16 243
传送门：http://acm.hdu.edu.cn/search.php?field=problem&key=2013ACM-ICPC%BA%BC%D6%DD%C8%FC%C7%F8%C8%AB%B9%FA%D1%FB%C7%EB%C8%FC&source=1&searchmode=source A. Robot 题意：一个n个格子的圈，编号1-n。一个机器人初始在位置1。然后给出m条指令。每条指令，可以控制机器人移动x个位置。但是方向是任意的，可以顺时针可以逆时针。概率相同。问最后
HDU acm-PPT课件
06-07
HDU（杭州电子科技大学）在ACM竞赛领域有着丰富的经验和优秀的教育资源，其提供的ACM课件是初学者入门和进阶的良好资源。 一、基础篇：编程语言与基础算法 在ACM竞赛中，C++是最常用的语言，因其高效、灵活和丰富...
hdu.zip_ACM_hdu
最新发布
09-20
在ACM(国际大学生程序设计竞赛)中，数据结构与算法是至关重要的组成部分，而"杭电"(杭州电子科技大学)的在线判题系统HDU提供了许多这类问题供参赛者练习。"hdu.zip_ACM_hdu"这个压缩包很可能是包含了一道或若干道...
HDU_2010.rar_hdu 2010_hdu 20_hdu acm20
09-14
"hdu"可能是指杭州电子科技大学（Hangzhou Dianzi University），这所学校经常举办ACM编程竞赛，并有自己的在线判题系统——HDU Online Judge，供参赛者提交代码并测试解决方案。 【压缩包子文件的文件名称列表】中...
2013 ACM-ICPC杭州赛区全国邀请赛——题目重现
Bright-XL
08-10 1394
2013 ACM-ICPC杭州赛区全国邀请赛
hdu 4585 Shaolin 2013 ACM-ICPC杭州赛区全国邀请赛——题目重现
逍遥小章
10-02 1194
题意大意：很多人想进少林寺，少林寺最开始只有一个和尚，每个人有有一个武力值，若这个人想进少林，必须和比他先进去的人比武并且武力值最接近他的比武，如果有相同的则选择武力值比他小的，问当他进去的时候要和哪个和尚比武。   思路：正常情况暴力，n=100000         n*n的操作肯定不行的，最少要优化到nlogn，联想到map，map里的数会自动排序，它是由红黑树实现的，能实
Transformations ACM
jimzhai的专栏
02-24 654
Transformations A square pattern of size N x N (1 <= N <= 10) black and white square tiles is transformed into another square pattern. Write a program that will recognize the minimum transformation
HDU 4578 Transformation（2013 杭州赛区全国邀请赛 1003）
ddyyxx的程序员之路
08-14 2131
http://acm.hdu.edu.cn/showproblem.php?pid=4578 Transformation Time Limit: 15000/8000 MS (Java/Others)    Memory Limit: 65535/65536 K (Java/Others) Total Submission(s): 1483    Accepted Subm
hdu 4576 robot 2013 ACM-ICPC杭州赛区全国邀请赛——题目重现-1001-robot
柯安的专栏
08-10 1675
题意：有一个机器人他的初始位置是1，在一个环形轨道上(轨道的长度n 既然决定暴力，那么就是一步一步的算机器人走到轨道的某处的情况数了，可以一步步递推每次接到指令后就枚举轨道的位置，复杂度是n*m=2*10^8很凶残啊。。。。 输入格式是先输入 n,m,l,r,接下来再有m行输入指令的w值，若n,m,l,r,都是0则输入结束 需要注意的是所有的情况数将会非常庞大可以爆掉int，会WA，但是用l
HDU4577(2013年ACM杭州赛区邀请赛B题)
ACdreamer
08-10 2350
题目：X-Boxes   import java.math.BigInteger; import java.util.*; public class Main { public static void main(String[] args) { Scanner cin=new Scanner(System.in); long[] f = new long[30];
HDU ACM 4578 Transformation->段树-间隔的变化
weixin_34380948的博客
10-06 90
分析：复杂的经营分部树。只有一个查询操作，这是要求[l，r]的数量之间p钍总和。并不是所有的查询所有节点，会议TLE。最好的是查询部件[a。b]。所有这个区间值我们是平等的，即能返回（b-a+1）*val 的值。区间内全部值都同样的情况的区间。对于置初值和加乘操作。分两种情况：1、当为置初值操作。直接覆盖区间就可以。并把标记的加乘操作赋为初始值。2、当为加乘操作时。先推断当前区间段是否为同...
杭电1041 Computer Transformation
hui1570的专栏
07-28 458
Computer Transformation Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Submission(s): 5856    Accepted Submission(s): 2122 Problem Description A sequ
hdu 4578 Transformation（线段树+多种操作）
Darkarts
08-01 563
TransformationProblem DescriptionYuanfang is puzzled with the question below: There are n integers, a1, a2, …, an. The initial values of them are 0. There are four kinds of operations. Operation 1:

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值