Written with StackEdit.

Description

\(Z\)城市居住着很多只跳蚤。在\(Z\)城市周六生活频道有一个娱乐节目。一只跳蚤将被请上一个高空钢丝的正中央。钢丝很长,可以看作是无限长。节目主持人会给该跳蚤发一张卡片。卡片上写有\(N+1\)个自然数。其中最后一个是\(M\),而前\(N\)个数都不超过\(M\),卡片上允许有相同的数字。跳蚤每次可以从卡片上任意选择一个自然数\(S\),然后向左,或向右跳\(S\)个单位长度。而他最终的任务是跳到距离他左边一个单位长度的地方,并捡起位于那里的礼物。比如当\(N=2,M=18\)时,持有卡片\((10, 15, 18)\)的跳蚤,就可以完成任务:他可以先向左跳\(10\)个单位长度,然后再连向左跳\(3\)次,每次\(15\)个单位长度,最后再向右连跳\(3\)次,每次\(18\)个单位长度。而持有卡片\((12, 15, 18)\)的跳蚤,则怎么也不可能跳到距他左边一个单位长度的地方。当确定\(N\)和\(M\)后,显然一共有\(M^N\)张不同的卡片。现在的问题是,在这所有的卡片中,有多少张可以完成任务。

Input

输入文件有且仅有一行,包括用空格分开的两个整数\(N\)和\(M\)。

Output

输出文件有且仅有一行,即可以完成任务的卡片数.\(1≤M≤10^8,1≤N≤M,\)且\(M^N≤10^16.\)

Sample Input

2 4

Sample Output

12

Solution

  • 显然,根据裴蜀定理,满足\(gcd(a_1,a_2,a_3,......a_n,M)=1\)的即为一组解.我们要对其进行计数.
  • 考虑用所有的方案\(M^N\)减去\(gcd\)不为\(1\)的数.
  • 因为最后一个数为\(M\),所以\(gcd\)一定是\(M\)的约数.
  • 将\(M\)分解质因数后,二进制枚举质因数的组合情况.
  • 若当前组合出来的数是\(prod\),那么是\(prod\)倍数的有\(\frac{p}{prod}\)个.
  • 那么这样的方案数为\((\frac{p}{prod})^N.\)
  • 由于会产生重复计数,所以需要容斥处理.奇数个质因子\(+\),偶数个\(-\).

这样看来,\(M^N\)也是一种特殊情况,\(prod=1\),个数为\(0\),所以用加.

#include<bits/stdc++.h>
using namespace std;
typedef long long LoveLive;
inline int read()
{
int out=0,fh=1;
char jp=getchar();
while ((jp>'9'||jp<'0')&&jp!='-')
jp=getchar();
if (jp=='-')
{
fh=-1;
jp=getchar();
}
while (jp>='0'&&jp<='9')
{
out=out*10+jp-'0';
jp=getchar();
}
return out*fh;
}
int n,m;
LoveLive fpow(LoveLive a,LoveLive b)
{
LoveLive res=1;
while(b)
{
if(b&1)
res*=a;
a*=a;
b>>=1;
}
return res;
}
LoveLive factor[50];
int main()
{
n=read(),m=read();
LoveLive ans=fpow(m,n);
int p=m,siz=0;
for(int i=2;i*i<=p;++i)
if(p%i==0)
{
factor[++siz]=i;
while(p%i==0)
p/=i;
}
if(p>1)
factor[++siz]=p;
LoveLive lim=(1<<siz)-1;//质因数的组合方式,二进制枚举
for(LoveLive i=1;i<=lim;++i)
{
LoveLive prod=1,cnt=0;//乘积与这种方案的质因数个数
for(int j=1;j<=siz;++j)
if((i>>(j-1))&1)
{
prod*=factor[j];
++cnt;
}
if(cnt&1)
ans-=fpow(m/prod,n);
else
ans+=fpow(m/prod,n);
}
printf("%lld\n",ans);
return 0;
}

