精度问题
查看原帖
精度问题
383395
WaterSun楼主2024/2/22 23:33

精度被卡了,答案能算对,但是转成分数就会被卡精度/ll

#include <bits/stdc++.h>
#define re register
#define fst first
#define snd second
#define double long double

using namespace std;

typedef pair<int,int> pii;
const int N = 4e4 + 10;
double eps = 1e-10;
int n,k;
int f[N];

inline int read(){
    int r = 0,w = 1;
    char c = getchar();
    while (c < '0' || c > '9'){
        if (c == '-') w = -1;
        c = getchar();
    }
    while (c >= '0' && c <= '9'){
        r = (r << 3) + (r << 1) + (c ^ 48);
        c = getchar();
    }
    return r * w;
}

inline pii get(double x){
    int a = 0,b = 1;
    for (re int i = 1;i <= n;i++){
        int t = x * i;
        // cout << i << " " << t << " " << a << " " << b << " ???\n";
        if (fabs(1.0 * a / b - x) > fabs(1.0 * t / i - x)) a = t,b = i;
    }
    return {a,b};
}

inline bool check(double x){
    int num = 0;
    for (re int i = 1;i <= n;i++){
        f[i] = i * x;
        for (re int j = 2;j * j <= x;j++){
            if (i % j == 0){
                f[i] -= f[j];
                if (j != i / j) f[i] -= f[i / j];
            }
        }
        num += f[i];
    }
    return (num <= k);
}

int main(){
    n = read(),k = read();
    double l = 0.0,r = 1.0;
    while (r - l > eps){
        double mid = (l + r) / 2.0;
        if (check(mid)) l = mid;
        else r = mid;
    }
    // printf("%.3Lf\n",l);
    pii ans = get(l);
    printf("%d %d",ans.fst,ans.snd);
    return 0;
}
2024/2/22 23:33
加载中...