QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#570104#9309. GraphjeffqiRE 2ms4792kbC++172.5kb2024-09-17 13:50:262024-09-17 13:50:33

Judging History

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

  • [2024-09-17 13:50:33]
  • 评测
  • 测评结果:RE
  • 用时:2ms
  • 内存:4792kb
  • [2024-09-17 13:50:26]
  • 提交

answer

#include<bits/stdc++.h>
#include<ext/pb_ds/assoc_container.hpp>
#include<ext/pb_ds/hash_policy.hpp>
#define ll long long
#define all(v) v.begin(),v.end()
#define sz(v) ((ll)v.size())
#define V vector
#define vi V<int>
#define vll V<ll>
#define eb emplace_back
#define pii pair<int,int>
#define pll pair<ll,ll>
#define fi first
#define se second
#define A array
#define pb push_back
#define mset multiset
#define gpumap __gnu_pbds::gp_hash_table
#define ccumap __gnu_pbds::cc_hash_table
#define ui unsigned int
#define ull unsigned ll
#define i128 __int128
#define cerr if (test) cerr
#define freopen if (test) freopen
#define whtest if (test)
using namespace std;

const int test = 0;

namespace jhsy {
	constexpr int P = 998244353,N = 1e6;
	vi p;
	void init(const int n = N) {
		A<bool,N+1> isp; isp.fill(1);
		for (int i = 2; i <= n; i++) {
			if (isp[i]) {
				p.eb(i);
			}
			for (auto x:p) {
				if (i*x > n) {
					break;
				}
				isp[i*x] = 0;
				if (!(i%x)) {
					break;
				}
			}
		}
	}
	ll power(ll b,ll k) {
		ll a = 1;
		b = (b%P+P)%P;
		k = (k%(P-1)+(P-1))%(P-1);
		for (; k; b = b*b%P, k >>= 1) {
			if (k&1) {
				a = a*b%P;
			}
		}
		return a;
	}
	void main() {
		ll n;
		cin >> n;
		ll ans = 1;
		const ll B = sqrt(1.L*n+1);
		vll mp(B*2);
		auto id = [&](ll x) {
			return x <= B ? x : B+n/x;
		};
		{
			vll vec;
			for (ll lst = 0,cur; lst < n; lst = cur) {
				ll k = n/(lst+1); cur = n/k;
				vec.eb(cur);
				assert(id(cur) < sz(mp));
			}
			reverse(all(vec));
			for (auto x:vec) {
				mp[id(x)] = x-1;
			}
			for (int k = 0; k < sz(p); k++) {
				ll t = 1ll*p[k]*p[k];
				if (t > n) {
					break;
				}
				for (auto x:vec) {
					if (t > x) {
						break;
					}
					mp[id(x)] -= mp[id(x/p[k])]-k;
				}
			}
		}
		auto get = [&](ll x) {
			ll res = mp[id(x)]-mp[id(x/2)];
			if (x/2 < 2 && x >= 2) {
				res--;
			}
			return res;
		};
		auto calc = [&](ll x) -> ll {
			if (x == 1) {
				return 1;
			}
			ll p = get(x)+1,q = x-p;
			return q%P*power(p+q,p-1)%P;
		};
		for (ll lst = 0,cur; lst < n; lst = cur) {
			ll k = n/(lst+1); cur = n/k;
			ans = ans*power(calc(k),cur-lst)%P;
		}
		cout << ans << '\n';
	}
}

int main() {
	freopen("test.in", "r", stdin);
	freopen("test.out", "w", stdout);
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	cout << fixed << setprecision(20);

	jhsy::init();

	int T = 1;
//	cin >> T;
	while (T--) {
		jhsy::main();
	}

	return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 2ms
memory: 4792kb

input:

4

output:

8

result:

ok answer is '8'

Test #2:

score: -100
Runtime Error

input:

2

output:


result: