题目描述 Description###

\(φ\) 函数是数论中非常常用的函数。对于正整数 \(x\) ,\(φ(x)\) 表示不超过 \(x\) 的所有正整数与 \(x\) 互质的个数。

现在我们对它进行一次拓展:对于正整数 \(x\) ,\(y\) ,定义 \(φ(x,y)\) 表示不超过 \(y\) 的所有正整数与 \(x\) 互质的个数。

现在我们给定正整数 \(n\) 和 \(m\) ,对于所有不超过 \(n\) 的正整数 \(i\) ,求 \(φ(i,m)\) 。

输入描述 Input Description###

输入仅一行两个正整数 \(n\) 和 \(m\) 。

输出描述 Output Description###

输出 \(n\) 行,每行一个整数。第 \(i\) 行表示 \(φ(i,m)\) 。

样例输入 Sample Input###

11 10

样例输出 Sample Output###

10
5
7
5
8
3
9
5
7
4
10

数据范围及提示 Data Size & Hint###

\(n \leq 10^5,m \leq 10^{15}\)

之前的一些废话###

题解###

正难则反,我们考虑\(m\) 以内能被\(n\) 整除的所有数,含有所有\(n\) 的因子的很显然是可以的,然后我们就可以通过容斥原理来加加减减来求得这个值,找找规律可以发现,含有平方因子的数不会参与其中的讨论,含有奇数个质因子的数应该是加,含有偶数个质因子的数应该是减,所以我们来一个素数筛来预处理每一个数的符号,然后$\sqrt{i} $ 的时间枚举每一个i的因数,然后进行计算。复杂度$ O(预处理+n\sqrt{n} )$

代码###

#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
typedef long long LL;
typedef pair<int,int> PII;
#define X first
#define Y second
#define mp make_pair
#define mem(a,b) memset(a,b,sizeof(a))
inline LL read()
{
LL x=0,f=1;char c=getchar();
while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
while(isdigit(c)){x=x*10+c-'0';c=getchar();}
return x*f;
}
const int maxn=100000;
int n,prime[maxn],end,sign[maxn+20];
LL m,sum;
bool vis[maxn+20];
void init()
{
for(int i=2;i<=maxn;i++)
if(!vis[i])
{
prime[++end]=i;
for(LL j=(LL)i*(LL)i;j<=(LL)maxn;j+=i)vis[j]=1;
}
for(int i=2;i<=maxn;i++)
{
if(!vis[i])sign[i]=1;
if(!sign[i])continue;
for(int j=1;j<=end && (LL)i*(LL)prime[j]<=maxn;j++)
{
if(i%prime[j]==0)sign[i*prime[j]]=0;
else sign[i*prime[j]]=-sign[i];
}
}
}
int main()
{
init();
n=(int)read();m=read();
printf("%lld\n",m);
for(int i=2;i<=n;i++)
{
int s=(int)sqrt(i);
sum=0ll;
for(int j=1;j<=s;j++)
if(i%j==0)
{
if(j!=1)sum+=(m/j)*sign[j];
if(j*j!=i)sum+=(m/(i/j))*sign[i/j];
}
printf("%lld\n",m-sum);
}
return 0;
}

总结###

找找规律画画图这种题还是比较简单的。

事后发现这玩意居然是莫比乌斯函数,我居然自己推出了这么高深的东西!

