結果

問題 No.2262 Fractions
コンテスト
ユーザー chro_96
提出日時 2023-04-08 00:51:39
言語 C
(gcc 13.3.0)
結果
TLE  
実行時間 -
コード長 4,119 bytes
コンパイル時間 238 ms
コンパイル使用メモリ 30,848 KB
実行使用メモリ 15,872 KB
最終ジャッジ日時 2024-10-02 21:52:26
合計ジャッジ時間 6,702 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample TLE * 1
other -- * 45
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>

int main () {
  int t = 0;
  
  int res = 0;
  
  int np[300001] = {};
  int nf[300001] = {};
  int pcnt[300001] = {};
  
  np[1] = 1;
  nf[1] = 1;
  for (int i = 2; i <= 300000; i++) {
    if (np[i] <= 0) {
      pcnt[i] = 1;
      for (int p = 2; i*p <= 300000; p++) {
        np[i*p] = 1;
        pcnt[i*p]++;
        if (p%i == 0) {
          nf[i*p] = 1;
        }
      }
    }
  }
  
  res = scanf("%d", &t);
  while (t > 0) {
    long long n = 0LL;
    long long k = 0LL;
    long long ans[2] = { 1LL, 0LL };
    long long l[2] = { 0LL, 1LL };
    res = scanf("%lld", &n);
    res = scanf("%lld", &k);
    while (l[0]+ans[0] <= n && l[1]+ans[1] <= n) {
      long long nxt[2] = { l[0]+ans[0], l[1]+ans[1] };
      long long cnt = 0LL;
      for (long long i = 1LL; i <= n; i += 1LL) {
        long long d = (nxt[0]*i)/nxt[1];
        if (d > n) {
          d = n;
        }
        cnt += d;
        if (nf[(int)i] <= 0) {
          long long tmp = 0LL;
          for (long long j = 1LL; i*j <= n; j += 1LL) {
            long long d = (nxt[0]*i*j)/nxt[1];
            if (d > n) {
              d = n;
            }
            tmp += d/i;
          }
          if (pcnt[(int)i]%2 == 0) {
            cnt += tmp;
          } else {
            cnt -= tmp;
          }
        }
      }
      if (cnt < k) {
        l[0] += ans[0];
        l[1] += ans[1];
      } else {
        ans[0] += l[0];
        ans[1] += l[1];
      }
      /*if (cnt < k) {
        long long t[2] = { 1LL, n };
        if (t[1] > (n-l[0])/ans[0]) {
          t[1] = (n-l[0])/ans[0];
        }
        if (ans[1] > 0LL && t[1] > (n-l[1])/ans[1]) {
          t[1] = (n-l[1])/ans[1];
        }
        t[1] += 1LL;
        while (t[1]-t[0] > 1LL) {
          long long nxt_v = (t[0]+t[1])/2LL;
          long long nxt[2] = { l[0]+nxt_v*ans[0], l[1]+nxt_v*ans[1] };
          long long cnt = 0LL;
          for (long long i = 1LL; i <= n; i += 1LL) {
            long long d = (nxt[0]*i)/nxt[1];
            if (d > n) {
              d = n;
            }
            cnt += d;
            if (nf[(int)i] <= 0) {
              long long tmp = 0LL;
              for (long long j = 1LL; i*j <= n; j += 1LL) {
                long long d = (nxt[0]*i*j)/nxt[1];
                if (d > n) {
                  d = n;
                }
                tmp += d/i;
              }
              if (pcnt[(int)i]%2 == 0) {
                cnt += tmp;
              } else {
                cnt -= tmp;
              }
            }
          }
          if (cnt < k) {
            t[0] = nxt_v;
          } else {
            t[1] = nxt_v;
          }
        }
        l[0] += t[0]*ans[0];
        l[1] += t[0]*ans[1];
      } else {
        long long t[2] = { 1LL, n };
        if (l[0] > 0LL && t[1] > (n-ans[0])/l[0]) {
          t[1] = (n-ans[0])/l[0];
        }
        if (t[1] > (n-ans[1])/l[1]) {
          t[1] = (n-ans[1])/l[1];
        }
        t[1] += 1LL;
        while (t[1]-t[0] > 1LL) {
          long long nxt_v = (t[0]+t[1])/2LL;
          long long nxt[2] = { nxt_v*l[0]+ans[0], nxt_v*l[1]+ans[1] };
          long long cnt = 0LL;
          for (long long i = 1LL; i <= n; i += 1LL) {
            long long d = (nxt[0]*i)/nxt[1];
            if (d > n) {
              d = n;
            }
            cnt += d;
            if (nf[(int)i] <= 0) {
              long long tmp = 0LL;
              for (long long j = 1LL; i*j <= n; j += 1LL) {
                long long d = (nxt[0]*i*j)/nxt[1];
                if (d > n) {
                  d = n;
                }
                tmp += d/i;
              }
              if (pcnt[(int)i]%2 == 0) {
                cnt += tmp;
              } else {
                cnt -= tmp;
              }
            }
          }
          if (cnt >= k) {
            t[0] = nxt_v;
          } else {
            t[1] = nxt_v;
          }
        }
        ans[0] += t[0]*l[0];
        ans[1] += t[0]*l[1];
      }*/
    }
    if (ans[1] == 0LL) {
      printf("-1\n");
    } else {
      printf("%lld/%lld\n", ans[0], ans[1]);
    }
    t--;
  }
  
  return 0;
}
0