Bazinga 题解

第十四届浙江财经大学程序设计竞赛重现赛-B题

https://www.nowcoder.com/acm/contest/89/B

可能最近，脑子有问题，看见数论题都是秒，学弟问我这题怎么做，结果我沉思了几分钟，居然秒了

首先题意是求[l,r]内与k互质的数的乘积%k的值。

首先我知道两个比较简单的结论

结论1.：gcd(a+k,k)=gcd(a,k)

没错这就是更相减损法，如果你不认识的话，那应该认识这个算法的优化算法——辗转相除法。就是经常用来求GCD的那个算法。

这个算法给出了的一个重要的结论就是与k的互质的数，是以k为周期出现的。因为若gcd(a,k)=1,则gcd(a+k,k)还是1，也就是a与a+k.都与k互质

结论2：a*(a+k)*(a+2k)……*(a+mk)%k=a^(m+1)%k

这个结论比较好证，首先把每项都%k,你就会发现其实就是 a*a*a*a*a……*a=a^(m+1)%k

题解思路

枚举[1,k-1]内每个与互质的数a，然后算a,a+k,a+2k,a+3k,……中有多少个数落在[l,r].这个比较好求我就不详细讲了，假设最后算出来是m个。

则a的贡献就是a^m%k.。然后把所有数的贡献累乘一下就行了。

不论你是怎么实现的，写得再搓复杂度一般也在k*log(n)以内

因为时间给的比较多，所以我没加任何优化，直接硬上就过了。

#include <stdio.h>
#include<stack>
#include<algorithm>
using namespace std;
#define MAXN 2000005
#define  MAX 0
using namespace std;
long long qpow(long long a,long long n,long long mod)
{
    long long ans=1;
    while(n)
    {
        if(n&1)
            ans=ans*a%mod;
        a=a*a%mod;
        n>>=1;
    }
    return ans;
}
int cal(long long l,long long r,int k)
{
    int i;
    long long ans=1,m;
    for(i=2; i<k&&i<=r; i++)
    {
        if(__gcd(i,k)==1)
        {
            m=(r+k-i)/k-(l+k-i-1)/k;
            ans=ans*qpow(i,m,k)%k;
        }
    }
    return ans;
}
int main()
{
    long long l,r,k;
    int i,j,t,cas=1;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lld%lld%lld",&l,&r,&k);
        printf("Case #%d: ",cas++);
        if(k==1)
            puts("0");
        else
            printf("%d
",cal(l,r,k));

    }
    return 0;
}