题目

在一片美丽的大陆上有100000个国家,记为1到100000。这里经济发达,有数不尽的账房,并且每个国家有一个银行。某大公司的领袖在这100000个银行开户时都存了3大洋,他惜财如命,因此会不时地派小弟GFS清点一些银行的存款或者让GFS改变某个银行的存款。该村子在财产上的求和运算等同于我们的乘法运算,也就是说领袖开户时的存款总和为3100000。这里发行的软妹面额是最小的60个素数(p1=2,p2=3,…,p60=281),任何人的财产都只能由这60个基本面额表示,即设某个人的财产为fortune(正整数),则fortune=p1^k1*p2^k2*……p60^K60。

领袖习惯将一段编号连续的银行里的存款拿到一个账房去清点,为了避免GFS串通账房叛变,所以他不会每次都选择同一个账房。GFS跟随领袖多年已经摸清了门路,知道领袖选择账房的方式。如果领袖选择清点编号在[a,b]内的银行财产,他会先对[a,b]的财产求和(计为product),然后在编号属于[1,product]的账房中选择一个去清点存款,检验自己计算是否正确同时也检验账房与GFS是否有勾结。GFS发现如果某个账房的编号number与product相冲,领袖绝对不会选择这个账房。怎样才算与product不相冲呢?若存在整数x,y使得number*x+product*y=1,那么我们称number与product不相冲,即该账房有可能被领袖相中。当领袖又赚大钱了的时候,他会在某个银行改变存款,这样一来相同区间的银行在不同的时候算出来的product可能是不一样的,而且领袖不会在某个银行的存款总数超过1000000。

现在GFS预先知道了领袖的清点存款与变动存款的计划,想请你告诉他,每次清点存款时领袖有多少个账房可以供他选择,当然这个值可能非常大,GFS只想知道对19961993取模后的答案。

题解

首先看到题目要求$number*x+product*y=1$

这其实就是一个二元一次不定方程

在学习扩展欧几里得的时候我们知道

如果a、b是整数,那么一定存在整数x、y使得 $ax+by=gcd(a,b)$

那么$ax+by=1$有解当且仅当a,b互质,即要满足number与product互质

也就是说对于一次询问,求得就是product以内有多少与它互质的数,这个就是欧拉函数

但是product太大了,无法直接求,注意到题目说product可以表示成$p1^{k1}*p2^{k2}*……p60^{K60}$

我们知道欧拉函数有三条性质

1. 若a为质数,则$ψ(a)=a-1$

2. 若a,b互质,则$ψ(a*b)$=$ψ(a)*ψ(b)$

3. 若a为质数,b为a的倍数,则$ψ(a*b)=ψ(b)*a$

那么

$ψ(p1^{k1}*p2^{k2}*……p60^{K60})$

$=ψ(p1^{k1})*ψ(p2^{k2})……$(性质2)

$=ψ(p1)*p1^{k1-1}*ψ(p2)*p2^{k2-1}……$(性质3)

$=(p1-1)*p1^{k1-1}*(p2-1)*p2^{k2-1}……$(性质1)

$={(p1^{k1}*p2^{k2}*……)}/(p1*p2*……)*(p1-1)*(p2-1)…… $

$=product/(p1*p2*……)*(p1-1)*(p2-1)…… $

上式中的p都是被product包含的

除法可以用逆元处理

那么,我们可以用线段树维护每段区间的product(取模后),每个质数是否出现(可以用一个long long或者数组存储),询问时根据返回的数组计算后面那一坨就行了

代码

注意线段树右边界要代100000进去,而不是n(我在这里调了0.5h)

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
#define int long long
#define N 500010
#define lc id*2
#define rc id*2+1
#define mid (l+r)/2
#define mod 19961993
bool fact[N][61];
int prim[61]={0,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97,101,103,107,109,113,127,131,137,139,149,151,157,163,167,173,179,181,191,193,197,199,211,223,227,229,233,239,241,251,257,263,269,271,277,281};
int inv[61];
int sum[N];
void build(int id,int l,int r)
{
fact[id][2]=1;
if(l==r)
{
sum[id]=3;
return;
}
build(lc,l,mid);
build(rc,mid+1,r);
sum[id]=sum[lc]*sum[rc]%mod;
}
void modify(int id,int l,int r,int pos,int val)
{
if(l==r)
{
memset(fact[id],0,sizeof(fact[id]));
sum[id]=val;
for(int i=1;i<=60;i++)
if(val%prim[i]==0) fact[id][i]=1;
return;
}
if(pos<=mid) modify(lc,l,mid,pos,val);
else modify(rc,mid+1,r,pos,val);
for(int i=1;i<=60;i++) fact[id][i]=fact[lc][i]|fact[rc][i];
sum[id]=sum[lc]*sum[rc]%mod;
}
int query(int id,int l,int r,int tl,int tr,bool arr[])
{
if(l>=tl&&r<=tr)
{
for(int i=1;i<=60;i++) arr[i]|=fact[id][i];
return sum[id];
}
int ret=1;
if(tl<=mid) ret=query(lc,l,mid,tl,tr,arr);
if(tr>mid) ret*=query(rc,mid+1,r,tl,tr,arr);
return ret%mod;
}
int qpow(int a,int b)
{
int ans=1;
while(b)
{
if(b&1) ans*=a,ans%=mod;
a*=a;
a%=mod;
b>>=1;
}
return ans;
}
signed main()
{
int n;
cin>>n;
build(1,1,100000);
for(int i=1;i<=60;i++) inv[i]=qpow(prim[i],mod-2);
for(int i=1;i<=n;i++)
{
int a,b,c;
scanf("%lld%lld%lld",&a,&b,&c);
if(a==0)
{
bool arr[61]={0};
int tot=1,tot2=1,tot3=query(1,1,100000,b,c,arr);
for(int j=1;j<=60;j++) if(arr[j]) tot*=inv[j],tot2*=(prim[j]-1),tot%=mod,tot2%=mod;
printf("%lld\n", tot3*tot2%mod*tot%mod);
}
else modify(1,1,100000,b,c);
}
}

  

