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=;
     while(n)
     {
         if(n&)
             ans=ans*a%mod;
         a=a*a%mod;
         n>>=;
     }
     return ans;
 }
 int cal(long long l,long long r,int k)
 {
     int i;
     long long ans=,m;
     for(i=; i<k&&i<=r; i++)
     {
         if(__gcd(i,k)==)
         {
             m=(r+k-i)/k-(l+k-i-)/k;
             ans=ans*qpow(i,m,k)%k;
         }
     }
     return ans;
 }
 int main()
 {
     long long l,r,k;
     int i,j,t,cas=;
     scanf("%d",&t);
     while(t--)
     {
         scanf("%lld%lld%lld",&l,&r,&k);
         printf("Case #%d: ",cas++);
         if(k==)
             puts("");
         else
             printf("%d\n",cal(l,r,k));

     }
     return ;
 }