2019ICPC银川 Pot （线段树）

$\quad$这道题就是$1-n$个数$m$个操作，其中一个就是让区间的数都乘上$2-10$中的一个数，然后求区间的任意一个满足可以被$p^{m}$整除，不能被$p^{m+1}$整除，$p$是一个质数，求最大的$m$。

$\quad$因为只有乘$2-10$，$2-10$直接的质数只有$2,3,5,7$，要保证m最大我们肯定只考虑$2,3,5,7$，这样看每次乘的数然后分解为$2,3,5,7$，建立四颗线段树维护最值就行了。

AC代码：

#include <cstdio>
#include <vector>
#include <queue>
#include <cstring>
#include <cmath>
#include <map>
#include <set>
#include <string>
#include <iostream>
#include <algorithm>
#include <iomanip>
using namespace std;
#define sd(n) scanf("%d",&n)
#define sdd(n,m) scanf("%d%d",&n,&m)
#define sddd(n,m,k) scanf("%d%d%d",&n,&m,&k)
#define pd(n) printf("%d\n", n)
#define pc(n) printf("%c", n)
#define pdd(n,m) printf("%d %d", n, m)
#define pld(n) printf("%lld\n", n)
#define pldd(n,m) printf("%lld %lld\n", n, m)
#define sld(n) scanf("%lld",&n)
#define sldd(n,m) scanf("%lld%lld",&n,&m)
#define slddd(n,m,k) scanf("%lld%lld%lld",&n,&m,&k)
#define sf(n) scanf("%lf",&n)
#define sc(n) scanf("%c",&n)
#define sff(n,m) scanf("%lf%lf",&n,&m)
#define sfff(n,m,k) scanf("%lf%lf%lf",&n,&m,&k)
#define ss(str) scanf("%s",str)
#define rep(i,a,n) for(int i=a;i<=n;i++)
#define per(i,a,n) for(int i=n;i>=a;i--)
#define mem(a,n) memset(a, n, sizeof(a))
#define debug(x) cout << #x << ": " << x << endl
#define pb push_back
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define mod(x) ((x)%MOD)
#define gcd(a,b) __gcd(a,b)
#define lowbit(x) (x&-x)
typedef pair<int,int> PII;
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
const int MOD = 1e9 + 7;
const double eps = 1e-9;
const ll INF = 0x3f3f3f3f3f3f3f3fll;
const int inf = 0x3f3f3f3f;
inline int read()
{
    int ret = 0, sgn = 1;
    char ch = getchar();
    while(ch < '0' || ch > '9')
    {
        if(ch == '-')
            sgn = -1;
        ch = getchar();
    }
    while (ch >= '0' && ch <= '9')
    {
        ret = ret*10 + ch - '0';
        ch = getchar();
    }
    return ret*sgn;
}
inline void Out(int a)    //Êä³öÍâ¹Ò
{
    if(a>9)
        Out(a/10);
    putchar(a%10+'0');
}

int qpow(int m, int k, int mod)
{
    int res = 1, t = m;
    while (k)
    {
        if (k&1)
            res = res * t % mod;
        t = t * t % mod;
        k >>= 1;
    }
    return res;
}

ll gcd(ll a,ll b)
{
    return b==0?a : gcd(b,a%b);
}

ll lcm(ll a,ll b)
{
    return a*b/gcd(a,b);
}
#define rtl rt<<1
#define rtr rt<<1|1
using namespace std;
const int MAX_N=2e5+5;
struct node
{
    int l,r,max1,lazy;
} tr2[MAX_N<<2],tr3[MAX_N<<2],tr5[MAX_N<<2],tr7[MAX_N<<2];
inline void update(node *a,int rt)
{
    a[rt].max1=max(a[rtl].max1,a[rtr].max1);
}
void build(node *a,int rt,int l,int r)
{
    a[rt].l=l,a[rt].r=r;
    a[rt].lazy=a[rt].max1=0;
    if(l==r)
        return;
    int mid=(l+r)>>1;
    build(a,rtl,l,mid),build(a,rtr,mid+1,r);
    update(a,rt);
}

void pushdown(node *a,int rt)
{
    if(a[rt].l==a[rt].r)
    {
        a[rt].lazy=0;
        return ;
    }
    a[rtl].max1+=a[rt].lazy;
    a[rtr].max1+=a[rt].lazy;
    a[rtl].lazy+=a[rt].lazy;
    a[rtr].lazy+=a[rt].lazy;
    a[rt].lazy=0;
}

void add(node *a,int rt,int l,int r,int x)
{
    if(a[rt].l==l&&a[rt].r==r)
    {
        a[rt].max1+=x;
        a[rt].lazy+=x;
        return;
    }
    if(a[rt].lazy)
        pushdown(a,rt);
    int mid=(a[rt].l+a[rt].r)>>1;
    if(r<=mid)
        add(a,rtl,l,r,x);
    else if(l>mid)
        add(a,rtr,l,r,x);
    else
    {
        add(a,rtl,l,mid,x);
        add(a,rtr,mid+1,r,x);
    }
    update(a,rt);
}

int query(node *a,int rt,int l,int r)
{
    if(a[rt].lazy)
        pushdown(a,rt);
    if(a[rt].l==l&&a[rt].r==r)
        return a[rt].max1;
    int mid=(a[rt].l+a[rt].r)>>1;
    if(r<=mid)
        return query(a,rtl,l,r);
    if(l>mid)
        return query(a,rtr,l,r);
    return max(query(a,rtl,l,mid),query(a,rtr,mid+1,r));
}
int main()
{
    int n,q;
    sd(n);
    build(tr2,1,1,n);
    build(tr3,1,1,n);
    build(tr5,1,1,n);
    build(tr7,1,1,n);
    sd(q);
    while(q--)
    {
        char s[10];
        ss(s);
        if(s[1]=='U')
        {
            int l,r,x;
            sddd(l,r,x);
            if(x==2)
                add(tr2,1,l,r,1);
            else if(x==3)
                add(tr3,1,l,r,1);
            else if(x==4)
                add(tr2,1,l,r,2);
            else if(x==5)
                add(tr5,1,l,r,1);
            else if(x==6)
                add(tr2,1,l,r,1),add(tr3,1,l,r,1);
            else if(x==7)
                add(tr7,1,l,r,1);
            else if(x==8)
                add(tr2,1,l,r,3);
            else if(x==9)
                add(tr3,1,l,r,2);
            else if(x==10)
                add(tr2,1,l,r,1),add(tr5,1,l,r,1);
        }
        else
        {
            int l,r;
            sdd(l,r);
            int max1=max(query(tr2,1,l,r),query(tr3,1,l,r));
            int max2=max(query(tr5,1,l,r),query(tr7,1,l,r));
            printf("ANSWER %d\n",max(max1,max2));
        }
    }
    return 0;
}