[luogu4140] 奇数国的更多相关文章

  1. [BZOJ 3813]奇数国

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

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

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

  3. AC日记——【清华集训2014】奇数国 uoj 38

    #38. [清华集训2014]奇数国 思路: 题目中的number与product不想冲: 即为number与product互素: 所以,求phi(product)即可: 除一个数等同于在模的意义下乘 ...

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

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

  5. Bzoj 3813 奇数国 题解 数论+线段树+状压

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

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

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

  7. HYSBZ - 3813 奇数国 欧拉函数+树状数组(线段树)

    HYSBZ - 3813奇数国 中文题,巨苟题,巨无敌苟!!首先是关于不相冲数,也就是互质数的处理,欧拉函数是可以求出互质数,但是这里的product非常大,最小都2100000,这是不可能实现的.所 ...

  8. BZOJ3813: 奇数国

    传送门 欧拉函数+线段树 因为只有60个素数,所以把状态压成long long的形式.用线段树维护区间和和区间和中有多少个质数.然后xjb搞搞就行了,具体参见代码. //BZOJ 3813 //by ...

  9. 【线段树/数学/扩展欧几里得】 Bzoj 3913:奇数国

    Description 在一片美丽的大陆上有100000个国家,记为1到100000.这里经济发达,有数不尽的账房,并且每个国家有一个银行.某大公司的领袖在这100000个银行开户时都存了3大洋,他惜 ...

随机推荐

  1. 简单实用的办公软件导航网站,IT经理必备工具

    最近非常忙,因为公司上线了业财一体化系统.今天分享一个非常实用的办公软件导航网站,节省了我很多百度的时间. 快氪导航,让软件服务更简单. 一.办公软件导航 站长已经按照功能进行了分类:协同办公,流程审 ...

  2. 搭建sonarqube分析golang代码

    准备postgres 为什么不使用Mysql呢,因为从7.9就不支持了 docker启动postgres docker run -d \ --name sonar-postgres \ -e POST ...

  3. 9.CSMA_CD协议

    先听再说,边听边说 载波监听多点接入/碰撞检测CSMA/CD( carrier sense multiple access with collision detection) CD:碰撞检测(冲突检测 ...

  4. 怎么训练出一个NB的Prophet模型

    上篇<神器の争>主要是介绍Prophet的特点以及prophet入门的一些注意事项,但离真正的实际运用还有段距离.本篇主要讲解实际运用中Prophet调参的主要步骤以及一些本人实际经验. ...

  5. 学会Python除了不能生孩子,其他的都能做。

    随着人工智能的迅猛发展,相信大家对于it行业最熟悉的词莫过于 Python.那么,Python究竟可以做些什么呢?一个资深程序员说:“学会Python除了不能生孩子,其他的都能做.”加入3137821 ...

  6. 性能测试 -- docker安装influxdb

    一.前提 1.项目已经部署好 2.docker已经安装好 二.docker安装influxdb 1.下载influxdb镜像:docker pull tutum/influxdb    1)超时报错: ...

  7. sublime text 配置文件

    不慌. 上网上下的clang全是报错,所以就想先配置一下文件. 其实它本来的配置我觉得挺好的. 可能因人而异有些地方还是改了改. 原版有英语注释的,就当是学英语了 // While you can e ...

  8. Text多行文本框基本用法

    1.Text(root,width,height,bg) 主窗口,宽度,高度,背景色 2.使用 .insert() 方法添加内容 Text 对象.insert(几行.几列,"内容" ...

  9. PHP xpath() 函数

    定义和用法 xpath() 函数运行对 XML 文档的 XPath 查询.高佣联盟 www.cgewang.com 如果成功,该函数返回 SimpleXMLElements 对象的一个数组.如果失败, ...

  10. Canvas知识点补充

    Canvas笔记 复习 初识canvas <canvas> 是 HTML5 新增的,一个可以使用脚本(通常为 JavaScript) 在其中绘制图像的 HTML 元素.它可以用来制作照片集 ...