結果
| 問題 | No.2318 Phys Bone Maker |
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-05-26 22:22:27 |
| 言語 | C++17 (gcc 15.2.0 + boost 1.89.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 993 bytes |
| 記録 | |
| コンパイル時間 | 4,778 ms |
| コンパイル使用メモリ | 269,240 KB |
| 最終ジャッジ日時 | 2025-02-13 07:11:53 |
|
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 WA * 1 |
| other | AC * 6 WA * 39 |
ソースコード
#include<bits/stdc++.h>
using namespace std;
#include<atcoder/all>
using namespace atcoder;
using ll = long long;
ll x;
vector<long long> divisor(long long n) {
vector<long long> ret;
for (long long i = 1; i * i <= n; i++) {
if (n % i == 0) {
ret.push_back(i);
if (i * i != n) ret.push_back(n / i);
}
}
sort(ret.begin(), ret.end()); // 昇順に並べる
return ret;
}
using mit = modint998244353;
unordered_map<ll,ll> mp;
unordered_map<ll,mit> dp;
unordered_map<ll,vector<ll>> divisor_mp;
void solve(){
auto v = divisor(x);
for(auto &i:v){
divisor_mp[i] = divisor(i);
}
dp[1] = 1;
for(auto &i:v){
auto &v2 = divisor_mp[x/i];
for(auto &j:v2){
if(j==1)continue;
dp[i*j] += dp[i]*(divisor_mp[j/gcd(i,j)].size()-1);
}
/*
for(auto &j:v){
cout<<i<<':'<<j<<':'<<dp[j].val()<<endl;;
}
*/
}
cout<<dp[x].val()<<endl;
}
signed main(){
cin.tie(nullptr);
ios::sync_with_stdio(false);
cin >> x;
solve();
}