标签：space info span -- 推导 个数 end mamicode 大于
题意：
一个数论题，要求满足如下等式的x有多少个。
思路：
当时比赛是，队伍看到这个题，也没有做太多的思考，就是无从下手，几乎放弃。但是看到学校另外两支队伍都过了这个题，感觉自己还是好菜。
打表可以发现，当a为奇数的时候答案为1。当a为偶数的时候，x一定也是偶数，这个还是比较明显的。
对左边进行推导，因为a为偶数，设a=2*t，所以a^x=2^x*t^x,所以当x大于p时，这个求余之后一定为0。由于p很小，可以直接暴力求解，所以对于右边直接考虑x^a求余之后为0的x的情况。由于x为偶数。我们再次对x分解，设x=2^k*t,就是把x写成一个2的 次幂和t的乘积的形式，所以x^a=2^(ka)*t^a，所以只需要ka>=p,就可以满足条件，而且只有ka>=p可以满足条件，因为t^a不可能可以整除，所以现在要求k>=p/a，对p/a取上整即可，就是考察有多少个t，在1-up的范围内（up为x的上界），个数为up/ceil(p/a)/2^k。当然最后的答案要减去1-p范围内的数，因为这一部分是特殊考虑的。
1 #include <iostream> 2 #include<cmath> 3 #include<stdio.h> 4 //a^x 与x^a 5 // 1<=x<=2^p 6 using namespace std; 7 #define ll long long 8 ll ppow(int x,int n,ll mod) 9 { 10 ll res=1; 11 ll base=x; 12 while(n) 13 { 14 if(n&1) 15 res=(res*base)%mod; 16 base=(base*base)%mod; 17 n>>=1; 18 } 19 return res%mod; 20 } 21 int main() 22 { 23 int a,p; 24 int T; 25 while(~scanf("%d",&T)) 26 { 27 while(T--) 28 { 29 scanf("%d%d",&a,&p); 30 31 if(a%2) 32 { 33 cout<<1<<endl; 34 continue; 35 } 36 //x>=p 37 ll up=(1LL<<p); 38 ll ans=0; 39 for(int i=1;i<=p;i++){ 40 if(ppow(a,i,up)==ppow(i,a,up))ans++; 41 } 42 ll tem=ceil(p/(a+0.0)); 43 ans+=(up>>tem)-(p>>tem); 44 cout<<ans<<endl; 45 } 46 } 47 48 return 0; 49 }
2019山东省赛K - Happy Equation ZOJ - 4123 题解
标签：space info span -- 推导 个数 end mamicode 大于
原文地址：https://www.cnblogs.com/gzr2018/p/10883184.html
