[Fjwc2018]欧拉函数——欧拉筛、线段树、bitset

最新推荐文章于 2022-09-13 21:57:59 发布

偶耶XJX

最新推荐文章于 2022-09-13 21:57:59 发布

阅读量204

点赞数

分类专栏：信息竞赛解题文章标签：算法数据结构线段树数学 bitset

本文链接：https://blog.csdn.net/weixin_43960287/article/details/122482703

版权

信息竞赛解题专栏收录该内容

195 篇文章 5 订阅

订阅专栏

[Fjwc2018]欧拉函数

题解

首先对于第一种询问，易知和的范围在 $2 e 9$ 以内，我们直接枚举质因子求 $\phi$ 即可。

这么做看起来是 $O(\sqrt{S})$ 的，实际上由于数是随机的，期望下时间是 $O(\frac{\sqrt{S}}{\ln S})$ ，再加上预处理质数就可以达到 $O(\frac{\sqrt{S}}{\ln^2 S})\approx O(100)$ 。

然后对于第二种询问，显然乘积是极大的，我们需要知道乘积的所有质因子才能求欧拉函数。

注意到 $a_i\le 40000$ ，也就是说最大质因子不超过 40000，而 40000 以内的质数只有 4203 个，我们可以用 $\rm bitset$ 存下每个数的质因子集合，然后在线段树上并起来，这样可以在 $O(\frac{4203}{w}\log n)=O(66\log n)$ 时间以内求出答案的质因子集合。

然后我们需要提取出集合内所有质因子对应的系数的乘积，一个一个提取肯定是不行的。我们可以预处理出每连续13位质因子的所有 $2^{13}$ 种集合的系数乘积，预处理显然只需要 $O(66\times 5\times 2^{13})$ ，然后一次询问乘积就只需要 $O(66\times 5)=O(330)$ 而不是 $O (4203)$ 的时间。

显然算法的瓶颈在于线段树维护+ $\rm bitset$ 那里，所以总时间复杂度是 $O(66q\log n)$ 的，建议用用常数较小的线段树（如zkw）。

代码

#include<bits/stdc++.h>//JZM yyds!!
#define ll long long
#define uns unsigned
#define IF (it->first)
#define IS (it->second)
#define END putchar('\n')
using namespace std;
const int MAXN=100005;
const ll INF=1e18;
inline ll read(){
	ll x=0;bool f=1;char s=getchar();
	while((s<'0'||s>'9')&&s>0){if(s=='-')f^=1;s=getchar();}
	while(s>='0'&&s<='9')x=(x<<1)+(x<<3)+(s^48),s=getchar();
	return f?x:-x;
}
int ptf[30],lpt;
inline void print(ll x,char c='\n'){
	if(x<0)putchar('-'),x=-x;
	ptf[lpt=1]=x%10;
	while(x>9)x/=10,ptf[++lpt]=x%10;
	while(lpt)putchar(ptf[lpt--]^48);
	if(c>0)putchar(c);
}
inline ll lowbit(ll x){return x&-x;}

const ll MOD=1e9+7;
int pr[1000007],np;
bool nop[10000007];//预处理的越多越好
ll phi[10000007],cg[MAXN];
inline ll PHI(ll x){
	ll s=x,res=x;
	for(int i=1;1ll*pr[i]*pr[i]<=s;i++)
		if(s%pr[i]==0){
			res-=res/pr[i];
			while(s%pr[i]==0)s/=pr[i];
			if(s<=10000000)return res/s*phi[s]%MOD;
		}
	if(s>1)res-=res/s;
	return res%MOD;
}
inline ll ksm(ll a,ll b,ll mo){
	ll res=1;
	for(;b;b>>=1,a=a*a%mo)if(b&1)res=res*a%mo;
	return res;
}
const int w=64,M=66;
ll F[M][5][1<<13];
int lg[1<<13];
struct BS{
	uns ll a[M];BS(){}
	BS(uns ll x){memset(a,x,sizeof(a));}
	void OR(int x){a[x/w]|=(1ull<<x%w);}
	BS operator|(const BS&b)const{
		BS res;
		for(int i=0;i<M;i++)res.a[i]=a[i]|b.a[i];
		return res;
	}
	ll cgs(){
		ll res=1,c=(1<<13)-1;
		for(int i=0;i<M;i++){
			uns ll x=a[i];
			for(int j=0;j<5;j++,x>>=13)
				res=res*F[i][j][x&c]%MOD;
		}return res;
	}
}bc[40005];

