考试题 T1
题意分析
就是让你求
\]
那么接下来就是化简式子
\]
\]
\]
我们考虑枚举\(x\)的因子\(k=dy\)
那么贡献就是\(μ(\frac{k}{y})*s[i]\)
其中\(s[i]=\sum_{i|d}val[d]\)
这里直接让\(val[a[i]]=val[i]\)了
然后查询修改我们都可以\(\sqrt n\)了
CODE:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<string>
#include<queue>
#include<map>
#include<stack>
#include<list>
#include<set>
#include<deque>
#include<vector>
#include<ctime>
#define ll long long
#define inf 0x7fffffff
#define N 10000008
#define IL inline
#define M 1008611
#define D double
#define maxn 10000000
#define mod 20020303
#define R register
using namespace std;
template<typename T>IL void read(T &_)
{
T __=0,___=1;char ____=getchar();
while(!isdigit(____)) {if(____=='-') ___=0;____=getchar();}
while(isdigit(____)) {__=(__<<1)+(__<<3)+____-'0';____=getchar();}
_=___ ? __:-__;
}
/*-------------OI使我快乐-------------*/
int n,m,tot;
struct Node
{
int xx,yy;
}e[M];
int prime[N],mul[N],val[N];
ll sum[N];
bool mark[N];
IL void work()
{
mul[1]=1;
for(R int i=2;i<=maxn;++i)
{
if(!mark[i]) {prime[++tot]=i;mul[i]=-1;}
for(R int j=1;j<=tot&&prime[j]*i<=maxn;++j)
{
mark[prime[j]*i]=1;
if(i%prime[j]==0)
{
mul[prime[j]*i]=0;
break;
}
else mul[prime[j]*i]=-mul[i];
}
}
}
int main()
{
freopen("t1.in","r",stdin);
freopen("t1.out","w",stdout);
read(n);read(m);work();
for(R int i=1,x,y;i<=n;++i)
{
read(x);read(y);
val[x]=y;
}
for(R int i=1;i<=maxn;++i)
for(R int j=i;j<=maxn;j+=i)
sum[i]+=val[j];
while(m--)
{
int knd,x,y;ll ans=0;
read(knd);read(x);read(y);
if(knd==1)
{
for(R int i=1;i*i<=x;++i)
{
if(x%i) continue;
int cdy=i,wzy=x/i;
if(cdy%y==0)
{
int now=cdy/y;
ans=(ans+mul[now]*sum[cdy]%mod+mod)%mod;
}
if(cdy==wzy) continue;
if(wzy%y==0)
{
int now=wzy/y;
ans=(ans+mul[now]*sum[wzy]%mod+mod)%mod;
}
}
printf("%lld\n",ans);
}
else
{
if(val[x])
{
for(R int i=1;i*i<=x;++i)
{
if(x%i) continue;
int cdy=i,wzy=x/i;
sum[cdy]-=val[x];
if(cdy==wzy) continue;
sum[wzy]-=val[x];
}
val[x]=0;
}
else
{
val[x]=y;
for(R int i=1;i*i<=x;++i)
{
if(x%i) continue;
int cdy=i,wzy=x/i;
sum[cdy]+=val[x];
if(cdy==wzy) continue;
sum[wzy]+=val[x];
}
}
}
}
fclose(stdin);
fclose(stdout);
return 0;
}
HEOI 2019 RP++
考试题 T1的更多相关文章
- kma 2019CSP前刷题记录
2019/10/25 \([LNOI2014]\ LCA\) \([Luogu\ P2774]\) 方格取数问题 \(Gauss\)消元板 \([JSOI2008]\)球形空间产生器 2019/10/ ...
- T1加权像(T1 weighted image,T1WI)
T1加权成像(T1-weighted imaging,T1WI)是指这种成像方法重点突出组织纵向弛豫差别,而尽量减少组织其他特性如横向弛豫等对图像的影响. 弛豫:物理用语,从某一个状态恢复到平衡态的过 ...
- 关于2016.12.12——T1的反思:凸包的意义与应用
2016.12.12 T1 给n个圆,保证圆圆相离,求将圆围起来的最小周长.n<=100 就像上图.考场上,我就想用切线的角度来做凸包.以圆心x,y排序,像点凸包一样,不过用两圆之间的下切线角度 ...
- T2 Func<in T1,out T2>(T1 arg)
委托调用方法的4种方式. using System; using System.Collections.Generic; namespace ConsoleApplication1 { delegat ...
- E1、T1链路
北美的24路脉码调制PCM简称T1 速率是1.544Mbit/s 北美使用的T1系统共有24个话路,每个话路采样脉冲用7bit编码,然后再加上1位信令码元,因此一个话路占用8bit. 帧同步码是在24 ...
- Action<T1, T2>委托
封装包含两个参数的方法委托,没有返回值. 语法 public delegate void Action<in T1, in T2>( T1 arg1, T2 arg2 ) 类型参数 in ...
- 有三个线程T1 T2 T3,如何保证他们按顺序执行-转载
T3先执行,在T3的run中,调用t2.join,让t2执行完成后再执行t3 在T2的run中,调用t1.join,让t1执行完成后再让T2执行 public class Test { // 1.现在 ...
- 现在有T1、T2、T3三个线程,怎样保证T2在T1执行完后执行,T3在T2执行完后执行?使用Join
public class TestJoin { public static void main(String[] args) { Thread t1 = new Thread(new T1(), &q ...
- 【测试】在hr用户下自行创建T1和T2表写一条SQL语句,(NL连接)
SQL> select t1.* from t1,t2 where t1.object_id=t2.object_id; rows selected. Execution Plan ------ ...
随机推荐
- MYSQL 存储过程通用
返回随机时间 函数 )) ) CHARSET utf8 BEGIN )); -- 随机天数 60天以内随机天数 )); -- 随机小时 ));-- 随机分 ));-- 随机秒 ); IF type = ...
- c++静态变量与菲静态变量
刚开始用C++写程序,遇到了一个很奇怪的问题,就是在类定义的面定义了一个静态成员变量,但在使用时编译器报错称变量不是类的成员. 文件a.h内容: /**********************A.h* ...
- Windows pip安装失败:no module named pkg_resources
通常是Setuptools安装出错,下载以下ez_setup.py文件后,先执行:ez_setup.py -U setuptools 重新安装setuptools 通过此ez_setup.py pip ...
- [OS]windows 2012 server-Local users and groups-Backup Operators
怎样找到windows 2012 server上的Backup Operators Press the Windows + R keys to open the Run dialog, type lu ...
- 10个实用的UX设计作品推销小窍门
以下内容由Mockplus(摹客)团队翻译整理,仅供学习交流,Mockplus是更快更简单的原型设计工具. 众所周知,产品用户体验很重要,即使是Google也知道这一点.但是,当真正涉及到UX 设计销 ...
- Android内核的编译与测试
1.下载Android内核 source.android.com/source->Downloading and Building Building Kernels 大概要花2个小时,其源码在培 ...
- C语言条件编译及编译预处理阶段(转)
一.C语言由源代码生成的各阶段如下: C源程序->编译预处理->编译->优化程序->汇编程序->链接程序->可执行文件 其中 编译预处理阶段,读取c源程序,对其中的 ...
- B-spline Curves 学习之B样条曲线的导数(8)
Derivatives of a B-spline Curve 本博客转自前人的博客的翻译版本,前几章节是原来博主的翻译内容,但是后续章节博主不在提供翻译,后续章节我在完成相关的翻译学习. (原来博客 ...
- 基于httpcore(httpclient component)搭建轻量级http服务器
下面是apache官网例子 服务器端接受请求,实现接收文件请求处理器 import java.io.File; import java.io.IOException; import java.io.I ...
- shell 脚本 查看班上同学的网络状态
#!/bin/base # a=192.168.100. //定义变量 c=('王浩' '谢云生' '黄科杨' '何星宇' '张宸兵' '胡燕' '刘桃') //定义数组 for b in {101. ...