QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#850942 | #9962. Diminishing Fractions | Zawos | RE | 0ms | 0kb | C++23 | 3.2kb | 2025-01-10 13:23:44 | 2025-01-10 13:23:45 |
answer
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#define FOR(i,a,b) for (int i = (a); i < (b); i++)
using namespace std;
using namespace __gnu_pbds;
using ll=long long;
using ld=long double;
using vi=vector<int>;
template<class T> using oset =tree<T, null_type, less<T>, rb_tree_tag,tree_order_statistics_node_update> ;
//上
const int N = 50001;
int wh[N];
ll F[5217][5217];
int euclid(int a, int b, ll &x, ll &y) {
if (!b) return x = 1, y = 0, a;
int d = euclid(b, a % b, y, x);
return y -= a/b * x, d;
}
int main(){
ios_base::sync_with_stdio(false);
cin.tie(NULL);
memset(wh,-1,sizeof(wh));
FOR(i,0,5217) FOR(j,0,5217) F[i][j] = 1;
vector<int> prime(N+1,1);
for(int i = 2; i <= N; i++){
if(prime[i] == 1) for(int j = i; j <=N;j+=i) prime[j] = i;
}
vector<pair<int,int>> pot;
for(int i = 2; i <= N;i++){
unordered_map<int,int> fac;
int x = i;
int r = 1;
while( x > 1){
fac[prime[x]]++;
r *= prime[x];
x /= prime[x];
if(fac.size() > 1) break;
}
if(fac.size() == 1){
for(auto &u: fac) pot.push_back({u.first,r});
}
}
for(int i = 0; i < 5217;i++) wh[pot[i].second] = i;
for(int i = 0; i <5217;i++){
for(int j = 0; j < 5217;j++){
if(pot[i].first == pot[j].first){
if(j >0) F[i][j] = F[i][j-1];
}else{
if(j >0){
F[i][j] = F[i][j-1]*pot[j].first;
}else{
F[i][j] *= pot[j].first;
}
}
}
}
ll y;
int t;
cin >>t;
unordered_map<int,int> mp;
vector<pair<int,int>> q;
for(int tt = 0; tt < t; tt++){
int n;
cin >> n;
q.push_back({n,tt});
}
sort(q.begin(),q.end());
vector<pair<int,string>> res;
int last = 1;
int lastp = 0;
for(int tt = 0; tt < t; tt++){
if(q[tt].first == 1){
res.push_back({q[tt].second,"1/1\n"});
continue;
}else if(q[tt].first == 2){
res.push_back({q[tt].second,"1/2\n"});
continue;
}
for(int i = last; i <= q[tt].first; i++){
if(wh[i] != -1){
mp[pot[wh[i]].first] = wh[i];
lastp = wh[i];
}
last++;
}
double ans = 0.0;
int co = 0;
res.push_back({q[tt].second,""});
for(auto &u: mp){
ll a = F[wh[u.second]][lastp];
euclid(a,pot[u.second].second,a,y);
if(a <0) a += pot[u.second].second;
if(co > 0)(res.back()).second+="+";
(res.back()).second+=to_string(a)+"/"+to_string(pot[u.second].second);
co++;
ans += (double)a/(double)pot[u.second].second;
}
ans += 0.5;
int xx = floor(ans);
(res.back()).second+="-"+to_string(xx)+"/1\n";
}
sort(res.begin(),res.end());
for(auto &u: res) for(auto&e:u.second) printf("%c",e);
}
詳細信息
Test #1:
score: 0
Runtime Error
input:
2 3 6