博客
关于我
[bzoj2818][莫比乌斯反演]Gcd
阅读量:91 次
发布时间:2019-02-26

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

Description

给定整数N,求1<=x,y<=N且Gcd(x,y)为素数的 数对(x,y)有多少对.

Input

一个整数N

Output

如题

Sample Input

4

Sample Output

4

HINT

对于样例(2,2),(2,4),(3,3),(4,2)

1<=N<=10^7

题解

重学莫反

这里写图片描述

#include
#include
#include
#include
#include
using namespace std;typedef long long LL;int mu[11000000],prime[11000000],pr,n;bool v[11000000];void getmu(){ memset(v,true,sizeof(v)); mu[1]=1;pr=0; for(int i=2;i<=10000000;i++) { if(v[i]==true) { prime[++pr]=i; mu[i]=-1; } for(int j=1;i*prime[j]<=10000000 && j<=pr;j++) { v[i*prime[j]]=false; if(i%prime[j]==0) { mu[i*prime[j]]=0; break; } else mu[i*prime[j]]=-mu[i]; } }}int main(){ scanf("%d",&n); getmu();LL ans=0; for(int i=1;prime[i]<=n;i++) { for(LL j=1;j<=(LL)n/prime[i];j++) ans+=(LL)mu[j]*((n/prime[i])/j)*((n/prime[i])/j); } printf("%lld\n",ans); return 0;}
你可能感兴趣的文章
MTD技术介绍
查看>>
MySQL
查看>>
MySQL
查看>>
mysql
查看>>
MTK Android 如何获取系统权限
查看>>
MySQL - 4种基本索引、聚簇索引和非聚索引、索引失效情况、SQL 优化
查看>>
MySQL - ERROR 1406
查看>>
mysql - 视图
查看>>
MySQL - 解读MySQL事务与锁机制
查看>>
MTTR、MTBF、MTTF的大白话理解
查看>>
mt_rand
查看>>
mysql -存储过程
查看>>
mysql /*! 50100 ... */ 条件编译
查看>>
mudbox卸载/完美解决安装失败/如何彻底卸载清除干净mudbox各种残留注册表和文件的方法...
查看>>
mysql 1264_关于mysql 出现 1264 Out of range value for column 错误的解决办法
查看>>
mysql 1593_Linux高可用(HA)之MySQL主从复制中出现1593错误码的低级错误
查看>>
mysql 5.6 修改端口_mysql5.6.24怎么修改端口号
查看>>
MySQL 8.0 恢复孤立文件每表ibd文件
查看>>
MySQL 8.0开始Group by不再排序
查看>>
mysql ansi nulls_SET ANSI_NULLS ON SET QUOTED_IDENTIFIER ON 什么意思
查看>>