QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#570638 | #9309. Graph | lprdsb | AC ✓ | 2842ms | 797780kb | C++14 | 4.2kb | 2024-09-17 16:50:25 | 2024-09-17 16:50:39 |
Judging History
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