跳蚤
Time Limit: 1000MS | Memory Limit: 10000K | |
Total Submissions: 7841 | Accepted: 2275 |
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张不同的卡片。现在的问题是,在这所有的卡片中,有多少张可以完成任务。
比如当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
Sample 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)
(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)
Source
本题可以转化为判断ax+by+cz+....pM=1有整数解的情况,根据扩展欧几里德定理,则a,b,c...M等数的公约数应该为1.其中1<=a...M<=M,我们可以求出不符合要求的数列,在不符合要求的数列中,前N个数的公约数中至少有一个是M的质因子。设公约数中包含有n个质因子的数列个数为t(n),则根据容斥原理我们可以得到最终的卡片数为M^N - t(1)+t(2)...+(-1)^(k-1)t(k).另外由于本题的测试数据较弱加上本题的重点是数论问题的应用,所以没使用高精度计算。#include <cstdio>
#define MAX_P 120000 + 10
__int64 prime[MAX_P],pchoose[MAX_P];
__int64 primesize;
int N;
__int64 M;
__int64 res = 0;
__int64 Pow(__int64 x, int y)
{
__int64 temp = 1;
for(int i = 0; i < y; i++)
temp = x * temp;
return temp;
}
void getPrimeFactor(__int64 m) //获取m的质因子
{
primesize = 0;
for(__int64 i = 2; i * i <= m; i++)
{
if(m % i == 0)
{
prime[primesize++] = i;
while(m % i == 0)
m = m / i;
}
}
if(m != 1)
prime[primesize++] = m;
}
void getCount(__int64 pos, __int64 curnum, __int64 totalnum)
{
if(curnum == totalnum)
{
__int64 temp = M;
for(__int64 i = 0; i < curnum; i++)
temp /= pchoose[i];
res = res + Pow(temp, N);
}
else
{
for(__int64 i = pos; i < primesize; i++)
{
pchoose[curnum] = prime[i];
getCount(i + 1, curnum + 1, totalnum);
}
}
}
void solve()
{
getPrimeFactor(M);
__int64 totalCount = Pow(M, N);
for(__int64 i = 1; i <= primesize; i++)
{
res = 0;
getCount(0, 0, i);
if(i % 2 == 1)
totalCount -= res;
else
totalCount += res;
}
printf("%I64d\n", totalCount);
}
int main()
{
scanf("%d %I64d", &N, &M);
solve();
return 0;
}