int n,Q;
ll a[MAXN];
struct itn{
	ll a,b;BS c;itn(){}
	itn(ll A,ll B,BS C){a=A,b=B,c=C;}
}f[MAXN<<1];
int p;
inline void update(int x){
	f[x].a=f[x<<1].a+f[x<<1|1].a;
	f[x].b=f[x<<1].b*f[x<<1|1].b%MOD;
	f[x].c=f[x<<1].c|f[x<<1|1].c;
}
inline void change(int x,ll a){
	for(f[p+x]=itn(a,a,bc[a]),x=(p+x)>>1;x;x>>=1)
		update(x);
}
inline void build(){
	for(int i=0;i<=n+2;i++)f[p+i]=itn(0,1,BS(0));
	for(int i=1;i<=n;i++)f[p+i]=itn(a[i],a[i],bc[a[i]]);
	for(int i=(p+n+1)>>1;i;i--)update(i);
}
inline ll sum(int l,int r){
	ll res=0;
	for(l=p+l-1,r=p+r+1;l^1^r;l>>=1,r>>=1){
		if(~l&1)res+=f[l^1].a;
		if(r&1)res+=f[r^1].a;
	}return res;
}
inline ll phimul(int l,int r){
	ll res=1;BS g=BS(0);
	for(l=p+l-1,r=p+r+1;l^1^r;l>>=1,r>>=1){
		if(~l&1)res=res*f[l^1].b%MOD,g=g|f[l^1].c;
		if(r&1)res=res*f[r^1].b%MOD,g=g|f[r^1].c;
	}
	return res*g.cgs()%MOD;
}
signed main()
{
	nop[0]=nop[1]=1,phi[0]=phi[1]=1;
	for(int a=2;a<=10000000;a++){
		if(!nop[a])pr[++np]=a,phi[a]=a-1;
		for(int i=1;i<=np&&pr[i]*a<=10000000;i++){
			nop[pr[i]*a]=1;
			if(a%pr[i]==0){
				phi[a*pr[i]]=phi[a]*pr[i];
				break;
			}else phi[a*pr[i]]=phi[a]*phi[pr[i]];
		}
	}
	for(int i=1;i<=np&&pr[i]<=40000;i++)
		for(int j=pr[i];j<=40000;j+=pr[i])bc[j].OR(i-1);
	for(int i=1;i<=50000;i++)
		cg[i]=ksm(i,MOD-2,MOD)*phi[i]%MOD;
	for(int i=0;i<13;i++)lg[1<<i]=i;
	for(int i=0;i<M;i++)
		for(int j=0;j<5;j++){
			int x=i*w+j*13+1;
			F[i][j][0]=1;
			for(int s=1;s<(1<<13);s++){
				int t=lg[lowbit(s)];
				F[i][j][s]=F[i][j][s^(1<<t)]*cg[pr[x+t]]%MOD;
			}
		}
	//预处理↑
	n=read(),Q=read();
	for(p=1;p<n+2;p<<=1);
	for(int i=1;i<=n;i++)a[i]=read();
	build();
	while(Q--){
		int op=read(),l=read(),r=read();
		if(op==0)a[l]=r,change(l,r);
		else if(op==1){
			ll s=sum(l,r);
			print(PHI(s));
		}else print(phimul(l,r));
	}
	return 0;
}