QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#570638#9309. GraphlprdsbAC ✓2842ms797780kbC++144.2kb2024-09-17 16:50:252024-09-17 16:50:39

Judging History

你现在查看的是最新测评结果

  • [2024-09-17 16:50:39]
  • 评测
  • 测评结果:AC
  • 用时:2842ms
  • 内存:797780kb
  • [2024-09-17 16:50:25]
  • 提交

answer

#include<bits/stdc++.h>
#define For(i, a, b) for(LL i = (a), en = (b); i <= en; ++i)
#define Rof(i, a, b) for(LL i = (a), en = (b); i >= en; --i)
#define Tra(u, i) for(int i = hd[u]; ~i; i = e[i].net)
#define cst const
#define LL long long
#define DD double
#define LD long double
#define pb push_back
#define mp make_pair
#define fir first
#define sec second
#define inf 0x3f3f3f3f
#define Inf 0x3f3f3f3f3f3f3f3f
#define eps 1e-12
#define maxn 1000000
using namespace std;

LL pri[maxn + 5], cnt = 0, vis[maxn + 5];

template<class T>
T fp(T x, LL y) {
    T asi = 1;
    while(y) {
        if(y & 1) asi *= x;
        x *= x;
        y >>= 1;
    }
    return asi;
}

template<int P>
struct Mint {
    int x;
    Mint() { x = 0; }
    Mint(LL _x) : x{ norm(_x % P) } {}
    int norm(LL _x) {
        if(_x < 0) _x += P;
        if(_x >= P) _x -= P;
        return _x;
    }
    Mint inv() { return fp(*this, P - 2); }
    Mint operator - () { return { P - x }; }
    Mint& operator += (Mint _x) & {
        x = norm(x + _x.x);
        return *this;
    }
    Mint& operator -= (Mint _x) & {
        x = norm(x - _x.x);
        return *this;
    }
    Mint& operator *= (Mint _x) & {
        x = 1ll * x * _x.x % P;
        return *this;
    }
    Mint& operator /= (Mint _x) & { return *this *= _x.inv(); }
    friend Mint operator + (Mint _x, Mint _y) { return _x += _y; }
    friend Mint operator - (Mint _x, Mint _y) { return _x -= _y; }
    friend Mint operator * (Mint _x, Mint _y) { return _x *= _y; }
    friend Mint operator / (Mint _x, Mint _y) { return _x /= _y; }
    friend bool operator == (Mint _x, Mint _y) { return _x.x == _y.x; }
};

const int P = 998244353;
using Z = Mint<P>;

class Sieve {
    public:
    LL to[maxn + 5], tot = 0, k_n[maxn + 5];
    LL N;
    LL sum_f_p[maxn + 5];
    unordered_map<LL, int> id;
    vector<bool> vG[maxn + 5];
    vector<LL> G[maxn + 5];
    LL get_g(LL x) {
        //
        return 1;
    }
    Sieve() {}
    LL get_G(int n_id, int k) {
        if(pri[k] > to[n_id]) return get_g(1);
        if(k == 1) return to[n_id];
        if(1ll * pri[k - 1] * pri[k - 1] > to[n_id]) return get_G(n_id, k_n[n_id] + 1) - (sum_f_p[k - 1] - sum_f_p[k_n[n_id]]);
        return G[n_id][k];
    }
    void init(LL _N) {
        N = _N;
        tot = 0;
        For(i, 1, N) {
            ++tot;
            id[to[tot] = N / i] = tot;
            i = N / (N / i);
        }
        int now = cnt;
        For(i, 1, tot) {
            while(1ll * pri[now] * pri[now] > to[i]) now--;
            G[i].clear();
            G[i].resize(now + 5, 0);
            vG[i].clear();
            vG[i].resize(now + 5, 0);
            k_n[i] = now;
        }
        For(i, 1, cnt) sum_f_p[i] = sum_f_p[i - 1] + get_g(pri[i]);
        For(i, 2, k_n[1] + 1) {
            For(j, 1, tot) {
                if(1ll * pri[i - 1] * pri[i - 1] > to[j]) break;
                vG[j][i] = 1;
                G[j][i] = get_G(j, i - 1) - get_g(pri[i - 1]) * get_G(id[to[j] / pri[i - 1]], i - 1);
            }
        }
    }
    LL get_s(LL n) {
        return k_n[id[n]] + get_G(id[n], k_n[id[n]] + 1);
    }
} sieve;

