建设部网站防排烟规范网络域名综合查询

Problem - E - Codeforces

思路：观察给定的函数，其实就是求与这个数互质的数的个数，即欧拉函数，我们发现一个数迭代欧拉函数不会很多，那么对于第一个操作来说我们可以直接暴力修改，而对于第二个操作来说，就是求l,r的最近公共祖先，那么我们可以用线段树维护区间的最近公共祖先，并且由于迭代的次数很少，所以并不需要建图，直接暴力跳跃求最近公共祖先即可，那么最总的答案就是用l到r的深度之和-最近公共祖先的深度，乘以区间长度，这就是把这每个点跳跃到最近公共祖先的花费

一定不要建图，建图会MLE

// Problem: E. Li Hua and Array
// Contest: Codeforces - Codeforces Round 864 (Div. 2)
// URL: https://codeforces.com/problemset/problem/1797/E
// Memory Limit: 512 MB
// Time Limit: 3000 ms#include<bits/stdc++.h>
#include<sstream>
#include<cassert>
#define fi first
#define se second
#define i128 __int128
using namespace std;
typedef long long ll;
typedef double db;
typedef pair<int,int> PII;
const double eps=1e-7;
const int N=5e6+7 ,M=6e7+7, INF=0x3f3f3f3f,mod=1e9+7,mod1=998244353;
const long long int llINF=0x3f3f3f3f3f3f3f3f;
inline ll read() {ll x=0,f=1;char c=getchar();while(c<'0'||c>'9') {if(c=='-') f=-1;c=getchar();}
while(c>='0'&&c<='9') {x=(ll)x*10+c-'0';c=getchar();} return x*f;}
inline void write(ll x) {if(x < 0) {putchar('-'); x = -x;}if(x >= 10) write(x / 10);putchar(x % 10 + '0');}
inline void write(ll x,char ch) {write(x);putchar(ch);}
void stin() {freopen("in_put.txt","r",stdin);freopen("my_out_put.txt","w",stdout);}
bool cmp0(int a,int b) {return a>b;}
template<typename T> T gcd(T a,T b) {return b==0?a:gcd(b,a%b);}
template<typename T> T lcm(T a,T b) {return a*b/gcd(a,b);}
void hack() {printf("\n----------------------------------\n");}int T,hackT;
int n,m,k;
int cost[N];
int pr[3000],cnt;
bool st[3000];
int vis[N];
int depth[N];
struct Node{int l,r;int start;int sum;int add;
};
Node tr[N*4];int lca(int a,int b) {while(a!=b) {if(a>b) a=vis[a];else b=vis[b];}return a;
}void pushup(int u) {tr[u].sum=tr[u<<1].sum+tr[u<<1|1].sum;tr[u].add=tr[u<<1].add+tr[u<<1|1].add;tr[u].start=lca(tr[u<<1].start,tr[u<<1|1].start);
}void build(int u,int l,int r) {if(l==r) {if(depth[cost[l]]==0) tr[u]={l,r,cost[l],depth[cost[l]],1};else tr[u]={l,r,cost[l],depth[cost[l]],0};}else {tr[u]={l,r,INF,0,0};int mid=l+r>>1;build(u<<1,l,mid),build(u<<1|1,mid+1,r);pushup(u);}
}void modify(int u,int l,int r) {if(tr[u].r-tr[u].l+1==tr[u].add) return ;if(tr[u].l==tr[u].r) {int tp=vis[tr[u].start];if(depth[tp]==0) tr[u]={tr[u].l,tr[u].r,tp,depth[tp],1};else tr[u]={tr[u].l,tr[u].r,tp,depth[tp],0}; }else {int mid=tr[u].l+tr[u].r>>1;if(l<=mid) modify(u<<1,l,r);if(r>mid) modify(u<<1|1,l,r);pushup(u);}
}int query_sum(int u,int l,int r) {if(tr[u].l>=l&&tr[u].r<=r) return tr[u].sum;else {int mid=tr[u].l+tr[u].r>>1;int res=0;if(l<=mid) res+=query_sum(u<<1,l,r);if(r>mid) res+=query_sum(u<<1|1,l,r);return res;}
}int query_lca(int u,int l,int r) {if(tr[u].l>=l&&tr[u].r<=r) return tr[u].start;else {int mid=tr[u].l+tr[u].r>>1;if(l<=mid&&r>mid) {int a=query_lca(u<<1,l,r);int b=query_lca(u<<1|1,l,r);return lca(a,b);}else if(l<=mid) return query_lca(u<<1,l,r);else return query_lca(u<<1|1,l,r);}
}void get(int x) {int tx=x;int res=x;for(int i=0;i<cnt&&pr[i]<=x/pr[i];i++) {if(x%pr[i]==0) {res=res/pr[i]*(pr[i]-1);while(x%pr[i]==0) x/=pr[i];}}if(x!=1) res=res/x*(x-1);vis[tx]=res;depth[tx]=depth[res]+1;
}void init() {int t=sqrt(5000000);for(int i=2;i<=t;i++) {if(!st[i]) pr[cnt++]=i;for(int j=0;pr[j]<=t/i;j++) {st[pr[j]*i]=true;if(i%pr[j]==0) break;}}depth[1]=0;for(int i=2;i<=5000000;i++) get(i);
}void solve() {n=read(),m=read();for(int i=1;i<=n;i++) cost[i]=read();init();build(1,1,n);while(m--) {int op=read();int l=read(),r=read();if(op==1) {modify(1,l,r);}else if(op==2) {int s=query_lca(1,l,r);int sum=query_sum(1,l,r);printf("%d\n",sum-depth[s]*(r-l+1));}}
}   int main() {// init();// stin();// ios::sync_with_stdio(false); // scanf("%d",&T);T=1; while(T--) hackT++,solve();return 0;       
}