NOIP模拟赛 拓展的更多相关文章

  1. NOIP模拟赛20161022

    NOIP模拟赛2016-10-22 题目名 东风谷早苗 西行寺幽幽子 琪露诺 上白泽慧音 源文件 robot.cpp/c/pas spring.cpp/c/pas iceroad.cpp/c/pas ...

  2. contesthunter暑假NOIP模拟赛第一场题解

    contesthunter暑假NOIP模拟赛#1题解: 第一题:杯具大派送 水题.枚举A,B的公约数即可. #include <algorithm> #include <cmath& ...

  3. NOIP模拟赛 by hzwer

    2015年10月04日NOIP模拟赛 by hzwer    (这是小奇=> 小奇挖矿2(mining) [题目背景] 小奇飞船的钻头开启了无限耐久+精准采集模式!这次它要将原矿运到泛光之源的矿 ...

  4. 大家AK杯 灰天飞雁NOIP模拟赛题解/数据/标程

    数据 http://files.cnblogs.com/htfy/data.zip 简要题解 桌球碰撞 纯模拟,注意一开始就在袋口和v=0的情况.v和坐标可以是小数.为保险起见最好用extended/ ...

  5. 队爷的讲学计划 CH Round #59 - OrzCC杯NOIP模拟赛day1

    题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的讲学计划 题解:刚开始理解题意理解了好半天,然后发 ...

  6. 队爷的Au Plan CH Round #59 - OrzCC杯NOIP模拟赛day1

    题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的Au%20Plan 题解:看了题之后觉得肯定是DP ...

  7. 队爷的新书 CH Round #59 - OrzCC杯NOIP模拟赛day1

    题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的新书 题解:看到这题就想到了 poetize 的封 ...

  8. CH Round #58 - OrzCC杯noip模拟赛day2

    A:颜色问题 题目:http://ch.ezoj.tk/contest/CH%20Round%20%2358%20-%20OrzCC杯noip模拟赛day2/颜色问题 题解:算一下每个仆人到它的目的地 ...

  9. CH Round #52 - Thinking Bear #1 (NOIP模拟赛)

    A.拆地毯 题目:http://www.contesthunter.org/contest/CH%20Round%20%2352%20-%20Thinking%20Bear%20%231%20(NOI ...

随机推荐

  1. git 代码管理

  2. 基于mysqld_multi实现MySQL 5.7.24多实例多进程配置

    学习环境: 操作系统 IP地址 主机名 软件包 备注 CentOS7.5 192.168.200.111 localhost       实验初始配置:所有主机关闭防火墙与selinux [root@ ...

  3. HDU1075 What Are You Talking About(map)

    传送门 题目大意:一个单词对应另一个单词 翻译一段文字 题解:stl map走一波 代码: #include<iostream> #include<map> #include& ...

  4. PHP 7.4.0 发布

    近日,PHP 7.4.0 发布了,此版本标志着 PHP 7 系列的第四次特性更新. PHP 7.4.0 进行了许多改进,并带来了一些新特性,包括: Typed Properties  类型属性 类属性 ...

  5. 资源推荐:PPT快闪资源合集附配套字体下载

    样例ppt下载 搜索公众号“拒收”或扫码关注公众号 回复关键字“快闪ppt”获取全部福利 本公众号只出精品,拒收劣质 或者点击菜单链接获取获取全部资源

  6. 英语阅读——Speaking Chinese in America

    这篇文章是<新视野大学英语>第四册的第五单元的文章,第一遍英语阅读完后对比中文,发现自己对作者的观点理解有些出入.作者反对的是认为中国说话客套而美国人直接的观点,利用自己的经历表达了中文也 ...

  7. jquery 图片缩放插件使用

    https://cdn.bootcss.com/d3-transition/1.1.1/d3-transition.js http://www.jq22.com/jquery-info9291

  8. ActiveMQ笔记之ConnectionFactory

    一.ActiveMQ原生的连接工程:ActiveMQConnectionFactory 默认的maxThreadPoolSize=1000,也就是每个connection的session线程池最大值为 ...

  9. MySQL(9)---纪录一次实际开发过程中用到的复杂存储过程

    Mysql(9)---纪录一次实际开发过程中用到的复杂存储过程 为了尽可能的还原当时为什么需要用到存储过程,下面我写了个详细的文档,我们可以从需求文档出发来分析. 有关存储过程之前也写了两篇文章来做铺 ...

  10. C# System.Timers Demo

    static void Main(string[] args) { SystemTimersDemo(); Console.ReadLine(); } static void SystemTimers ...