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
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)
解题思路
先分析题目给的样例,可以写成10x1+15x2+18x3=1.其中一组解释x1=1,x2=3,x3=-3;根据扩展欧几里得可以知道
a1x1+a2x2+a3x3+…+anxn+(an+1)(xn+1)=1; 也就是gcd(a1,a2,a3,a4,…an+1)=1;根据排列组合可以知道一共有mn种
我们要减去最大因数不是1的,先求出m的质因子,我们可以容易的求出1~m中有几个能被质因子整除的(要用到奇加偶减),
比如题上的列子,18的质因子释2,3。18/2=9,18/3=6,18/3*2=3。18能被质因子整除的就有9+6-3=12种。回归正题,现在N=2,M=18,也就是说第一个数是这12种的其中一个,那么另外一个数就可以是这12种的任意一个(这样的话gcd就不是1),
比如让第一个数是2 第二个数可以是2 ,第三个数是m=18,这样的gcd(2,2,18)!=1,我们可以求出有多少种情况不等于1,然后拿总的情况减去就是符合条件的。题上的列子就可以计算为:182-92-62+32
#include<stdio.h>
#include<iostream>
using namespace std;
typedef long long ll;
const int maxn=1e5;
bool prime[maxn+10];
ll p[maxn+10],k=0;
ll as[maxn];
void ss()
{
prime[0]=prime[1]=1;
for(ll i=2;i<maxn;i++)
{
if(!prime[i])
{
p[k++]=i;
for(ll j=i*i;j<maxn;j+=i)
{
prime[j]=1;
}//这样打表求得比较快
}
}
}
ll ksm(ll a,ll b)
{
ll res=1;
while(b)
{
if(b&1)
res=res*a;
b>>=1;
a=a*a;
}
return res;
}//快速幂
int main()
{
ss();
ll n,m;
scanf("%lld %lld",&n,&m);
int cnt=0;
ll q=m;
ll res=ksm(m,n);
for(ll i=0;i<k&&p[i]*p[i]<=m;i++)
{
if(m%p[i]==0)
{
cnt++;
as[cnt]=p[i]; //求质因子
}
while(m%p[i]==0)
{
m/=p[i];
}
}
if(m>1) cnt++,as[cnt]=m;
for(ll i=1;i<(1<<cnt);i++)
{
ll ans=1,dis=0;
for(ll j=0;j<cnt;j++)
{
if(1&(i>>j))
{
ans*=as[j+1];//注意这点是加1
dis++;
}
}
ll p=ksm(q/ans,n);
if(dis&1)
{
res-=p;//本来是奇加偶减,但因为我们是反着求所以奇减偶加
}
else
{
res+=p;
}
}
printf("%lld\n",res);
return 0;
}