HDU 2685 I won’t tell you this is about number theory 奇奇怪怪的数论结论

HDU 2685 I won't tell you this is about number theory

标签

• 数学——数论——奇奇怪怪的数论结论

前言

• 我的csdn和博客园是同步的，欢迎来访danzh-博客园~

简明题意

• 给定a,n,m,mod，求：
  \[gcd(a^m-1,a^n-1)\%mod\]

思路

• 没啥思路，就是关于gcd有一个神奇的结论:
  \[gcd(x^a -y^a,x^b-y^b)=x^{gcd(a,b)}-y^{gcd(a,b)}\]
• 把这个结论套进去，就能得到原式=\(a^{gcd(m,n)}-1\)，算出来就可以了。

注意事项

• 取模操作一定要当心减法。有减法，就要先取模，再加模，再取模

总结

• 神奇的结论：
  \[gcd(x^a-y^a,x^b-y^b)=x^{gcd(a,b)}-y^{gcd(a,b)}\]

AC代码

#include<cstdio>
#include<queue>
#include<cmath>
#include<ctime>
using namespace std;

const long long maxn = 1e5 + 10;

long long gcd(long long a, long long b)
{
    if (b == 0) return a;
    return gcd(b, a % b);
}

int ksm(int a, int b, int mod)
{
    int ans = 1, base = a;
    while (b)
    {
        if (b & 1)
            ans = 1ll * ans * base % mod;
        b >>= 1;
        base = 1ll * base * base % mod;
    }
    return ans;
}

void solve()
{
    int t;
    scanf("%d", &t);
    while (t--)
    {
        int a, n, m, mod;
        scanf("%d%d%d%d", &a, &n, &m, &mod);
        printf("%d\n", ((ksm(a, gcd(m, n), mod) - 1) % mod + mod) % mod);
    }
}

int main()
{
    //freopen("Testin.txt", "r", stdin);
    //freopen("Testout.txt", "w", stdout);
    solve();
    return 0;
}

转载于:https://www.cnblogs.com/danzh/p/11402247.html