【数论】挖掘机技术哪家强(jzoj 3858)

挖掘机技术哪家强

jzoj 3858

题目大意

给你多个n,对于每个n,x为n的因数,设s(x)为小于x且与x互质的数的和,让你求每一个x的s(x)总和

原题

有人问现实中为什么总是男生追求女生,反过来很少。实际上女生也是想主动追求男生的,但是世俗中对于主动追求男生的女生有种歧视,这样就使得女生不大敢主动追求男生。但是面对喜欢的男生,难道就不出手么?女生只能步步为营,挖坑来引诱男生往里跳。这时候问题就来了,挖掘机技术到底哪家强?
被热血沸腾的广告洗脑了若干天后,Matt终于下定决心,毅然登上了开往泉城的列车,决心寻找生活的希望。
来到布鲁谢特学院后,Matt逐渐地了解了各种型号的挖掘机。在这里我们可以认为有大挖掘机和小挖掘机两种。
今天Matt的任务很简单:首先他要用大挖掘机挖出恰好N单位体积的砂土。由于小挖掘机比较笨拙,它每次挖的砂土体积是固定的。也就是说,设每次挖x单位体积砂土,那么N需要被x整除。在挖出若干堆体积为x的砂土后,Matt需要计算x的“难挖指数”。体积x的“难挖指数”定义如下:对于某个不超过x的体积y,如果x与y的最大公约数为1,那我们认为体积y是“难挖的”,x的“难挖指数”就要加上y。
由于Matt之后还需要用小挖掘机处理被大挖掘机挖出的砂土,他希望知道所有可能的x的难挖指数的和,这样他好估算今天要干多久,不然作为布鲁谢特的高才生,他出门要被笑话的。

输入样例

3
2
3
4

输出样例

2
4
6

数据范围

对于30%的数据有T⩽20,N⩽104。T\leqslant 20,N\leqslant 10^4。T20N104
对于60%的数据有T⩽100,N⩽107。T\leqslant 100,N\leqslant 10^7。T100N107
对于100%的数据有1⩽T⩽1000,1⩽N⩽109。1\leqslant T\leqslant 1000,1\leqslant N\leqslant 10^9。1T10001N109

解题思路

挖掘技术肯定是蓝翔啦(滑稽)
我们可以用欧拉函数x/φ(x)∗2x/φ(x)*2x/φ(x)2来求s(x)s(x)s(x)的值,具体可以上百度学习一下
然后我们把n质因数分解,然后再dfs枚举每一个质数的指数,同时计算φ(x)φ(x)φ(x)的值就行了

