#
# L Euler Function
题目大意:
给定 n 个数字,每一个数字不超过 100,m 次询问
- 把一个区间的所有数字乘以 w (w<=100)
- 求一个区间所有数的欧拉函数和 (mod 998244353)
首先明白一个性质:
而 n 也可以表示为
所以
那么 的值怎么求?
把 w 质因子分解,可以发现如果 n 和 w 有同一个质因子 c,那么,如果存在一个质因子 w 有,n 没有,那么,发现 w 的质因子在要乘的区间中都有,那么这道题就变成了简单的区间乘法,区间询问,但是区间中不一定包含 w 的所有质因子,所以就要去找到那些不包含 w 的某个质因子的位置,把这些位置单独拿出来乘以 (c-1),如何找到这些位置呢?
可以考虑开一个 vis 数组,vis [x] 标记一个区间是否都存在质因子 x,那么每次做区间乘法时,就可以把 w 质因子分解,对于每一个质因子 c,都去线段树中找,如果一个区间被 vis 标记了,那么这个区间都存在这个质因子 c,就可以直接进行区间修改,否则如果这个区间没有被标记,就向左右子树都找,直到找到这个位置
这里的 vis 合并时需要遍历 25 个质因子,进行与的操作,总时复:(O (mlogn*25)),时间快要超时,所以可以把 vis 数组改成 bitset,除以一个 32 的常数,稳过
注意区间乘法lazy数组初始化为1!!!
#include <bits/stdc++.h> | |
#define int long long | |
#define endl '\n' | |
#define ls u<<1 | |
#define rs u<<1|1 | |
#define ios ios::sync_with_stdio(0);cin.tie(0);cout.tie(0) | |
using namespace std; | |
typedef long long ll; | |
const ll mod=998244353; | |
const int N=2e5+100; | |
typedef long long ll; | |
int n,m; | |
int arr[N]; | |
struct node{ | |
int l,r; | |
ll val; | |
bitset<100> bit; | |
}tr[N<<2]; | |
ll lazy[N<<2]; | |
int getphi(int x){ | |
int res=x; | |
for(int i=2;i<=sqrt(x);i++){ | |
if(x%i==0){ | |
res=res*(i-1)/i; | |
while(x%i==0) x/=i; | |
} | |
} | |
if(x>1) res=res*(x-1)/x; | |
return res; | |
} | |
int pr[100],tot; | |
bitset<100> st[N]; | |
void init(){ | |
for(int i=2;i<=100;i++){ | |
int flag=0; | |
for(int j=2;j<=sqrt(i);j++){ | |
if(i%j==0){ | |
flag=1; | |
break; | |
} | |
} | |
if(!flag) pr[++tot]=i; | |
} | |
for(int i=1;i<=n;i++){ | |
for(int j=1;j<=tot;j++){ | |
if(arr[i]%pr[j]==0){ | |
st[i][pr[j]]=1; | |
} | |
} | |
} | |
} | |
void pushup(int u){ | |
tr[u].bit=tr[ls].bit&tr[rs].bit; // 维护标记 | |
tr[u].val=(tr[ls].val+tr[rs].val)%mod; // 维护区间欧拉函数和 | |
} | |
void pushdown(int u){ | |
if(lazy[u]!=1){ | |
lazy[ls]=lazy[u]*lazy[ls]%mod; | |
lazy[rs]=lazy[u]*lazy[rs]%mod; | |
tr[ls].val=tr[ls].val*lazy[u]%mod; | |
tr[rs].val=tr[rs].val*lazy[u]%mod; | |
lazy[u]=1; | |
} | |
} | |
void build(int u,int l,int r){ | |
lazy[u]=1; | |
if(l==r){ | |
tr[u]={l,r,getphi(arr[l]),st[l]}; | |
} | |
else{ | |
int mid=l+r>>1; | |
tr[u]={l,r}; | |
build(ls,l,mid); | |
build(rs,mid+1,r); | |
pushup(u); | |
} | |
} | |
void update(int u,int l,int r,int c){ | |
if(tr[u].l>=l && tr[u].r<=r && tr[u].bit[c]){ | |
lazy[u]=lazy[u]*c%mod; | |
tr[u].val=tr[u].val*c%mod; | |
} | |
else if(tr[u].l==tr[u].r){ | |
tr[u].val=tr[u].val*(c-1)%mod; | |
tr[u].bit[c]=1; | |
} | |
else{ | |
pushdown(u); | |
int mid=tr[u].l+tr[u].r>>1; | |
if(l<=mid) update(ls,l,r,c); | |
if(r>mid) update(rs,l,r,c); | |
pushup(u); | |
} | |
} | |
ll query(int u, int l, int r) { | |
if(l<=tr[u].l && tr[u].r<=r) return tr[u].val; | |
pushdown(u); | |
int mid=(tr[u].l+tr[u].r)>>1; | |
ll res=0; | |
if(l<=mid) res=(res+query(ls, l, r))%mod; | |
if(r>mid) res=(res+query(rs, l, r))%mod; | |
return res; | |
} | |
signed main() | |
{ | |
ios; | |
cin>>n>>m; | |
for(int i=1;i<=n;i++) cin>>arr[i]; | |
init(); | |
build(1,1,n); | |
while(m--){ | |
int op,l,r,w; | |
cin>>op; | |
if(op==0){ | |
cin>>l>>r>>w; | |
for(int i=1;i<=25;i++){ | |
if(w%pr[i]==0){ | |
while(w%pr[i]==0){ | |
w/=pr[i]; | |
update(1,l,r,pr[i]); | |
} | |
} | |
} | |
} | |
else{ | |
cin>>l>>r; | |
cout<<query(1,l,r)<<endl; | |
} | |
} | |
return 0; | |
} | |
/* | |
5 5 | |
5 1 6 2 13 | |
0 5 5 25 | |
0 5 5 18 | |
1 3 5 | |
0 1 3 24 | |
1 3 4 | |
*/ |
第二种写法是用并查集维护一个区间中是否所有数字都可以被质因子 c 整除, 表示第 i 个质因子第 j 个位置的数字往右最多可以延申到哪一个位置,假如可以延伸到 k,那么 [当前位置,k-1] 的区间的答案可以直接乘以,然后把第 k 个位置单独乘以,如此往复,去进行区间修改即可
#include <bits/stdc++.h> | |
#define int long long | |
#define endl '\n' | |
#define ls u<<1 | |
#define rs u<<1|1 | |
#define ios ios::sync_with_stdio(0);cin.tie(0);cout.tie(0) | |
using namespace std; | |
typedef long long ll; | |
const ll mod=998244353; | |
const int N=2e5+100; | |
typedef long long ll; | |
int n,m; | |
int arr[N],fa[26][N]; | |
struct node { | |
int l,r; | |
ll val; | |
} tr[N<<2]; | |
ll lazy[N<<2]; | |
int getphi(int x) { | |
int res=x; | |
for(int i=2; i<=sqrt(x); i++) { | |
if(x%i==0) { | |
res=res*(i-1)/i; | |
while(x%i==0) x/=i; | |
} | |
} | |
if(x>1) res=res*(x-1)/x; | |
return res; | |
} | |
void pushup(int u) { | |
tr[u].val=(tr[ls].val+tr[rs].val)%mod; | |
} | |
void pushdown(int u) { | |
if(lazy[u]!=1) { | |
lazy[ls]=lazy[u]*lazy[ls]%mod; | |
lazy[rs]=lazy[u]*lazy[rs]%mod; | |
tr[ls].val=tr[ls].val*lazy[u]%mod; | |
tr[rs].val=tr[rs].val*lazy[u]%mod; | |
lazy[u]=1; | |
} | |
} | |
void build(int u,int l,int r) { | |
lazy[u]=1; | |
if(l==r) tr[u]= {l,r,getphi(arr[l])}; | |
else { | |
int mid=l+r>>1; | |
tr[u]= {l,r}; | |
build(ls,l,mid); | |
build(rs,mid+1,r); | |
pushup(u); | |
} | |
} | |
ll query(int u, int l, int r) { | |
if(l<=tr[u].l && tr[u].r<=r) return tr[u].val; | |
pushdown(u); | |
int mid=(tr[u].l+tr[u].r)>>1; | |
ll res=0; | |
if(l<=mid) res=(res+query(ls, l, r))%mod; | |
if(r>mid) res=(res+query(rs, l, r))%mod; | |
return res; | |
} | |
int pr[100],tot; | |
void init() { | |
for(int i=2; i<=100; i++) { | |
int flag=0; | |
for(int j=2; j<=sqrt(i); j++) { | |
if(i%j==0) flag=1; | |
} | |
if(!flag) pr[++tot]=i; | |
} | |
for(int i=1; i<=tot; i++) { // 初始化 n 个位置的 fa 数组 | |
for(int j=1; j<=n; j++) { | |
if(arr[j]%pr[i]==0) fa[i][j]=j+1; | |
} | |
} | |
} | |
void mul(int u,int l,int r,int c) { // 线段树区间乘法 | |
if(tr[u].l>=l && tr[u].r<=r) { | |
lazy[u]=lazy[u]*c%mod; | |
tr[u].val=tr[u].val*c%mod; | |
} else { | |
pushdown(u); | |
int mid=tr[u].l+tr[u].r>>1; | |
if(l<=mid) mul(ls,l,r,c); | |
if(r>mid) mul(rs,l,r,c); | |
pushup(u); | |
} | |
} | |
int find(int c,int x) { | |
if(x==fa[c][x]) return x; | |
return fa[c][x]=find(c,fa[c][x]); | |
} | |
void merge(int l,int r,int c) { //[l,r] 区间乘以 pr [c] | |
int k; | |
for(int i=l; ; i=k+1) { | |
k=find(c,i); | |
if(k>r) { | |
mul(1,i,r,pr[c]); // 最后一个区间特判一下 | |
break; | |
} else { | |
if(k>i) mul(1,i,k-1,pr[c]); | |
fa[c][i]=k+1; // 并查集合并 | |
mul(1,k,k,pr[c]-1); // 单独相乘 | |
fa[c][k]=k+1; // 并查集合并 | |
} | |
} | |
} | |
void update(int l,int r,int w) { | |
for(int i=1; i<=25; i++) { // 第 i 个质因子有几个乘几次 | |
while(w%pr[i]==0) { | |
merge(l,r,i); | |
w/=pr[i]; | |
} | |
} | |
} | |
signed main() { | |
ios; | |
cin>>n>>m; | |
for(int i=1; i<=n; i++) { | |
cin>>arr[i]; | |
for(int j=1; j<=25; j++) fa[j][i]=i; | |
} | |
for(int j=1; j<=25; j++) fa[j][n+1]=n+1; // 因为 fa 数组指向的位置是 (最后一个包含质因子 c 的下一个位置),所以 n+1 也要进行初始化 | |
init(); | |
build(1,1,n); | |
while(m--) { | |
int op,l,r,w; | |
cin>>op; | |
if(op==0) { | |
cin>>l>>r>>w; | |
update(l,r,w); | |
} else { | |
cin>>l>>r; | |
cout<<query(1,l,r)<<endl; | |
} | |
} | |
return 0; | |
} | |
/* | |
5 5 | |
5 1 6 2 13 | |
0 5 5 25 | |
0 5 5 18 | |
1 3 5 | |
0 1 3 24 | |
1 3 4 | |
*/ |
# M Addition
可以借位的模拟题
#include <bits/stdc++.h> | |
#define endl '\n' | |
#define ios ios::sync_with_stdio(0);cin.tie(0);cout.tie(0) | |
using namespace std; | |
typedef long long ll; | |
const int N=100; | |
int s[N]; | |
int a[N],b[N]; | |
int ans[N]; | |
int main() { | |
ios; | |
int n; | |
cin>>n; | |
for(int i=1;i<=n;i++) cin>>s[i]; | |
for(int i=1;i<=n;i++) cin>>a[i]; | |
for(int i=1;i<=n;i++) cin>>b[i]; | |
int flag=0,pre; | |
for(int i=1;i<=n;i++){ | |
if(!flag){ | |
if(a[i] && b[i]){ | |
ans[i]=0; | |
flag=1; | |
pre=s[i]; | |
} | |
else if((a[i] && !b[i]) || (!a[i] && b[i])) ans[i]=1; | |
else ans[i]=0; | |
} | |
else{ | |
if(s[i]==pre){ | |
if(!a[i] && !b[i]) flag=0; | |
if((a[i] && b[i]) || (!a[i] && !b[i])) ans[i]=1; | |
else if((a[i] && !b[i]) || (!a[i] && b[i])) ans[i]=0; | |
} | |
else{ | |
if(a[i] || b[i]) flag=0; | |
if(a[i] && b[i]) ans[i]=1; | |
else if((a[i] && !b[i]) || (!a[i] && b[i])) ans[i]=0; | |
else ans[i]=1; | |
} | |
} | |
} | |
for(int i=1;i<=n;i++) { | |
if(i!=n) cout<<ans[i]<<" "; | |
else cout<<ans[i]<<endl; | |
} | |
return 0; | |
} | |
/* | |
*/ |