Master of GCD(分块大法)--优雅的暴力

题目： Master of GCD

题意：

 数据：

2
5 3
1 3 2
3 5 2
1 5 3
6 3
1 2 2
5 6 2
1 6 2

题意：角标1~n之间的数默认为1，通过m次修改区间[l, r]乘以2或者3，最后输出坐标1~n之间最大公约数。

思路：更新乘以2或者3的次数，分别找到最小更新2的次数p和更新3的次数q，然后2^p*3^q就是最后的结果。

代码：

#include<iostream>
#include<cstdio>
#include<cctype>
#include<cmath>
using namespace std;
typedef long long ll;
const int MAXN=100005;
const ll mod = 998244353;
const ll INF = 1e18;
ll r[MAXN],l[MAXN],inc1[MAXN];
ll a1[MAXN],a2[MAXN], inc2[MAXN],belong[MAXN],block;
int n,m,num;
void build()//初始化
{
    block=sqrt(n);
    num=n/block;
    if(n%block) num++;
    for(int i=1; i<=num; i++)
    {
        l[i]=(i-1)*block+1;
        r[i]=i*block;
    }
    for(int i=1; i<=n; i++)
        belong[i]=(i-1)/block+1;
    r[num]=n;
}
void updata2(int x,int y)//找2
{
    if(belong[x]==belong[y])
    {
        for(int i=x; i<=y; i++) a1[i]++;
        return ;
    }
    for(int i=x; i<=r[belong[x]]; i++)
        a1[i]++;
    for(int i=belong[x]+1; i<=belong[y]-1; i++)
        inc1[i]++;
    for(int i=l[belong[y]]; i<=y; i++)
        a1[i]++;
}
 
ll Quick_Power(ll a,ll b)//快速幂
{
    ll res = 1;
    a %= mod;
    while(b){
        if(b&1){
            res = (res * a) % mod;
        }
        a = (a * a) % mod;
        b >>= 1;
    }
    return res%mod;
}
 
ll query1(int x)//查找2的次数
{
    return (a1[x]+inc1[belong[x]]);
}
void updata3(int x,int y)//更新3的次数
{
    if(belong[x]==belong[y])
    {
        for(int i=x; i<=y; i++) a2[i]++;
        return ;
    }
    for(int i=x; i<=r[belong[x]]; i++)
        a2[i]++;
    for(int i=belong[x]+1; i<=belong[y]-1; i++)
        inc2[i]++;
    for(int i=l[belong[y]]; i<=y; i++)
        a2[i]++;
}
 
ll query2(int x)//查找3的次数
{
    return (a2[x]+inc2[belong[x]]);
}
int main()
{
    int T;
    scanf("%d", &T);
    while(T--)
    {
        scanf("%d %d", &n, &m);
        for(int i=1; i<=n; i++)
        {
            a1[i]=0;
            inc1[i] = 0;
            a2[i] = 0;
            inc2[i] = 0;
        }
        build();
        for(int i=1; i<=m; i++)
        {
            int l, r, w;
            scanf("%d %d %d", &l, &r, &w);
            if(w == 2)
            updata2(l, r);
            else updata3(l, r);
        }
        ll ma1 = INF, t1, t2, ma2 = INF;
        for(int i = 1; i <= n; i++)
        {
            t1 = query1(i);
            t2 = query2(i);
           // cout<<t<<endl;
            if(ma2 > t2){
                ma2 = t2;
            }
            if(ma1 > t1)
            {
                ma1 = t1;
            }
        }
        ll k = Quick_Power(2, ma1);
        ll f = Quick_Power(3, ma2);
        printf("%lld\n", (k*f)%mod);
    }
    return 0;
}