bzoj 4026 dC Loves Number Theory
把我写吐了 太弱了
首先按照欧拉函数性质 我只需要统计区间不同质数个数就好了
一眼主席树
其次我被卡了分解质因数这里
可以通过质数筛时就建边解决
不够灵性啊,不知道如何改
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 1e6+1000;
const int H = 5e4+5;
const int M = H*100;
const int mod = 1e6+777;
int isprime[N];
int prime[N]; int prnum[N]; int cc = 0;
int rev[N], has[N];
struct Node{
int to,nex;
}E[N*3];
int head[N], eot = 0;
int n,q;
int a[H];
int T[H], ls[M], rs[M], v[M]; int tot;
int puf[H];
void add(int x,int y) {
E[eot].to = y; E[eot].nex = head[x]; head[x] = eot++;
}
int build(int l,int r) {
int rt = tot++;
v[rt] = 1;
if(l == r) return rt;
int mid = (l+r)>>1;
ls[rt] = build(l,mid); rs[rt] = build(mid+1,r);
return rt;
}
int upd(int pos,int num,int l,int r,int pr) {
int rt = tot++;
v[rt] = 1ll*v[pr]*num%mod;
if(l == r) return rt;
int mid = (l+r)>>1;
if(pos <= mid) {
ls[rt] = upd(pos,num,l,mid,ls[pr]); rs[rt] = rs[pr];
}else {
ls[rt] = ls[pr]; rs[rt] = upd(pos,num,mid+1,r,rs[pr]);
}
return rt;
}
int query(int L,int R,int l,int r,int rt) {
if(L <= l && r <= R) return v[rt];
int mid = (l+r)>>1;
int ans = 1;
if(L <= mid) ans = 1ll*ans*query(L,R,l,mid,ls[rt]) %mod;
if(R > mid) ans = 1ll*ans*query(L,R,mid+1,r,rs[rt]) %mod;
return ans;
}
void debug(int l,int r,int rt) {
printf("%d %d %d\n",l,r,v[rt]);
if(l == r) return;
int mid = (l+r)>>1;
debug(l,mid,ls[rt]); debug(mid+1, r, rs[rt]);
}
int main(){
memset(head,-1,sizeof(head));
for(int i = 2; i < N; ++i) {
if(!isprime[i]) {
prime[++cc] = i;
for(int j = i; j < N; j += i) {
isprime[j] = 1;
add(j,i);
}
}
}
rev[0] = rev[1] = 1;
for(int i=2;i< N;++i) {
rev[i]= 1ll*(mod-mod/i)* rev[mod%i]%mod;
}
for(int i = 1; i <= cc; ++i) {
prnum[prime[i]] = 1ll*(prime[i]-1) * rev[prime[i]] %mod;
}
while(~scanf("%d %d",&n,&q)) {
memset(has,0,sizeof(has));
for(int i = 1; i <= n; ++i) scanf("%d",&a[i]);
tot = 0; puf[0] = 1;
T[0] = build(1,n);
for(int i = 1; i <= n; ++i) {
puf[i] = 1ll*puf[i-1]*a[i]%mod;
int tmp = T[i-1];
int tt = 1;
for(int j = head[a[i]]; j != -1; j = E[j].nex) {
int x = E[j].to;
if(has[x])
tmp = upd(has[x],rev[prnum[x]],1,n,tmp);
tt = 1ll*tt*prnum[x]%mod;
has[x] = i;
}
T[i] = upd(i,tt,1,n,tmp);
}
int ans = 0;
for(int i = 0; i < q; ++i) {
int l,r; scanf("%d %d",&l,&r);
l^=ans; r^=ans;
ans = query(l,r,1,n,T[r]);
ans = 1ll*ans*puf[r]%mod *rev[puf[l-1]] %mod;
printf("%d\n",ans);
}
}
return 0;
}
bzoj 4026 dC Loves Number Theory的更多相关文章
- bzoj 4026 dC Loves Number Theory 主席树+欧拉函数
题目描述 dC 在秒了BZOJ 上所有的数论题后,感觉萌萌哒,想出了这么一道水题,来拯救日益枯竭的水题资源.给定一个长度为 n的正整数序列A,有q次询问,每次询问一段区间内所有元素乘积的φ(φ(n)代 ...
- BZOJ 4026 dC Loves Number Theory (主席树+数论+欧拉函数)
题目大意:给你一个序列,求出指定区间的(l<=i<=r) mod 1000777 的值 还复习了欧拉函数以及线性筛逆元 考虑欧拉函数的的性质,(l<=i<=r),等价于 (p[ ...
- BZOJ 4026: dC Loves Number Theory 可持久化线段树 + 欧拉函数 + 数学
Code: #include <bits/stdc++.h> #define ll long long #define maxn 50207 #define setIO(s) freope ...
- 【BZOJ4026】dC Loves Number Theory 分解质因数+主席树
[BZOJ4026]dC Loves Number Theory Description dC 在秒了BZOJ 上所有的数论题后,感觉萌萌哒,想出了这么一道水题,来拯救日益枯竭的水题资源. 给 ...
- BZOJ4026: dC Loves Number Theory
Description dC 在秒了BZOJ 上所有的数论题后,感觉萌萌哒,想出了这么一道水题,来拯救日益枯 竭的水题资源. 给定一个长度为 n的正整数序列A,有q次询问,每次询问一段区间内所 ...
- 【bzoj4026】dC Loves Number Theory 可持久化线段树
题目描述 dC 在秒了BZOJ 上所有的数论题后,感觉萌萌哒,想出了这么一道水题,来拯救日益枯竭的水题资源. 给定一个长度为 n的正整数序列A,有q次询问,每次询问一段区间内所有元素乘积的φ(φ(n ...
- [BZOJ4026]dC Loves Number Theory(线段树)
根据欧拉函数的定义式可知,可以先算出a[l]*a[l+1]*...*a[r]的值,然后枚举所有存在的质因子*(p-1)/p. 发现这里区间中一个质因子只要计算一次,所以指计算“上一个同色点在区间外”的 ...
- [BZOJ4026]dC Loves Number Theory 欧拉函数+线段树
链接 题意:给定长度为 \(n\) 的序列 A,每次求区间 \([l,r]\) 的乘积的欧拉函数 题解 考虑离线怎么搞,将询问按右端点排序,然后按顺序扫这个序列 对于每个 \(A_i\) ,枚举它的质 ...
- [bzoj4026]dC Loves Number Theory_主席树_质因数分解_欧拉函数
dC Loves Number Theory 题目大意:dC 在秒了BZOJ 上所有的数论题后,感觉萌萌哒,想出了这么一道水题,来拯救日益枯竭的水题资源. 给定一个长度为 n的正整数序列A,有q次询问 ...
随机推荐
- Shell与脚本
shell是Linux操作系统的用户接口,我们经常需要编写脚本让操作系统自动执行一系列指令的需求,本文将简单介绍开发shell脚本的所需的语言特性. shell脚本是指令序列,其指令可以直接在终端中执 ...
- POJ 3683 Priest John's Busiest Day[2-SAT 构造解]
题意: $n$对$couple$举行仪式,有两个时间段可以选择,问是否可以不冲突举行完,并求方案 两个时间段选择对应一真一假,对于有时间段冲突冲突的两人按照$2-SAT$的规则连边(把不冲突的时间段连 ...
- chrome_options
用法 from selenium.webdriver.chrome.options import Options chorme_option=Options() chorme_option.add ...
- Material使用11 核心模块和共享模块、 如何使用@angular/material
1 创建项目 1.1 版本说明 1.2 创建模块 1.2.1 核心模块 该模块只加载一次,主要存放一些核心的组件及服务 ng g m core 1.2.1.1 创建一些核心组件 页眉组件:header ...
- LeetCode - 492. Construct the Rectangle
For a web developer, it is very important to know how to design a web page's size. So, given a speci ...
- 腾讯IVWEB前端工程化工具feflow思考与实践
本篇文章主要介绍腾讯IVWEB团队从0到1在工程化的思考和实践.feflow的全称是Front-end flow(前端工作流),致力于提升研发效率和规范的工程化解决方案.愿景是通过feflow,可以使 ...
- MySQL事务处理案例演示
<?php $mysqli=new mysqli("localhost", "root", "123456", "xsphp ...
- 使用scp从远程服务器下载文件到本地
[下载远程文件到本地] scp -P 6008 root@192.168.1.123:/usr/data/1.zip /Users/abc/www [上传本地文件到远程] scp -P 6008 ...
- Python基础总结
刚学习Python时,边学边总结的,采用思维导图的形式, 适合回顾使用.内容参考<Python:从入门到实践>一书. 再给出一张Datacamp网站上的一张关于Python基础的总 ...
- Python自动化--语言基础7--操作日志、加密、发送邮件、线程、生产者消费者
1.操作日志 logging.basicConfig:日志的统一处理器,对日志的输出格式和方式做配置日志级别等级CRITICAL > ERROR > WARNING > INFO & ...