代码:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
ll t, n, w, ans, prime[32000], s[32000], l[32000], p[32000];
void dfs(ll x, ll num, ll phi)
{if (x > w){ans += num * phi / 2;//计算总数就要这样if (num == 1) ans++;//1的情况return;}dfs(x + 1, num, phi);//指数为0ll sum = 1;for (ll i = 1; i <= l[x]; ++i){sum *= s[x];//总数dfs(x + 1, num * sum, phi * (sum / s[x]) * (s[x] - 1));//sum是计算n的一部分然后就是/s[x]*(s[x] - 1)}
}
int main()
{for (ll i = 2; i <= 32000; ++i)if (!p[i]){prime[++prime[0]] = i;//计算质因子for (ll j = i * i; j <= 32000; j += i)p[j] = 1;}scanf("%lld", &t);while(t--){ans = 0;w = 0;memset(l, 0, sizeof(l));scanf("%lld", &n);for (ll i = 1; i <= prime[0]; ++i)if (n % prime[i] == 0)//看看有没有这个质因子{s[++w] = prime[i];while(n % prime[i] == 0){l[w]++;n /= prime[i];}}if (n > 1) s[++w] = n, l[w] = 1;//大于sqrt(n)的最多有一个dfs(1, 1, 1);printf("%lld\n", ans);}return 0;
}

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

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

相关文章

.net core grpc consul 实现服务注册 服务发现 负载均衡(二)

在上一篇 .net core grpc 实现通信(一) 中&#xff0c;我们实现的grpc通信在.net core中的可行性&#xff0c;但要在微服务中真正使用&#xff0c;还缺少 服务注册&#xff0c;服务发现及负载均衡等&#xff0c;本篇我们将在 .net core grpc 通信 的基础上加上 服务注册&#x…

【BZOJ3218】a+b problem (最小割 + 主席树)

传送门 继续优化&#xff1a;把a[ ]离散化 #include<bits/stdc.h> using namespace std; const int inf1000000007; const int N200010; const int M1000010;struct Edge{int u,v,f,next; }edge[M]; int head[N],cnt; int s,t,flow,level[N];struct Node{int x,id; }e[…

ST表讲解

ST表主要用于解决RMQ问题&#xff08;区间最值问题&#xff09; 当然你可以用线段树等&#xff0c;但今天用一种ST表&#xff08;倍增算法&#xff09; ST表是倍增算法的一个典型应用 暴力做RMQ问题&#xff0c;往往会超时&#xff0c;ST表利用对其进行优化 给定一段序列A&am…

P3346-[ZJOI2015]诸神眷顾的幻想乡【广义SAM】

正题 题目链接:https://www.luogu.com.cn/problem/P3346 题目大意 一棵树&#xff0c;求树上所有路径构成的字符串有多少种。&#xff08;叶子不超过303030个&#xff09; 解题思路 如果是根节点到一些节点的路径的话很好做&#xff0c;直接建广义SAMSAMSAM即可&#xff0c;…

纪中B组模拟赛总结(2020.2.1)

成绩 rankrankranknamenamenamescorescorescoreT1T1T1T2T2T2T3T3T3T4T4T4414141lyflyflyf909090000909090000000 总结 T1T1T1直接爆炸&#xff0c;没想到是数学问题 T2T2T2想了出来&#xff0c;但忘记开mapmapmap了 T3T3T3直接懵逼 T4T4T4没有太了解欧拉函数导致翻车 题解 T1…

微软技术直通车(第三期)现场实录

微软技术直通车&#xff08;第三期&#xff09;于本月18日14时&#xff0c;在中国微软总部圆满举行。本次活动得到了微软最佳有价值专家&#xff08;Microsoft MVP&#xff09;中国区项目组的鼎力支持。同时&#xff0c;Microsoft MVP中国区项目组负责人也亲临现场致辞&#xf…

[NOI2019] 序列(模拟费用流)

原先自己想的建图&#xff1a; 正确建图&#xff1a; 但是 n 太大了&#xff0c;所以考虑模拟费用流&#xff1a; 注意&#xff1a; 在1中&#xff0c; 若选的两个位置相同&#xff0c;则为情况2&#xff0c;不用减 f&#xff1b; 若选的位置在另一序列中已被选&#xff0c;…

【结论】游戏(jzoj 5536)

游戏 jzoj 5536 题目大意&#xff1a; 给一个范围和一个矩形&#xff0c;让你在这个范围内放若干个这样的矩形&#xff08;不能改变方向&#xff0c;不能重叠&#xff09;&#xff0c;让你求出最少放多个矩形可以使范围内无法再放矩形 输入样例#1 11 4 3 2输入样例#2 10 …

[2020.11.26NOIP模拟赛]询问【字符串hash】

正题 题目链接:https://www.luogu.com.cn/problem/U142342?contestId37784 题目大意 一个字符串&#xff0c;定义两个字符串相似为用一些字母代替相同的字母后可以相同。 如urbbrgroorurbbrgroorurbbrgroor&#xff0c;apple≠abcdeapple\neq abcdeapple​abcde 要求支持…

LCA总结

文章目录LCA介绍解决方法概括&#xff1a;倍增法&#xff1a;Tarjan欧拉序树剖解法&#xff1a;看了很多关于LCA的文章&#xff0c;这里是一个总结我们通过这个题洛谷P3379 【模板】最近公共祖先来讲LCA LCA介绍 lca是啥&#xff1f;最近公共祖先 就是&#xff1a;两个点在这…

ASP.NET Core Web API + Ng6 实战视频 Day 2

第一天课程&#xff1a;ASP.NET Core Web API Identity Server 4 Angular 6 实战小项目视频Day 2 第一部分:Day 2 第二部分:视频专辑持续更新中....地址请点击原文链接.原文地址: http://v.qq.com/vplus/4cfb00af75c16eb8d198c58fb86eb4dc/foldervideos/8hk0029019k2fft.NET…

浮沉子

浮沉子制作及其原理 浮沉子的下沉过程&#xff1a;

纪中B组模拟赛总结(2020.2.3)

成绩 rankrankranknamenamenamescorescorescoreT1T1T1T2T2T2T3T3T3333lyflyflyf252.6252.6252.610010010010010010052.652.652.6 总结 T1以前做过&#xff0c;十分嗨皮,就直接切掉了 T2发现是一道傻逼题,不多说 T3比赛时理解题意错了&#xff0c;以为一个字符串属于一个组别就…

[2020.11.26NOIP模拟赛]勇者的后缀【SA,RMQ,主席树,二分】

正题 题目链接:https://www.luogu.com.cn/problem/U142356?contestId37784 题目大意 一个字符串&#xff0c;询问给出(x,l,r)(x,l,r)(x,l,r)表示询问在[l,r][l,r][l,r]中作为起点找一个后缀它与xxx作为起点的后缀的LCPLCPLCP最长&#xff0c;且满足最长的情况下字典序最小。…

.NET Core 项目指定SDK版本

一. 版本里的坑自从 .NET Core 2.1.0版本发布以后&#xff0c;近几个月微软又进行了几次小版本的发布&#xff0c;可见 .NET Core 是一门生命力非常活跃的技术。经过一段时间的实践&#xff0c;目前做 ASP.NET Core 开发时&#xff0c;使用的 Nuget 包&#xff0c;比如 Microso…

牛客网【每日一题】4月17日题目精讲 华华给月月准备礼物

文章目录题目描述题解&#xff1a;代码&#xff1a;推荐例题&#xff1a;试题链接时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 32768K&#xff0c;其他语言65536K 64bit IO Format: %lld 题目描述 二月中旬虐狗节前夕&#xff0c;华华决定给…

[BZOJ 3173] [TJOI 2013] 最长上升子序列(splay)

根据题意很容易得到最原始的做法&#xff1a; #include<iostream> #include<cstdio> #include<cstring> #include<vector> using namespace std; const int N100005; int n,m,f[N]; int main(){scanf("%d",&n);int x;for(int i1;i<n…

P4300-[AHOI2006]上学路线【网络流,最短路】

正题 题目链接:https://www.luogu.com.cn/problem/P4300 题目大意 nnn个点mmm条边的无向图。求1∼n1\sim n1∼n的最短路和删除cic_ici​和最小的边使得最短路变长。 解题思路 显然我们需要跑一次最短路。 之后考虑如何求第二问&#xff0c;我们发现我们要割掉最短路上的边&a…

Mathematics(jzoj 1747)

Mathematics jzoj 1747 题目大意 有n堆石子&#xff0c;总和为2k2^k2k&#xff0c;现在对于两堆石子&#xff0c;你可以从a中取b的分量到b&#xff08;a要大于b&#xff09;&#xff0c;问合成一堆大小为2k2^k2k的石子要怎么做&#xff08;输出a,b&#xff09; 输入样例 2…

Bing.com在.NET Core 2.1上运行

相关知识请参考.netCore开发团队博客 https://blogs.msdn.microsoft.com/dotnet/ Bing.com是一种云服务&#xff0c;运行在遍布全球许多数据中心的数千台服务器上。Bing服务器每秒处理来自全球消费者的数千个用户查询&#xff0c;通过他们的浏览器&#xff0c;使用Microsoft认…