bzoj 1220 跳蚤的更多相关文章

  1. bzoj 4310: 跳蚤

    Description 很久很久以前,森林里住着一群跳蚤.一天,跳蚤国王得到了一个神秘的字符串,它想进行研究. 首先,他会把串分成不超过 k 个子串,然后对于每个子串 S,他会从S的所有子串中选择字典 ...

  2. ●BZOJ 4310 跳蚤

    ●赘述题目 给出一个字符串,要求分成k个子串,然后求出每个子串的字典序最大的子串(我称它为子子串),要使这k个子子串中的字典序最大的那个串(即魔力串)最小.输出该魔力串. (本题个人感觉很好,比较综合 ...

  3. bzoj 4310 跳蚤——后缀数组+二分答案+贪心

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4310 答案有单调性? 二分出来一个子串,判断的时候需要满足那些字典序比它大的子串都不出现! ...

  4. bzoj 4310 跳蚤 二分答案+后缀数组/后缀树

    题目大意 给定\(k\)和长度\(\le10^5\)的串S 把串分成不超过\(k\)个子串,然后对于每个子串\(s\),他会从\(s\)的所有子串中选择字典序最大的那一个,并在选出来的\(k\)个子串 ...

  5. bzoj 4310 跳蚤 —— 后缀数组+二分答案+贪心

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4310 二分答案——在本质不同的子串中二分答案! 如果二分到的子串位置是 st,考虑何时必须分 ...

  6. bzoj 4310: 跳蚤【后缀数组+st表+二分+贪心】

    先求一下SA 本质不同的子串个数是\( \sum n-sa[i]+1-he[i] \),按字典序二分子串,判断的时候贪心,也就是从后往前扫字符串,如果当前子串串字典序大于二分的mid子串就切一下,然后 ...

  7. 后缀数组 hash求LCP BZOJ 4310: 跳蚤

    后缀数组的题博客里没放进去过..所以挖了一题写写 充实下博客 顺便留作板子.. 一个字符串S中 内容不同的子串 有 sigma{n-sa[i]+1-h[i]}   (噢 这里的h[]就是大家熟知的he ...

  8. 跳蚤 BZOJ 4310

    跳蚤 [问题描述] 很久很久以前,森林里住着一群跳蚤.一天,跳蚤国王得到了一个神秘的字符串,它想进行研究. 首先,他会把串分成不超过 k 个子串,然后对于每个子串 S,他会从S的所有子串中选择字典序最 ...

  9. bzoj千题计划157:bzoj1220:[HNOI2002]跳蚤

    扩展欧几里得:ax+by=gcd(a,b) 一定有解 能跳到左边一格,即ax+by=-1 若a,b的gcd=1,则一定有解 所以问题转化为 求n个不大于m的数,他们与m的gcd=1 的方案数 容斥原理 ...

随机推荐

  1. Python面试题之Python正则表达式指南

    1. 正则表达式基础 1.1. 简单介绍 正则表达式并不是Python的一部分.正则表达式是用于处理字符串的强大工具,拥有自己独特的语法以及一个独立的处理引擎,效率上可能不如str自带的方法,但功能十 ...

  2. clipbrd剪切板查看器

    本文,我们来学习一下简单的概念,即,如何查看系统剪贴版里面有什么内容?   如果要想看.或者验证系统剪贴版里面都有什么内容,最为简单的方法就是通过"粘贴"的操作来验证!   但是, ...

  3. 第十节课-RNN介绍

    2017-08-21 这次的课程介绍了RNN的相关知识: 首先是RNN的几种模型: 分别又不同的应用场景,包括机器翻译,视频的分类... RNN的解释: 主要的特点就是用到了上一个隐含状态的信息,所以 ...

  4. BZOJ:(270,300]

    9/30 BZOJ3038:线段树,不带lazy标记,直接修改叶子. BZOJ3211:同3038 BZOJ1406:将式子转换成[(x-1)*(x+1)%n==0]然后枚举i.当i=x-1时,i*( ...

  5. java基础(3)--8种基本类型

    八种基本类型 数据类型分为两大类:基本类型,引用类型 基本类型有8种,除了基本类型(8种)以外任何类型都是引用类型.如: String 是引用类型 基本类型有8种(首字母是小写): 整数:byte, ...

  6. Hibernate的一级缓存、二级缓存和查询缓存。

    Hibernate的Session提供了一级缓存的功能,默认总是有效的,当应用程序保存持久化实体.修改持久化实体时,Session并不会立即把这种改变提交到数据库,而是缓存在当前的Session中,除 ...

  7. Android -- ContentProvider 内容提供者,创建和调用

    1. 概述 ContentProvider 在android中的作用是对外共享数据,也就是说你可以通过ContentProvider把应用中的数据共享给其他应用访问,其他应用可以通过ContentPr ...

  8. 十二道MR习题 - 2 - 多文件保存

    题目: 需要将MR的执行结果保存到3个文件中,该怎么做. 又是一个送分题. 对于Hadoop的MapReduce来说只需要设置一下reduce任务的数量即可.MR的Job默认reduce数量是1,需要 ...

  9. selenium学习笔记(智能等待)

    博主在尝试对百度首页用selenium完成自动登录的功能 反复多次尝试元素定位方法也未写错.最后发现问题原因: 脚本运行速度快于页面加载速度 如百度首页登录例子.脚本已经开始寻找登录弹窗 但是页面仍在 ...

  10. JavaEE 技术体系

    JavaEE 技术体系总结: 一:常见模式与工具 设计模式,流行的框架与组件 常见的设计模式,编码必备 Spring5,做应用必不可少的最新框架 MyBatis,玩数据库必不可少的组件 二:工程化与工 ...