HDU 4335はN何ですか?問題の(インデックスループ部)ソリューション

質問の意味:

どのようにいくつか尋ねる\(N- \)を満足\({!} N-N- ^ \当量のB \ MOD P \土地\ n-は\で\ [1、M]) データ範囲:\(M \ ^ {64} leq2 -1、P \ leq1e5 \)

アイデア:

この質問降順、オイラーを使用することは明らかである(\!N-)\未満です\(\ varphi(P)\ ) 直接の暴力、考えられている\(N-!\ NEQ 0 \ MODの\ varphi(P-)\)にも直接的な暴力。
(!N \当量0 \ \ MOD \ varphi(P)\) 明らかに、この場合には、素数の定数である\(\ varphi(P)\)ピジョンホール定理から調製し、\(X \)時定数、(\ 1 ^のx、2 ^ X、モッズのp \ \ \ドット、N ^ X、\ドット) として循環部\(\ varphi(P)\ ) を押しサイクルセクションではどのようなことができます行います。場合、注意\(B = 0、M = 2 ^ {64} -1、p = 1 \) 、応答バースト\(ロングロング\)


コード:

#include<map>
#include<set>
#include<queue>
#include<cmath>
#include<stack>
#include<ctime>
#include<vector>
#include<cstdio>
#include<string>
#include<cstring>
#include<sstream>
#include<iostream>
#include<algorithm>
typedef long long ll;
typedef unsigned long long ull;
using namespace std;
const int maxn = 1e5 + 5;
const int MAXM = 3e6;
const ll MOD = 998244353;
const ull seed = 131;
const int INF = 0x3f3f3f3f;
ull euler(ull n){
    ull res = n, a = n;
    for(int i = 2; i * i <= a; i++){
        if(a % i == 0){
            res = res / i * (i - 1);
            while(a % i == 0) a/= i;
        }
    }
    if(a > 1) res = res / a * (a - 1);
    return res;
}
ull ppow(ull a, ull b, ull mod){
    ull ret = 1;
    while(b){
        if(b & 1) ret = ret * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return ret;
}
ull rec[maxn];
int main(){
    int T, ca = 1;
    scanf("%d", &T);
    while(T--){
        ull b, p, m;
        scanf("%I64u%I64u%I64u", &b, &p, &m);
        if(b == 0 && p == 1){
            if(m == 18446744073709551615ULL)
                printf("Case #%d: 18446744073709551616\n", ca++);
            else
                printf("Case #%d: %I64u\n", ca++, m + 1);
            continue;
        }
        ull phi = euler(p);
        ull ans = 0, fac = 1;
        ull i = 1;
        if(b == 0) ans++;
        for(i = 1; i <= m; i++){
            if(fac * i >= phi) break;
            fac = fac * i;
            if(ppow(i, fac, p) == b) ans++;
        }
        for(; i <= m; i++){
            if(fac * i % phi == 0) break;
            fac = fac * i % phi;
            if(ppow(i, fac + phi, p) == b) ans++;
        }
        if(i <= m){
            ull cnt = 0;
            for(int j = 1; j <= p; j++){
                rec[j] = ppow(j, phi, p);
                if(rec[j] == b) cnt++;
            }
            for(; i <= p && i <= m; i++){
                if(rec[i] == b) ans++;
            }
            if(i <= m){
                ull rest = m - p;
                ans += rest / p * cnt;
                rest -= rest / p * p;
                for(i = 1; i <= rest; i++){
                    if(rec[i] == b) ans++;
                }
            }
        }
        printf("Case #%d: %I64u\n", ca++, ans);
    }
    return 0;
}

おすすめ

転載: www.cnblogs.com/KirinSB/p/11536602.html