题面: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(N <= 15 , M <= 100000000)。
Output
可以完成任务的卡片数。
Sample
Input
2 4
Output
12
Hint
这12张卡片分别是:
(1, 1, 4), (1, 2, 4), (1, 3, 4), (1, 4, 4), (2, 1, 4), (2, 3, 4),
(3, 1, 4), (3, 2, 4), (3, 3, 4), (3, 4, 4), (4, 1, 4), (4, 3, 4)
题意:已知m,且a1 ~ an 均<=m, 求系数x1,x2, …, xn,x(n + 1), 时的不定方程:x1 * a1 + c2 * a2 + … + x(n + 1) * m = 1有解
思路:(看了很多博客,只是简单提了一下结论,在这里我帮忙梳理一下)
首先,我们得知道裴蜀定理:对于任意整数a, b, 存在一对整数x, y, 满足ax + by = gcd(a, b);
进而我们可以知道其推论:(n个整数间的裴蜀定理)
设a1,a2,a3…an为n个整数,d是它们的最大公约数,那么存在整数x1…xn使得x1a1+x2a2+…xn*an=d
证明过程在这里
观察上下两式,发现题目求的是d = 1情况下的裴蜀定理
然后,正难则反,将m分解质因数,通过容斥原理来排除
这里的思路可以参考这里
最后是如何实现n个素数组合的枚举,这里可以用到二进制来进行枚举,即每个素数相当与二进制数的一位,1代表枚举这个素数,0则代表不枚举
具体代码实现可以参考这里
AC代码
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
typedef long long LL;
int prime[11];
LL n, m, mm, ans, tot, cnt;
LL pow(LL x, LL y)//x^y
{
LL res = 1;
for(int i = y; i; i >>= 1)
{
if(i & 1) res *= x;
x = x * x;
}
return res;
}
int main()
{
scanf("%lld%lld", &n, &m);
mm = m; tot = 0;
for(int i = 2; i * i <= mm; ++i)
{
if(mm % i == 0)
{
prime[++tot] = i;
while(mm % i == 0) mm /= i;
}
}
if(mm > 1) prime[++tot] = mm;
for(int i = 1; i < 1 << tot; ++i)
{
cnt = 0;
LL tmp = 1;
for(int j = 1; j <= tot; ++j)
{
if(i & 1 << (j - 1))
{
++cnt;
tmp *= prime[j];
}
}
if(cnt & 1) ans += pow(m / tmp, n);
else ans -= pow(m / tmp, n);
}
ans = pow(m, n) - ans;
printf("%lld\n", ans);
return 0;
}