A/B【费马小定理】

 今天是数论周的最后一天,记得一周前从一个什么是数论?费马定理有什么用都不会的萌新,现在已经能开始写欧拉公式的模版、欧几里得、拓展欧几里得之类的模版了。

要求(A/B)%9973,但由于A很大,我们只给出n(n=A%9973)(我们给定的A必能被B整除,且gcd(B,9973) = 1)。

Input

数据的第一行是一个T,表示有T组数据。 
每组数据有两个数n(0 <= n < 9973)和B(1 <= B <= 10^9)。

Output

对应每组数据输出(A/B)%9973。

Sample Input

2
1000 53
87 123456789

Sample Output

7922
6060

  对于这道题,题目给出gcd(B, 9973)=1即互质,又9973为质数,所以可以看出此处用到费马小定理。

由a=1(mod p),可以知道a^(p-1)=1(mod p)于是有a^(-1)=a^(p-2) (mod p)两边同除以a得到。

完整代码:

#include 
#include 
#include 
#include 
#include 
using namespace std;
typedef long long ll;
ll a,b;
ll s_mi(ll x, ll y, ll mod)
{
    x%=mod;
    //int k=0;
    ll ans=0;
    /*while(y)
    {
        int t=y&1;      //第k位上是否有1?
        ll temp=1;
        y=y>>1;
        if(t)
        {
            for(int i=1; i<=k; i++)
            {
                temp=(temp*(x*x)%mod)%mod;
            }
        }
        ans=(ans+temp)%mod;
        k++;        //放在最尾使得幂次从0开始
    }*/
    ans=x;
    for(int i=1; i

 

你可能感兴趣的:(数论)