Rikka with Phi 【HDU - 5634】【线段树+欧拉函数】

题目链接

---

  很好的一道题，也算是开阔了我的思维，一开始，我的想法是既然是区间求欧拉，那么到一定地步的时候，数的欧拉值就会降到1，那我们维护区间值等于区间长度不就是可以了吗？然后T了。

  然后我再想，是不是哪里可以优化？毕竟还有另外一个条件没用优化（区间赋予相同值），那我们需要怎么做？区间值为“1”就是意味着区间上每个值都相等，那么我们让区间上如果每个值都相等，那么直接赋予其等价欧拉，然后利用延迟标记lazy[]向下传递值不是更好吗？那么问题来了，怎么维护区间以及其子区间的值是否相等，我们从叶子节点出发，向上递推，lazy[]标记也随着向上更新而更新，如果节点的左右子树是相等的lazy[]那么，岂不是可以赋给目前节点它们的lazy[]不是吗？那么，对应左右子树不想等的情况，我们赋予其lazy[]值=0就是不相等，不用传递了的lazy[]。

  这样的做法有没有问题呢？我们利用这个做法模拟一下我们的题目，遇到区间更新欧拉函数（区间每个值变成自己的欧拉函数），那么，我们向下找区间，遇到刚好包含在查询区间内的节点区间，我们就判断它的lazy[]，如果不为0，说明它的子节点都是相等的，那么，我们就可以直接更新lazy[]标记为本身的欧拉函数，不然，如果我们查到叶子节点，他必定有一个确定的欧拉函数。

  接下来的下一个查询是区间所有值都变成X，那么就好解了，我们向下pushdown()传递lazy[]的同时，我们遇到包含的区间就是直接赋值，更新lazy[]以及线段树的区间和不就是可以了吗？

  最后说的是查询操作，查询倒也不难，我们看到要查询的是个区间，就如同平时线段树区间查询一样即可，所以也就不难了。

  还有一件事，欧拉函数记得预处理哦，不然时间复杂度还是很高的。

  细节！哦，对了，最后有一件事，我当时WA了半个多小时，有的同学（其实就是我了……）喜欢用long long * int，但是在这道题里会爆精度，尤其是(r-l+1)以及遇上(mid-l+1)或者(r-mid)求长度的时候，会爆精度，大家注意了！！！

---

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#define lowbit(x) ( x&(-x) )
#define INF 0x3f3f3f3f
#define pi 3.141592653589793
#define e 2.718281828459045
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
const int MAXX = 1e7+1;
const int maxN = 3e5+5;
int N, M, euler[MAXX];
ll tree[maxN<<2], lazy[maxN<<2];
void Euler()
{
    for(int i=1; i<MAXX; i++) euler[i] = i;
    for(int i=2; i<MAXX; i++)
    {
        if(euler[i] == i)
        {
            for(int j=i; j<MAXX; j+=i)
            {
                euler[j] = euler[j]/i*(i-1);
            }
        }
    }
}
void pushup(int rt)
{
    tree[rt] = tree[rt<<1] + tree[rt<<1|1];
    if(lazy[rt<<1] == lazy[rt<<1|1]) lazy[rt] = lazy[rt<<1];
    else lazy[rt] = 0;
}
void buildTree(int rt, int l, int r)
{
    if(l == r)
    {
        scanf("%lld", &tree[rt]);
        lazy[rt] = tree[rt];
        return;
    }
    int mid = (l + r)>>1;
    buildTree(rt<<1, l, mid);
    buildTree(rt<<1|1, mid+1, r);
    pushup(rt);
}
void pushdown(int rt, int l, int r)
{
    if(lazy[rt])
    {
        lazy[rt<<1] = lazy[rt<<1|1] = lazy[rt];
        int mid = (l + r)>>1;
        tree[rt<<1] = (ll)(mid - l + 1)*lazy[rt];
        tree[rt<<1|1] = (ll)(r - mid)*lazy[rt];
    }
}
void update_Euller(int rt, int l, int r, int ql, int qr)
{
    if(lazy[rt] && ql<=l && qr>=r)
    {
        tree[rt] = euler[lazy[rt]]*(ll)(r - l + 1);
        lazy[rt] = euler[lazy[rt]];
        return;
    }
    pushdown(rt, l, r);
    int mid = (l + r)>>1;
    if(ql>mid) update_Euller(rt<<1|1, mid+1, r, ql, qr);
    else if(qr<=mid) update_Euller(rt<<1, l, mid, ql, qr);
    else
    {
        update_Euller(rt<<1, l, mid, ql, qr);
        update_Euller(rt<<1|1, mid+1, r, ql, qr);
    }
    pushup(rt);
}
void update_val(int rt, int l, int r, int ql, int qr, ll val)
{
    if(ql<=l && qr>=r)
    {
        lazy[rt] = val;
        tree[rt] = val*(ll)(r - l + 1);
        return;
    }
    pushdown(rt, l, r);
    int mid = (l + r)>>1;
    if(ql>mid) update_val(rt<<1|1, mid+1, r, ql, qr, val);
    else if(qr<=mid) update_val(rt<<1, l, mid, ql, qr, val);
    else
    {
        update_val(rt<<1, l, mid, ql, qr, val);
        update_val(rt<<1|1, mid+1, r, ql, qr, val);
    }
    pushup(rt);
}
ll query(int rt, int l, int r, int ql, int qr)
{
    if(ql<=l && qr>=r) return tree[rt];
    pushdown(rt, l, r);
    int mid = (l + r)>>1;
    if(ql>mid) return query(rt<<1|1, mid+1, r, ql, qr);
    else if(qr<=mid) return query(rt<<1, l, mid, ql, qr);
    else return query(rt<<1, l, mid, ql, mid) + query(rt<<1|1, mid+1, r, mid+1, qr);
}
int main()
{
    Euler();
    int T;  scanf("%d", &T);
    while(T--)
    {
        scanf("%d%d", &N, &M);
        buildTree(1, 1, N);
        int op, L, R;
        ll X;
        while(M--)
        {
            scanf("%d", &op);
            if(op == 1)
            {
                scanf("%d%d", &L, &R);
                update_Euller(1, 1, N, L, R);
            }
            else if(op == 2)
            {
                scanf("%d%d%lld", &L, &R, &X);
                update_val(1, 1, N, L, R, X);
            }
            else if(op == 3)
            {
                scanf("%d%d", &L, &R);
                printf("%lld\n", query(1, 1, N, L, R));
            }
        }
    }
    return 0;
}