template <class T>
void read(T &x) {
    char ch;
    bool ok;
    for(ok = 0, ch = getchar(); !isdigit(ch); ch = getchar()) if(ch == '-') ok = 1;
    for(x = 0; isdigit(ch); x = x * 10 + ch - '0', ch = getchar());
    if(ok) x = -x;
}

int main() {
    // freopen("in.txt", "r", stdin);
    For(i, 2, maxn) {
        if(!vis[i]) pri[++cnt] = i;
        for(int j = 1; j <= cnt && i * pri[j] <= maxn; ++j) {
            vis[i * pri[j]] = 1;
            if(i % pri[j] == 0) break;
        }
    }
    LL n = 4;
    read(n);
    sieve.init(n);
    // if(n == 84848352911ll) {
    //     cout << "asd" << endl;
    //     return 0;
    // }
    Z res = 1;
    for(LL i = 1; i <= n; ++i) {
        LL now = n / i, to = n / now;
        if(now <= 2) break;
        LL cntp = sieve.get_s(now) - sieve.get_s(now / 2);
        LL cnt = cntp + 1 + (now - cntp - 1 != 0);
        res *= fp(fp(Z{ now }, cnt - 2) * Z { (now - cntp - 1) ? (now - cntp - 1) : 1 }, to - i + 1);
        i = to;
    }
    cout << res.x << endl;
    return 0;
}


这程序好像有点Bug,我给组数据试试?

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 14ms
memory: 78972kb

input:

4

output:

8

result:

ok answer is '8'

Test #2:

score: 0
Accepted
time: 16ms
memory: 81440kb

input:

2

output:

1

result:

ok answer is '1'

Test #3:

score: 0
Accepted
time: 12ms
memory: 83492kb

input:

123

output:

671840470

result:

ok answer is '671840470'

Test #4:

score: 0
Accepted
time: 8ms
memory: 81404kb

input:

233

output:

353738465

result:

ok answer is '353738465'

Test #5:

score: 0
Accepted
time: 24ms
memory: 80444kb

input:

5981

output:

970246821

result:

ok answer is '970246821'

Test #6:

score: 0
Accepted
time: 8ms
memory: 81528kb

input:

86422

output:

897815688

result:

ok answer is '897815688'

Test #7:

score: 0
Accepted
time: 8ms
memory: 79324kb

input:

145444

output:

189843901

result:

ok answer is '189843901'

Test #8:

score: 0
Accepted
time: 16ms
memory: 79748kb

input:

901000

output:

819449452

result:

ok answer is '819449452'

Test #9:

score: 0
Accepted
time: 7ms
memory: 82632kb

input:

1000000

output:

113573943

result:

ok answer is '113573943'

Test #10:

score: 0
Accepted
time: 14ms
memory: 83432kb

input:

23333333

output:

949849384

result:

ok answer is '949849384'

Test #11:

score: 0
Accepted
time: 20ms
memory: 89028kb

input:

102850434

output:

604886751

result:

ok answer is '604886751'

Test #12:

score: 0
Accepted
time: 73ms
memory: 112952kb

input:

998244353

output:

0

result:

ok answer is '0'

Test #13:

score: 0
Accepted
time: 75ms
memory: 117572kb

input:

1000000007

output:

318420284

result:

ok answer is '318420284'

Test #14:

score: 0
Accepted
time: 120ms
memory: 134268kb

input:

2147483547

output:

688759898

result:

ok answer is '688759898'

Test #15:

score: 0
Accepted
time: 266ms
memory: 177220kb

input:

5120103302

output:

116870489

result:

ok answer is '116870489'

Test #16:

score: 0
Accepted
time: 842ms
memory: 323072kb

input:

19834593299

output:

523663743

result:

ok answer is '523663743'

Test #17:

score: 0
Accepted
time: 1726ms
memory: 546628kb

input:

52500109238

output:

195086665

result:

ok answer is '195086665'

Test #18:

score: 0
Accepted
time: 2476ms
memory: 724684kb

input:

84848352911

output:

107959260

result:

ok answer is '107959260'

Test #19:

score: 0
Accepted
time: 2828ms
memory: 795056kb

input:

99824435322

output:

0

result:

ok answer is '0'

Test #20:

score: 0
Accepted
time: 2842ms
memory: 797780kb

input:

99999999354

output:

316301711

result:

ok answer is '316301711'

Test #21:

score: 0
Accepted
time: 2816ms
memory: 795788kb

input:

100000000000

output:

396843576

result:

ok answer is '396843576'

Extra Test:

score: 0
Extra Test Passed