没做完的题

news/2025/10/10 21:29:32/文章来源:https://www.cnblogs.com/Vinson-/p/19133911

放在这里待办
定义

\[g(n)=2\sum ^n_{i=2}\left \lfloor \frac{n}{i} \right \rfloor^2 \phi(i)+n^2 \]

\(g(n)\)

定义

\[f(n)=\frac{g(n)-n^2}{2} \]

\[f(n)=\sum ^n_{i=2}\left \lfloor \frac{n}{i} \right \rfloor^2 \phi(i) \]

\[f(n-1)=\sum ^{n-1}_{i=2}\left \lfloor \frac{n-1}{i} \right \rfloor^2 \phi(i) \]

\[f(n)-f(n-1)=\sum ^{}_{i|n}((\frac {n}{i})^2-(\frac {n}{i}-1)^2)\phi(i) \]

\[=2\sum ^{}_{i|n}\frac {n}{i}\phi(i) - n \]

定义

\[g(n)=2f(n)+n^2 \]

\[sum_g=2sum_f+\frac{n(n+2)(2n+1)}{6} \]

#include<bits/stdc++.h>//记得开int128
#define int long long
#define ll __int128
using namespace std;
const int N=1e7+5;
int p[N],cnt,phi[N],s[N],low[N],T,sum[N],sums[N];
bool isp[N];
unordered_map<int,int>ans;
void Euler(int n){s[1]=low[1]=phi[1]=isp[1]=sum[1]=1;for(int i=2;i<=n;i++){if(!isp[i])p[++cnt]=i,phi[i]=i-1,s[i]=2*i-1,low[i]=i;sum[i]=sum[i-1]+phi[i],sums[i]=sums[i-1]+2*s[i]-i;for(int j=1;j<=cnt&&i*p[j]<=n;j++){isp[i*p[j]]=1;if(i%p[j]==0){low[i*p[j]]=low[i]*p[j],phi[i*p[j]]=phi[i]*p[j];if(low[i]==i){s[i*p[j]]=s[i]*p[j]+phi[i*p[j]];}else{s[i*p[j]]=s[i/low[i]]*s[low[i]*p[j]];}break;}low[i*p[j]]=p[j],s[i*p[j]]=s[i]*s[p[j]],phi[i*p[j]]=phi[i]*phi[p[j]];}}for (int i=1;i<=10;i++) cout<<s[i]<<" "; cout<<"\n";
}
int djs(int n){if(n<=N-5)return sum[n];if(ans[n])return ans[n];__int128 res=(__int128)n*(n+1)/2;for(int l=2,r;l<=n;l=r+1){r=n/(n/l);res-=(r-l+1)*djs(n/l);}return ans[n]=res;
}
signed main(){freopen("1.in","r",stdin);// freopen("1.out","w",stdout);cin>>T;Euler(N-5);while(T--){int n;cin>>n;if(n<=N-5){cout<<2*s[n]+n*n<<endl;continue;}int res=0;for(int l=1,r;l<=n;l=r+1){r=n/(n/l);res+=(djs(r)-djs(l-1))*(n/l)*(n/l);}cout<<2*res-n*n<<endl;}
}

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

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

相关文章

JavaScriptDay1

1.什么是JavaScript 1.1 概述 JavaScript是一门世界上最流行的脚本语言 Java、JavaScript 一个合格的后端人员,必须要精通JavaScript 2.快速入门 2.1、引入javaScript内部标签 <script>//... </script> 外…

淘宝NPM镜像地址https://registry.npm.taobao.org不可用

1 地址已经更新了,https://registry.npm.taobao.org已经停服2 新地址是 https://registry.npmmirror.com

星星充电一面

1、你是否了解ThreadLocal吗 2、如果有父线程和子线程,那么子线程是否能访问到父线程ThreadLocal里面的内容吗 3、那你的项目中只是一个单体项目吗,没有子线程?那如果我想要子线程访问父线程里面的内容怎么办 4、布…

6 CF1034 div3 题解

CF1034 div3 打了一下虚拟赛,感觉自己写代码不够认真,前面三个题都是能秒的,结果还寄了几发 这波直接写不会的题的题解吧,节省时间 F 题面 给定一个 \(n\) ,对于一个长度为 \(n\) 的排列,称满足下面条件的排列为…

5 ABC413 题解

ABC 413 题解 D 题面 给定一个长度为 \(N\) 的整数序列 \(A=(A_1,A_2,\ldots,A_N)\) 。 保证对于任何 \(i\ (1\le i\le N)\) , \(A_i\) 不为 \(0\) 。 判断是否存在 \(A\) 的一个排列 \(B=(B_1,B_2,\ldots,B_N)\) 使得…

