博客
关于我
hdu6434 Problem I. Count(数论)(好题)
阅读量:395 次
发布时间:2019-03-05

本文共 1325 字,大约阅读时间需要 4 分钟。


Problem Description
Multiple query, for each n, you need to get
n i-1
∑ ∑ [gcd(i + j, i - j) = 1]
i=1 j=1

Input
On the first line, there is a positive integer T, which describe the number of queries. Next there are T lines, each line give a positive integer n, as mentioned above.
T<=1e5, n<=2e7

Output
Your output should include T lines, for each line, output the answer for the corre- sponding n.

Sample Input
4
978
438
233
666

Sample Output
194041
38951
11065
89963

Source
2018 Multi-University Training Contest 10
思路:在这里插入图片描述
变换差不多就是这样,于是对于每一个i就变成了求2i的欧拉函数除2了,求一下前缀和就行。

#include<bits/stdc++.h>using namespace std;typedef long long ll;const int N=2e7+5;ll euler[N<<1],sum[N],prim[N<<1];void get_euler(int n){    //欧拉函数    memset(euler, 0, sizeof euler);  euler[1] = 1;  int id = 0;  for(int i = 2; i < n; ++i) {       if(!euler[i]) {         euler[i] = i - 1;      prim[id++] = i;    }    for(int j = 0; j < id && prim[j]*i <n; ++j) {         if(i % prim[j]) {           euler[i*prim[j]] = euler[i] * (prim[j]-1);      } else {           euler[i*prim[j]] = euler[i] * prim[j];        break;      }    }  }}int main(){       int n,T;    get_euler(N<<1);    sum[1]=0;    for(int i=2;i<=N;++i) sum[i]=euler[2*i]/2,sum[i]+=sum[i-1];    scanf("%d",&T);    while(T--)    {           scanf("%d",&n);        printf("%lld\n",sum[n]);    }}

转载地址:http://sqewz.baihongyu.com/

你可能感兴趣的文章
pipeline配置与管理Job企业级实战
查看>>
pipeline项目配置实战
查看>>
Pipenv 与 Conda?
查看>>
QVGA/HVGA/WVGA/FWVGA分辨率屏含义及大小//Android虚拟机分辨率
查看>>
pipreqs : 无法将“pipreqs”项识别为 cmdlet、函数、脚本文件或可运行程序的名称。请检查名称的拼写,如果包括路径,请确保路径 正确,然后再试一次。
查看>>
pipy国内镜像的网址
查看>>
quiver绘制python语言
查看>>
pip下载缓慢
查看>>
PIP使用SSH从BitBucket安装自定义软件包,无需输入SSH密码
查看>>
pip命令提示unknow or unsupported command install解决方法
查看>>
pip在安装模块时提示Read timed out
查看>>
pip更换源
查看>>
SpringBoot之Banner源码深度分解
查看>>
Pix2Pix如何工作?
查看>>
QuickBI助你成为分析师——搞定数据源
查看>>
pkl来存储python字典
查看>>
quick sort | 快速排序 C++ 实现
查看>>
pkpmbs 建设工程质量监督系统 Ajax_operaFile.aspx 文件读取漏洞复现
查看>>
pkpmbs 建设工程质量监督系统 文件上传漏洞复现
查看>>
pku 2400 Supervisor, Supervisee KM求最小权匹配+DFS回溯解集
查看>>