CodeForces - 1114F 线段树+欧拉函数+状压

如果你想初步了解这种状压的思想可以看看：牛客 15667 统计颜色 线段树+状压

一开始不想求逆元 想着用乘法维护  及    其中pj是ai积所包含的质因子 k是每个质因子出现的次数，所有我开了63*2个树状数组来维护每个质因子的个数，然后果断超时了 。。 

后来看了题解才知道  就是从欧拉函数的原始公式下手        那我们要维护的就是区间l到r的积以及 区间出现的质因数种数（对应到值） 第一个很好维护 那第二个怎么维护呢   需要状压   我们可以发现 300以内 有62个素数  接近ll的范围  那我可以用一个ll来表示 这个位置有哪些素数  然后用或运算来维护区间有哪些素数 

1. 注意乘法标记要初始化为1 并且 下推以后也要置为一
2. 1<<i大于int范围时应该用  1ll<<i

#include<bits/stdc++.h>
using namespace std;
const int N = 4e5+100;
typedef long long ll;
const ll mod = 1e9+7;
ll sum[N<<2],pi[N<<2],tag[N<<2],ptag[N<<2],inv[320],las[320],pri[320];
int  cnt = -1,vis[320];
#define pa pair<ll,ll>
void init(){
	inv[1]=1;
	for(int i = 2; i <= 300; i++){
		inv[i]=(mod-mod/i)*inv[mod%i]%mod;
		if(!vis[i]){
		pri[++cnt]=i;
		for(int j = i+i; j <= 300; j+=i)
			vis[j]=1;
		}
	} 
	//for(int i = 0; i <= cnt; i++) printf("i=%d pri[i]=%lld\n",i,pri[i]);
	for(int i = 0; i <= cnt; i++) las[i]=(pri[i]-1)*inv[pri[i]]%mod; 
}
ll qpow(ll a,ll b){
	ll ret = 1;
	while(b){
		if(b&1) ret=ret*a%mod;
		a=a*a%mod;
		b>>=1; 
	}return ret%mod;
}
void pushup(int id){
	sum[id]=sum[id<<1]*sum[id<<1|1]%mod;
	pi[id]=pi[id<<1]|pi[id<<1|1];
}
void modify(int id,ll len,ll val,ll f){
	sum[id]=sum[id]*qpow(val,len)%mod;
	pi[id]|=f;
	tag[id]=tag[id]*val%mod;
	ptag[id]|=f;
}
void pushdown(int id,int l,int r){
	int mid = l+r>>1;
	if(tag[id]!=1||ptag[id]){
		modify(id<<1,(ll)(mid-l+1),tag[id],ptag[id]);
		modify(id<<1|1,(ll)(r-mid),tag[id],ptag[id]);
		ptag[id]=0;
		tag[id]=1ll;
	}
}
void build(int id,int l,int r){
	tag[id]=1ll;
	if(l==r){
		scanf("%lld",&sum[id]);
		for(int i = 0; i <= cnt; i++)
		if(sum[id]%pri[i]==0) 
		pi[id]|=(1ll<<i); 
		//printf("bsum[id]=%lld\n",sum[id]);
		return ;
	}
	int mid = l+r>>1;
	build(id<<1,l,mid);build(id<<1|1,mid+1,r);
	pushup(id);
}
void update(int id,int l,int r,int L,int R,ll val,ll f){
	if(L<=l&&R>=r){
		modify(id,(ll)(r-l+1),val,f);
		//printf("tag[id]=%lld usum[id]=%lld\n",tag[id],sum[id]);
		return;
	}
	pushdown(id,l,r);
	int mid = l+r>>1;
	if(L<=mid) update(id<<1,l,mid,L,R,val,f);
	if(R>mid) update(id<<1|1,mid+1,r,L,R,val,f);
	pushup(id);
}
pa query(int id,int l,int r,int L,int R){
	if(L<=l&&R>=r){
		return make_pair(sum[id],pi[id]);		
	}	
	int mid = l+r>>1;
	pa ans=make_pair(1ll,0),q;
	pushdown(id,l,r);
	if(L<=mid){
		q=query(id<<1,l,mid,L,R);
		ans.first*=q.first,ans.second|=q.second;
		ans.first%=mod;
	}
	if(R>mid){
		q=query(id<<1|1,mid+1,r,L,R);
		ans.first*=q.first,ans.second|=q.second;
		ans.first%=mod;
	}
	return ans;
}
int main(){
	init();
	int n,m;
	scanf("%d%d",&n,&m);
	build(1,1,n);
	char op[30];
	int l,r;
	ll x;
	for(int i = 1; i <= m; i++){
		scanf("%s%d%d",op,&l,&r);
		if(op[0]=='T'){
			pa ret = query(1,1,n,l,r);
			ll ans = ret.first;
			for(int i = 0; i <= cnt; i++)
			if((ret.second>>i)&1) 
			ans=ans*las[i]%mod;
			printf("%lld\n",ans);
		}else{
			scanf("%lld",&x);
			ll c=0;
			for(int i = 0; i <= cnt; i++)
				if(x%pri[i]==0) 
				c|=(1ll<<i);
			update(1,1,n,l,r,x,c);
		}
	}
	return 0;
}