[luogu4140] 奇数国
题目
在一片美丽的大陆上有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] 奇数国的更多相关文章
- [BZOJ 3813]奇数国
3813: 奇数国 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 736 Solved: 416[Submit][Status][Discuss] ...
- 【BZOJ3813】奇数国 线段树+欧拉函数
[BZOJ3813]奇数国 Description 给定一个序列,每次改变一个位置的数,或是询问一段区间的数的乘积的phi值.每个数都可以表示成前60个质数的若干次方的乘积. Sample Input ...
- AC日记——【清华集训2014】奇数国 uoj 38
#38. [清华集训2014]奇数国 思路: 题目中的number与product不想冲: 即为number与product互素: 所以,求phi(product)即可: 除一个数等同于在模的意义下乘 ...
- 【bzoj3813】: 奇数国 数论-线段树-欧拉函数
[bzoj3813]: 奇数国 题意:给定一个序列,每个元素可以分解为最小的60个素数的形式.(x=p1^k1*p2^k2*......p60^k60)(p1=2,p2=3,…,p60=281) 支持 ...
- Bzoj 3813 奇数国 题解 数论+线段树+状压
3813: 奇数国 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 748 Solved: 425[Submit][Status][Discuss] ...
- [BZOJ3813] 奇数国 - 线段树
3813: 奇数国 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 912 Solved: 508[Submit][Status][Discuss] ...
- HYSBZ - 3813 奇数国 欧拉函数+树状数组(线段树)
HYSBZ - 3813奇数国 中文题,巨苟题,巨无敌苟!!首先是关于不相冲数,也就是互质数的处理,欧拉函数是可以求出互质数,但是这里的product非常大,最小都2100000,这是不可能实现的.所 ...
- BZOJ3813: 奇数国
传送门 欧拉函数+线段树 因为只有60个素数,所以把状态压成long long的形式.用线段树维护区间和和区间和中有多少个质数.然后xjb搞搞就行了,具体参见代码. //BZOJ 3813 //by ...
- 【线段树/数学/扩展欧几里得】 Bzoj 3913:奇数国
Description 在一片美丽的大陆上有100000个国家,记为1到100000.这里经济发达,有数不尽的账房,并且每个国家有一个银行.某大公司的领袖在这100000个银行开户时都存了3大洋,他惜 ...
随机推荐
- Java多线程原理+基础知识(超级超级详细)+(并发与并行)+(进程与线程)1
Java多线程 我们先来了解两个概念!!!! 1.什么是并发与并行 2.什么是进程与线程 1.什么是并发与并行 1.1并行:两个事情在同一时刻发生 1.2并发:两个事情在同一时间段内发生 并发与并行的 ...
- 更优雅的在 Xunit 中使用依赖注入
Xunit.DependencyInjection 7.0 发布了 Intro 上次我们已经介绍过一次大师的 Xunit.DependencyInjection https://www.cnblogs ...
- 使用truncate ,截断有外键约束的父表
此时有两种方法,解决1.删除外键约束,删除该表,在重建外键约束--查询外键约束select TABLE_NAME,CONSTRAINT_NAME,CONSTRAINT_TYPE,R_CONSTRAIN ...
- mdk/iar汇编区别
在代码移植中,经常遇到iar的代码转换问题,在此不间断记录一些,个人感觉还是IAR的更接近C一些,备查: 1. #ifdef的使用 // IAR #ifdef MACRO_XX #endif // M ...
- 谷歌浏览器又隐藏的HTTPS和WWW前缀
谷歌工程师 Emily Schechter 曾在 Chromium 反馈页面中提到:Chrome 团队一直将简易性.可用性.安全性作为衡量 UI 的标准.为了让 URL 能更好地被用户理解.移除那些容 ...
- java-把生成的随机数,指定范围(如:100-200),指定打印次数(如:50次),并进行去重。
package main.demo; public class Demo4 { /** * 随机指定范围内N个不重复的数 * 最简单最基本的方法 * @param min 指定范围最小值 * @par ...
- 使用ImpromptuInterface反射库方便的创建自定义DfaGraphWriter
在本文中,我为创建的自定义的DfaGraphWriter实现奠定了基础.DfaGraphWriter是公开的,因此您可以如上一篇文章中所示在应用程序中使用它,但它使用的所有类均已标记为internal ...
- vscode 无法自动补全第三方库
点击Settings 找到“Extentions”下的“Python”,点击“Auto Completes: Extra Paths”的“Edit in settings.json”,如下图: 在se ...
- 初步知道scss 简化css复杂层级
简介:今天在调试前端样式的时候,el-button组件位置需要调整并且 需要改字体大小 .直接上了一个 font-text:20px; 发现没作用,谷歌调试发现并未作用到组件里的<span> ...
- 启动扫描闪退,因为忘了在manifest里申请手机镜头使用许可了。
启动扫描闪退,因为忘了在manifest里申请手机镜头使用许可了.