Codeforces Round 37 F. SUM and REPLACE【线段树】

Let D(x) be the number of positive divisors of a positive integer x. For example, D(2) = 2 (2 is divisible by 1 and 2), D(6) = 4 (6 is divisible by 123 and 6).

You are given an array a of n integers. You have to process two types of queries:

  1. REPLACE l r — for every  replace ai with D(ai);
  2. SUM l r — calculate .

Print the answer for each SUM query.

Input

The first line contains two integers n and m (1 ≤ n, m ≤ 3·105) — the number of elements in the array and the number of queries to process, respectively.

The second line contains n integers a1a2, ..., an (1 ≤ ai ≤ 106) — the elements of the array.

Then m lines follow, each containing 3 integers tiliri denoting i-th query. If ti = 1, then i-th query is REPLACE li ri, otherwise it's SUM liri (1 ≤ ti ≤ 21 ≤ li ≤ ri ≤ n).

There is at least one SUM query.

Output

For each SUM query print the answer to it.

Example
input
7 6
6 4 1 10 3 2 4
2 1 7
2 4 5
1 3 5
2 4 4
1 5 7
2 1 7
output
30
13
4
22

思路:

区间修改线段树。但是这一题不能使用懒惰标记,因为需要精确修改区间内所有点,然后向上更新。但是为了降低时间消耗,减少不必要的更新操作,我们引入了一个ok数组。
当d[i]等于1或2时就没有必要进行d操作了,ok[i]=1.
当区间内的所有点都是1或2时,那么这个区间就不用更新了。ok数组向上更新即可。
d数组可以用筛法求。第一层循环枚举因子,第二层循环枚举因子的倍数。

代码:

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<iostream>
#include<queue>
using namespace std;
int n,m;
int d[1000005],ok[4*300005];
long long sum[4*300005];
void init()
{
    memset(sum,0,sizeof sum);
    memset(ok,0,sizeof ok);
    for(int i=1;i<=1000000;i++)
    {
        for(int j=i;j<=1000000;j+=i)
            d[j]++;
    }
}
void build(int v,int l,int r)
{
    if(l==r)
    {
        long long x;
        scanf("%lld",&sum[v]);
        if(sum[v]==1 || sum[v]==2) ok[v]=1;
        return ;
    }
    int mid=(l+r)/2;
    build(v*2,l,mid);
    build(v*2+1,mid+1,r);
    sum[v]=sum[v*2]+sum[v*2+1];
    ok[v]=ok[v*2] & ok[v*2+1];
}

void updata(int v,int l,int r,int ql,int qr)
{
    if(ok[v]) return ;
    if(l==r)
    {
        sum[v]=d[sum[v]];
        if(sum[v]==1 || sum[v]==2) ok[v]=1;
        return ;
    }
    int mid=(l+r)/2;
    if(ql<=mid) updata(v*2,l,mid,ql,qr);
    if(qr>mid) updata(v*2+1,mid+1,r,ql,qr);
    sum[v]=sum[v*2]+sum[v*2+1];
    ok[v]=ok[v*2] & ok[v*2+1];
}
long long query(int v,int l,int r,int ql,int qr)
{
    if(ql<=l && r<=qr)
        return sum[v];
    int mid=(l+r)/2;
    long long s=0;
    if(ql<=mid) s+=query(v*2,l,mid,ql,qr);
    if(qr>mid) s+=query(v*2+1,mid+1,r,ql,qr);
    return s;
}
int main()
{
    scanf("%d%d",&n,&m);
    init();
    build(1,1,n);
    int op,l,r;
    for(int i=0;i<m;i++)
    {
        scanf("%d%d%d",&op,&l,&r);
        if(op==1) updata(1,1,n,l,r);
        else printf("%lld\n",query(1,1,n,l,r));
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值