P2350-[HAOI2012]外星人【线性筛】

正题

题目链接:https://www.luogu.com.cn/problem/P2350


题目大意

给出NNN质因数分解之后的结果,求每次N=φ(N)N=\varphi(N)N=φ(N),多少次后N=1N=1N=1

N=∏i=1mpiqi,1≤m≤2000,1≤pi≤105,1≤qi≤109N=\prod_{i=1}^mp_i^{q_i},1\leq m\leq 2000,1\leq p_i\leq 10^5,1\leq q_i\leq 10^9N=i=1mpiqi,1m2000,1pi105,1qi109


解题思路

我是傻逼((((

首先φ(∏i=1mpiqi)=∏i=1m(pi−1)piqi−1\varphi(\prod_{i=1}^mp_i^{q_i})=\prod_{i=1}^m(p_i-1)p_i^{q_i-1}φ(i=1mpiqi)=i=1m(pi1)piqi1
开始想太多了,以为传递的延迟影响很大,后来发现设因为pi−1p_i-1pi1一定是偶数,所以把pi−1p_i-1pi1分解后一定有222这个质因子,所以只有在NNN没有222这个质因子的情况下传递的延迟会造成一次的影响,否则都没有影响,因为222还没有抵消延迟就已经结束了。

所以我们只需要考虑每个质数能传递出多少个222,这个用线性筛可以解决因为最后更新一个质数的一定是新增的一个质因子。


code

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1e5+10,M=1e5+10;
int T,n,cnt,pri[N],f[M],v[M];
void Prime(){f[1]=1;for(int i=2;i<M;i++){if(!v[i])pri[++cnt]=i,f[i]=f[i-1];for(int j=1;j<=cnt&&i*pri[j]<M;j++){v[i*pri[j]]=1;f[i*pri[j]]=f[i]+f[pri[j]];if(i%pri[j]==0)break;}}return;
}
int main()
{Prime();scanf("%d",&T);while(T--){scanf("%d",&n);if(!n)puts("0");long long ans=1;for(int i=1;i<=n;i++){int x,w;scanf("%d%d",&x,&w);if(x==2)ans--;ans+=1ll*f[x]*w;}printf("%lld\n",ans);}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318773.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

牛客题霸 [ 最长回文子串] C++题解/答案

牛客题霸 [ 最长回文子串] C题解/答案 题目描述 对于一个字符串&#xff0c;请设计一个高效算法&#xff0c;计算其中最长回文子串的长度。 给定字符串A以及它的长度n&#xff0c;请返回最长回文子串的长度。 题解&#xff1a; 两个方法&#xff1a; 一个是经典暴力&#…

强连通分量:洛谷P3387 模板:缩点

传送门 顾名思义&#xff0c;模板awa #include <cstdio> #include <cstring> #include <cmath> #include <algorithm> #include <iostream> #include <string> #include <queue> #include <vector> #define ll long long usi…

2021牛客暑期多校训练营4 E-Tree Xor(异或+思维+区间交 or Trie树)

E-Tree Xor 首先不考虑区间限制条件&#xff0c;我们给定其中一个点的权值后&#xff0c;那么其他点的权值也就确定。比如 val10\text{val}_10val1​0&#xff0c;即可通过变得限制求出其他点valu\text{val}_uvalu​&#xff0c;而且不难发现如果val10⊕a\text{val}_10\oplus …

【.NET Core项目实战-统一认证平台】第八章 授权篇-IdentityServer4源码分析

上篇文章我介绍了如何在网关上实现客户端自定义限流功能&#xff0c;基本完成了关于网关的一些自定义扩展需求&#xff0c;后面几篇将介绍基于IdentityServer4&#xff08;后面简称Ids4&#xff09;的认证相关知识&#xff0c;在具体介绍ids4实现我们统一认证的相关功能前&…

[科技]Loj#6564-最长公共子序列【bitset】

正题 题目链接:https://loj.ac/p/6564 题目大意 给两个序列a,ba,ba,b求它们的最长公共子序列。 1≤n,m,ai,bi≤71041\leq n,m,a_i,b_i\leq 7\times 10^41≤n,m,ai​,bi​≤7104 解题思路 无意间看到的一个bitsetbitsetbitset科技。 首先设fi,jf_{i,j}fi,j​表示aaa串匹配到…

牛客题霸 [找到字符串的最长无重复字符子串] C++题解/答案

牛客题霸 [找到字符串的最长无重复字符子串] C题解/答案 题目描述 给定一个数组arr&#xff0c;返回arr的最长无的重复子串的长度(无重复指的是所有数字都不相同)。 题解&#xff1a; i和j两个指针分别指不重复子串的两端&#xff0c;根据情况各自移动&#xff0c;记录最大…

并查集:P1196 [NOI2002] 银河英雄传说

传送门 带权并查集最棒了 dep记录深度&#xff0c;size记录集合大小 转移见代码 #include<cstdio> #include<algorithm> #include<iostream> using namespace std; const int N31000; int n,m; int t; int dep[N],fa[N],size[N]; char s; int a,b; int find…

2021牛客暑期多校训练营4 D-Rebuild Tree(prufer序列+树形dp)

D-Rebuild Tree Prufer 是这样建立的&#xff1a;每次选择一个编号最小的叶结点并删掉它&#xff0c;然后在序列中记录下它连接到的那个结点。重复n−2n-2n−2次后就只剩下两个结点&#xff0c;算法结束。&#xff08;为什么不是n−1n-1n−1次呢&#xff1f;因为第n−1n-1n−1…

牛客题霸 [ 求二叉树的层序遍历] C++题解/答案

牛客题霸 [ 求二叉树的层序遍历] C题解/答案 题目描述 给定一个二叉树&#xff0c;返回该二叉树层序遍历的结果&#xff0c;&#xff08;从左到右&#xff0c;一层一层地遍历&#xff09; 例如&#xff1a; 给定的二叉树是{3,9,20,#,#,15,7}, 该二叉树层序遍历的结果是 [ [3…

T183637-变异距离(2021 CoE III C)【单调栈】

正题 题目链接:https://www.luogu.com.cn/problem/T183637 题目大意 给出nnn个二元组(xi,yi)(x_i,y_i)(xi​,yi​)&#xff0c;求最大的 ∣xi−xj∣min{∣yi∣,∣yj∣}|x_i-x_j|\times min\{|y_i|,|y_j|\}∣xi​−xj​∣min{∣yi​∣,∣yj​∣} 1≤n≤2106,−106≤xi≤106,−…

dfs剪枝:洛谷P2809 hzwer爱折纸

传送门 解析 dfs暴力枚举即可 这题的重点是如何剪枝 不难发现&#xff0c;随着不断处理&#xff0c;纸条只会越来越短&#xff0c;且所有数字总加和不变 我一开始想到了2个条件&#xff1a; 1.当前长度比理想纸条小&#xff0c;return&#xff1b; 2.总加和与理想纸条不等&a…

2021牛客暑期多校训练营2 B-Cannon(组合+推式子)

B-Cannon 首先nnn个炮在一行操作一次的方案数为2(n−2)2(n-2)2(n−2):前面两个炮只能向右吃&#xff0c;最后两个跑只能向左吃&#xff0c;而其余的炮既可以向左也可以向右&#xff0c;于是有42(n−4)42(n-4)42(n−4)种 于是操作mmm次的操作排列的方案数为2m(n−2)(n−3)...(…

[译]聊聊C#中的泛型的使用

写在前面今天忙里偷闲在浏览外文的时候看到一篇讲C#中泛型的使用的文章&#xff0c;因此加上本人的理解以及四级没过的英语水平斗胆给大伙进行了翻译&#xff0c;当然在翻译的过程中发现了一些问题&#xff0c;因此也进行了纠正&#xff0c;当然&#xff0c;原文的地址我放在最…

牛客题霸 [ 最长递增子序列] C++题解/答案

牛客题霸 [ 最长递增子序列] C题解/答案 题目描述 给定数组arr&#xff0c;设长度为n&#xff0c;输出arr的最长递增子序列。&#xff08;如果有多个答案&#xff0c;请输出其中字典序最小的&#xff09; 题意&#xff1a; 直接暴力会超时 应该用二分贪心 题解&#xff1a…

P5405-[CTS2019]氪金手游【树形dp,容斥,数学期望】

前言 话说在LojLojLoj下了个数据发现这题的名字叫fgofgofgo 正题 题目链接:https://www.luogu.com.cn/problem/P5405 题目大意 nnn张卡的权值为1/2/31/2/31/2/3的概率权重分别是px,1/2/3p_{x,1/2/3}px,1/2/3​&#xff0c;然后按照权值每次获得一张未获得的卡&#xff0c;然后…

模板:线段树

文章目录引言思想模板建树单点修改 / 查询区间修改/查询总结练习引言 有一类题目:要求在区间上维护信息&#xff0c;比如带修改区间求和问题。考虑到枚举求和肯定会超时&#xff0c;我们可以通过一些数据结构来维护信息&#xff0c;例如线段树。 它功能强大&#xff0c;支持区…

Docker最全教程——从理论到实战(三)

容器是应用走向云端之后必然的发展趋势&#xff0c;因此笔者非常乐于和大家分享我们这段时间对容器的理解、心得和实践。本篇教程持续编写了2个星期左右&#xff0c;只是为了大家更好地了解、理解和消化这个技术&#xff0c;能够搭上这波车。你可以关注我们的公众号“magiccode…

2021牛客暑期多校训练营4 H-Convolution(数学)

H-Convolution x⨂yxygcd⁡(x,y)x\bigotimes y\frac{xy}{\gcd(x,y)}x⨂ygcd(x,y)xy​ 下面大佬题解推&#x1f981; Ultraman-Ace题解 #include<bits/stdc.h> using namespace std; using lllong long; template <class Tint> T rd() {T res0;T fg1;char chgetcha…

牛客题霸 [ 排序] C++题解/答案

题目描述 给定一个数组&#xff0c;请你编写一个函数&#xff0c;返回该数组排序后的形式。 题意&#xff1a; 排序的方式有很多 二分呀&#xff0c;桶排呀归并等等 stl里的sort快速排序方便好用 题解&#xff1a; class Solution { public:/*** 代码中的类名、方法名、参…

ARC106E-Medals【hall定理,高维前缀和】

正题 题目链接:https://atcoder.jp/contests/arc106/tasks/arc106_e 题目大意 nnn个员工&#xff0c;第iii个在[1,Ai][1,A_i][1,Ai​]工作&#xff0c;[Ai1,2Ai][A_i1,2\times A_{i}][Ai​1,2Ai​]休息&#xff0c;[2Ai1,3Ai][2\times A_i1,3\times A_i][2Ai​1,3Ai​]工作…以…