题面

传送门

思路

这题目是真的难读......阅读理解题啊......

但是理解了以后就发现,题目等价于:

给你一个区间,支持单点修改,以及查询一段区间的乘积的欧拉函数值,这个答案对19961993取模

这里是欧拉函数的原因显然,题目中的那个不相冲实际上就是扩展欧几里得里面的那个定理,要满足不相冲(也就是方程有解),$product$和$number$必须互质

序列当中,每个元素大小不超过1e6,质因数都是前60个

那么我们显然可以开一棵线段树来维护这个区间乘积,但是怎么处理欧拉函数呢?$O(\sqrt{n})$的复杂度求吗?但是这题可以到$1000000^{100000}$诶......

没关系,我们来看一个神秘小技巧

设一个数$x=\prod_{i=1}{k}p_i{a_i}$,那么:

$\varphi(x)=\prod_{i=1}{k}(p_i-1)p_i{a_i-1}=x\prod_{i=1}^{k}\frac{p_i-1}{p_i}$

那么我们再开一棵线段树,把60个质因数在对应区间里的出现情况压进一个long long里面

每次查询的时候,查询出来取模过的乘积,再对每个出现过的质因数乘上模意义下的$\frac{p_i-1}{p_i}$,就是答案了

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cassert>
#define ll long long
#define mp make_pair
using namespace std;
inline int read(){
int re=0,flag=1;char ch=getchar();
while(ch>'9'||ch<'0'){
if(ch=='-') flag=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9') re=(re<<1)+(re<<3)+ch-'0',ch=getchar();
return re*flag;
}
const ll MOD=19961993;
ll qpow(ll a,ll b){
ll re=1ll;
while(b){
if(b&1) re=re*a%MOD;
a=a*a%MOD;b>>=1;
}
return re;
}
int vis[310],pri[70],cntp,inv[70];
void init(){
int i,j,k;vis[1]=1;
for(i=2;i<=281;i++){
if(!vis[i]) pri[++cntp]=i,inv[cntp]=qpow(i,MOD-2);
for(j=1;j<=cntp;j++){
k=i*pri[j];if(k>281) break;
vis[k]=1;
if(i%pri[j]==0) break;
}
}
}
ll a[400010],bit[400010];//a是乘积,b是压位的质因数状态
void update(int num){
int son=num<<1;
a[num]=a[son]*a[son+1]%MOD;
bit[num]=bit[son]|bit[son+1];
}
void build(int l,int r,int num){
int mid=(l+r)>>1;
if(l==r){
a[num]=3;bit[num]=2;return;
}
build(l,mid,num<<1);build(mid+1,r,(num<<1)+1);
update(num);
}
void change(int l,int r,int num,int pos,ll val){
int mid=(l+r)>>1,i;
if(l==r){
a[num]=val;bit[num]=0;
for(i=1;i<=60;i++) if(val%pri[i]==0) bit[num]|=(1ll<<(i-1));
return;
}
if(mid>=pos) change(l,mid,num<<1,pos,val);
else change(mid+1,r,(num<<1)+1,pos,val);
update(num);
}
pair<ll,ll> query(int l,int r,int ql,int qr,int num){
int mid=(l+r)>>1;pair<ll,ll>re=mp(1,0),tmp;
if(l>=ql&&r<=qr) return mp(a[num],bit[num]);
if(mid>=ql){
tmp=query(l,mid,ql,qr,num<<1);
re.first=re.first*tmp.first%MOD;
re.second|=tmp.second;
}
if(mid<qr){
tmp=query(mid+1,r,ql,qr,(num<<1)+1);
re.first=re.first*tmp.first%MOD;
re.second|=tmp.second;
}
return re;
}
int main(){
int n=read(),i,t1,t2,t3;build(1,100000,1);pair<ll,ll>tmp;
init();
while(n--){
t1=read();t2=read();t3=read();
if(t1) change(1,100000,1,t2,t3);
else{
tmp=query(1,100000,t2,t3,1);
for(i=1;i<=60;i++)
if(tmp.second&(1ll<<(i-1)))
tmp.first=tmp.first*(pri[i]-1)%MOD*inv[i]%MOD;
printf("%lld\n",tmp.first);
}
}
}

[bzoj3813] 奇数国 [线段树+欧拉函数]的更多相关文章

  1. 【BZOJ3813】奇数国 线段树+欧拉函数

    [BZOJ3813]奇数国 Description 给定一个序列,每次改变一个位置的数,或是询问一段区间的数的乘积的phi值.每个数都可以表示成前60个质数的若干次方的乘积. Sample Input ...

  2. BZOJ 3813--奇数国(线段树&欧拉函数&乘法逆元&状态压缩)

    3813: 奇数国 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 755  Solved: 432[Submit][Status][Discuss] ...

  3. 【bzoj3813】: 奇数国 数论-线段树-欧拉函数

    [bzoj3813]: 奇数国 题意:给定一个序列,每个元素可以分解为最小的60个素数的形式.(x=p1^k1*p2^k2*......p60^k60)(p1=2,p2=3,…,p60=281) 支持 ...

  4. [BZOJ3813] 奇数国 - 线段树

    3813: 奇数国 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 912  Solved: 508[Submit][Status][Discuss] ...

  5. Please, another Queries on Array?(Codeforces Round #538 (Div. 2)F+线段树+欧拉函数+bitset)

    题目链接 传送门 题面 思路 设\(x=\prod\limits_{i=l}^{r}a_i\)=\(\prod\limits_{i=1}^{n}p_i^{c_i}\) 由欧拉函数是积性函数得: \[ ...

  6. 线段树+欧拉函数——cf1114F

    调了半天,写线段树老是写炸 /* 两个操作 1.区间乘法 2.区间乘积询问欧拉函数 欧拉函数计算公式 phi(mul(ai))=mul(ai) * (p1-1)/p1 * (p2-1)/p2 * .. ...

  7. Please, another Queries on Array? CodeForces - 1114F (线段树,欧拉函数)

    这题刚开始看成求区间$\phi$和了........先说一下区间和的做法吧...... 就是说将题目的操作2改为求$(\sum\limits_{i=l}^{r}\phi(a[i]))\%P$ 首先要知 ...

  8. BZOJ4869 六省联考2017相逢是问候(线段树+欧拉函数)

    由扩展欧拉定理,a^(a^(a^(……^x)))%p中x作为指数的模数应该是φ(φ(φ(φ(……p)))),而p取log次φ就会变为1,也即每个位置一旦被修改一定次数后就会变为定值.线段树维护区间剩余 ...

  9. BZOJ 4026: dC Loves Number Theory 可持久化线段树 + 欧拉函数 + 数学

    Code: #include <bits/stdc++.h> #define ll long long #define maxn 50207 #define setIO(s) freope ...

随机推荐

  1. java后台poi根据模板导出excel

    public class ExcelUtils { private static final String INSPECTIONRECORD_SURFACE_TEMPLET_PATH = " ...

  2. Apache Maven(一):快速入门

    Maven 是什么? Maven 是一个项目管理和整合工具.Maven 为开发者提供了一套完整的构建生命周期框架.开发团队几乎不用花多少时间就能够自动完成工程的基础构建配置,因为 Maven 使用了一 ...

  3. Pandas 数据结构Series:基本概念及创建

    Series:"一维数组" 1. 和一维数组的区别 # Series 数据结构 # Series 是带有标签的一维数组,可以保存任何数据类型(整数,字符串,浮点数,Python对象 ...

  4. python——集合

    在python中,字典的亲戚就是集合,集合就是无映射关系的字典,花括号并不是字典的特权.如下面程序所示: >>> num = {} >>> type(num) &l ...

  5. Balance POJ - 1837

    Description Gigel has a strange "balance" and he wants to poise it. Actually, the device i ...

  6. 大话循环神经网络(RNN)

      在上一篇文章中,介绍了 卷积神经网络(CNN)的算法原理,CNN在图像识别中有着强大.广泛的应用,但有一些场景用CNN却无法得到有效地解决,例如: 语音识别,要按顺序处理每一帧的声音信息,有些结果 ...

  7. Apache 多端口配置方法

    首先修改httpd.conf配置文件. 添加8080端口 Listen 8080 打开虚拟配置文件 # Virtual hosts Include conf/extra/httpd-vhosts.co ...

  8. DOS程序员手册(一)

    当今MS-Windows横扫大江南北,让我们这就来研究一下它的祖宗——MS-DOS! 这本书很难得,希望读者好好学习! DOS程序员手册(一) DOS教程 (以下内容全部为原作者的阐述,照样保留) 这 ...

  9. Pascal小游戏 随机函数

    一个被人写滥了的小程序,新手学习,Pascal By Chaobs 初学者可以用它来学习随机函数的运用,当然你完全可以自己写一个随机函数. var   player1,player2:longint; ...

  10. USACO Section1.1 Broken Necklace 解题报告

    beads解题报告 —— icedream61 博客园(转载请注明出处)---------------------------------------------------------------- ...