题目链接:http://47.96.116.66/problem.php?id=1911
解题思路:
①关键在于求欧拉函数
参考自 liuzibujian 的博客 浅谈欧拉函数 https://blog.csdn.net/liuzibujian/article/details/81086324
几个关键结论:
(1)若n是质数,f(n)= n-1
证明:质数和任何小于它的数都互为质数
(2)若m,n互质,则f(m*n)=f(m)*f(n)
证明:m,n互质,则各自的质因子都不相同。
②线段树区间区间暴力修改,另开一个数组记录是否维护的值都是1,都是1的话就不用往下修改,因为f(1)= 1
而一个很大的数不需要很多次修改就会变为1,所以这么一剪枝就可以过了。
也可以不多开这个数组,只要区间和=r-l+1(区间长度) 就行了。
代码:
#include<cstdio>
#include<algorithm>
#include<cstring>
#define ll long long
#define lson rt<<1,l,m
#define rson rt<<1|1,m+1,r
#define mid int m = l+r>>1
#define tl tree[rt<<1]
#define tr tree[rt<<1|1]
using namespace std;
const int N = 1e7+5;
int f[N];
int prime[N>>2];
bool vis[N];
void getf(int n)
{
f[1] = 1;
int tot=0;
for (int i=2;i<=n;i++){
if (!vis[i]){
f[i] = i-1;
prime[++tot] = i;
}
for (int j=1;j<=tot && i*prime[j]<=n;j++){
vis[i*prime[j]] = true;
if (i%prime[j]) f[i*prime[j]] = f[i]*f[prime[j]];
else {
f[i*prime[j]] = f[i]*prime[j];
break;
}
}
}
}
const int MAXN = 5e5+5;
ll tree[MAXN<<2];
bool lower2[MAXN<<2];
void push_up(int rt)
{
tree[rt] = tl + tr;
lower2[rt] = lower2[rt<<1] && lower2[rt<<1|1];
}
void update(int L,int R,int rt,int l,int r)
{
if (lower2[rt]) return ;
if (l==r){
tree[rt] = f[tree[rt]];
if (tree[rt]<2) lower2[rt] = true;
return ;
}
mid;
if (L<=m) update(L,R,lson);
if (R>m) update(L,R,rson);
push_up(rt);
}
ll query(int L,int R,int rt,int l,int r)
{
if (L<=l && r<=R){
return tree[rt];
}
mid;
ll ans = 0;
if (L<=m) ans += query(L,R,lson);
if (R>m) ans += query(L,R,rson);
return ans;
}
void build(int rt,int l,int r)
{
if (l==r){
scanf("%lld",&tree[rt]);
lower2[rt] = tree[rt]<2? true:false;
return ;
}
mid;
build(lson);
build(rson);
push_up(rt);
}
int main()
{
getf(10000000);
int n;
while(~scanf("%d",&n)){
build(1,1,n);
int m,op,l,r;
scanf("%d",&m);
while (m--){
scanf("%d %d %d",&op,&l,&r);
if (op==1){
update(l,r,1,1,n);
}
else printf("%lld\n",query(l,r,1,1,n));
}
}
return 0;
}