4 CF 1032 div3 题解

31 CF 1032 div3 题解 C 题面 给你一个行数为 \(n\) 列数为 \(m\) 的整数矩阵。在第 \(i\) 行和第 \(j\) 列的交叉处的单元格中包含数字 \(a_{ij}\) 。 您可以执行以下操作次:选择两个数字 \(1 \leq r \leq n\) 和 \(…

3 ABC411 C ~ E题解

ABC411 C~E 题解 又是赤石的一天 C 有个长度为 \(N\) 的序列,初始序列中每个数为0 每次操作给定 \(pos\) ,将 pos^1 ,然后输出序列中有多少段不连续的 1 用小样例模拟一下可得,设当前颜色为 \(b\) ,左边颜色为 \(…

9 ABC408 D~F 题解

ABC408 D~F 题解 D 题面 给定一个长度为 \(n\) 的由 01 组成的字符串 \(S\),每次操作可以将某个 0 改成 1,或者将某个 1 改成 0 。 求字符串中至多有一个连续 1 串的最小操作次数。 题解 解法1 考场思路,将每个连续…

8 ABC425 G 题解

ABC425 G 题面 给定两个正整数 \(N,M\) 以及一个长度为 \(N\) 的非负整数序列 \(A = A_1, A_2, ..., A_N\)。求 \[\sum_{x = 0}^{M - 1} \min_{1 \le i \le N} (x \oplus A_i) \]\(1 \le N \le 2 \times 10^5\) \(1 \l…

智能防御,安全赋能:AI-FOCUS 滤海AI DLP 化解外部 AI 风险

在“员工把数据投喂给第三方AI”的常见场景下,企业需要一款可在输入与上传瞬时介入的产品。AI-FOCUS 滤海AI DLP 以流式网关为核心,在用户向大模型输入文本与提交文件/图片时进行敏感数据识别,按风险等级执行放行、…

VS code 中代码补全 自动补全函数括号

使用vscode的pylance插件补全python代码的时候无法补全函数的括号 需要在settings.json中添加这一行即可

学习ReAct并使用langgraph实现一个简单的ReAct AI Agent!!

ReAct介绍 要介绍ReAct最好要知道它是从哪来的。 ReAct这个概念出自《REACT : SYNERGIZING REASONING AND ACTING IN LANGUAGE MODELS》这篇很牛的论文。 论文地址:https://arxiv.org/abs/2210.03629 我们先来看下这篇…

23种设计模式之【策略模式】-核心原理与 Java 实践 - 详解

pre { white-space: pre !important; word-wrap: normal !important; overflow-x: auto !important; display: block !important; font-family: "Consolas", "Monaco", "Courier New", …

abc 408 d~f

这次做了一次 abc,d 做出来了,但是比较麻烦,又用正确方法写了一遍,整理一下 d,e,f,g 有一些超纲。 abc408d 考虑把区间 \(l,r\) 最后变成 1,然后尝试去表示这个时候的答案。 \(sum[i]\) 表示 \(i\) 位置以及之前…

RMQ与LCA学习笔记

在开始之前先提一下RMQ与LCA这两个东西有什么关系 对于一个序列,对它构建出一颗笛卡尔树之后,两个点的LCA就是原序列中这两个点之间的最大值/最小值(取决于建树时的比较方式) 而对于一棵树,求出来他的欧拉序之后,…

the charm of Chinese language

The charm of Chinese language is you use it to read books. when I can read the original books translation as Chinese. I found the original one lets people dizzy, maybe it needs some graphs to make the …

mamba-硬件感知算法

扫描操作由于A B C这些矩阵现在是动态的了,因此无法使用卷积表示来计算它们(卷积核是固定的),因此,我们只能使用循环表示,如此也就而失去了卷积提供的并行训练能力 Mamba通过并行扫描(parallel scan)算法使得最终并…

完整教程:lua代码解析1

完整教程:lua代码解析1pre { white-space: pre !important; word-wrap: normal !important; overflow-x: auto !important; display: block !important; font-family: "Consolas", "Monaco", &qu…

system表空间丢失部分文件恢复---惜分飞

system表空间丢失部分文件恢复---惜分飞联系:手机/微信(+86 17813235971) QQ(107644445) 标题:system表空间丢失部分文件恢复 作者:惜分飞©版权所有[未经本人同意,不得以任何形式转载,否则有进一步追究法律责任…

二维数点

介绍 给出一个二维平面內的若干个点,多次询问某个矩形区域內包含多少个点(边界也算)。又或者,给一个长为 n nn 的序列,多次询问区间 [ l , r ] [l,r][l,r] 中值在 [ x , y ] [x,y][x,y] 内的元素个数。 例题 P190…