题目:
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
Inputcopy | Outputcopy |
---|---|
2 4 | 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,分解M的质因子,比如当m为36时,质因子有2,3;在大于1小于等于m的,且是2的倍数的数有:2,4,6,8,10,12,14,16,18,20,22,24,26,28,30,32,34,36。是3的倍数的数有:3,6,9,12,15,18,21,24,27,30,33,36。其中6,12,18,24,30,36为重复数字。
所以设m一共有s个质因子,则答案为m^n - 所有是单个质因子倍数组成的n个数的情况数 + 所有是2个质因子乘积倍数组成的n个数的情况数 - ... + (-1)^s * 所有是s个质因子乘积倍数组成的n个数的情况数
代码:
#include<stdio.h>
#include<iostream>
#include<vector>
using namespace std;
typedef long long ll;
bool c[100];
ll ans = 0;
ll n,m;
ll yy[15];
int len;
ll T = 0;
vector<int> p;
ll pow(ll a,ll n)
{
ll ans = 1;
while(n)
{
if (n&1) ans *= a;
a *= a;
n >>= 1;
}
return ans;
}
void dfs(ll ss,int step,int stop)
{
if (step == stop)
{
if(ss <= m )
T += pow(m / ss,n);
}
else
{
for (int i = 0; i < len; i++)
{
if(!c[i])
{
c[i] = true;
dfs(ss * p[i],step + 1,stop);
c[i] = false;
}
}
}
}
int main()
{
for(int i = 1; i < 15; i++){
yy[i] = 1;
for(int j = 2; j <= i; j++)
{
yy[i] *= j;
}
}
scanf("%lld%lld",&n,&m);
ll tmp = m;
for (int i = 2; i * i <= m; i ++)
{
if (m % i == 0)
{
p.push_back(i);
ans -= pow(tmp / i,n);
while (m % i == 0)
{
m /= i;
}
}
}
if(m > 1) {
ans -= pow(tmp / m,n);
p.push_back(m);
}
len = p.size();
m = tmp;
ans += pow(m,n);
for (int i = 2; i <= len; i++){
T = 0;
dfs(1,0,i);
T /= yy[i];
if(i%2==0) ans += T;
else ans -= T;
}
cout <<ans;
return 0;
}