結果

問題 No.2262 Fractions
コンテスト
ユーザー chro_96
提出日時 2023-04-08 09:17:24
言語 C
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 2,295 bytes
コンパイル時間 150 ms
コンパイル使用メモリ 31,744 KB
実行使用メモリ 26,516 KB
最終ジャッジ日時 2024-10-03 05:26:30
合計ジャッジ時間 46,704 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample WA * 1
other WA * 45
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>

int main () {
  int t = 0;
  
  int res = 0;
  
  long long ps[300001][10] = {};
  int pcnt[300001] = {};
  long long fr[300000][2] = {};
  
  for (int i = 2; i <= 300000; i++) {
    if (pcnt[i] <= 0) {
      ps[i][0] = (long long)i;
      pcnt[i] = 1;
      for (int p = 2; i*p <= 300000; p++) {
        ps[i*p][pcnt[i*p]] = i;
        pcnt[i*p]++;
      }
    }
  }
  
  res = scanf("%d", &t);
  while (t > 0) {
    long long n = 0LL;
    long long k = 0LL;
    long long rcnt = 1LL;
    res = scanf("%lld", &n);
    res = scanf("%lld", &k);
    rcnt = n*n;
    for (long long i = 1LL; i <= n; i += 1LL) {
      for (int j = 1; j < (1<<pcnt[(int)i]); j++) {
        int bcnt = 0;
        long long p = 1LL;
        for (int l = 0; l < pcnt[(int)i]; l++) {
          if ((j&(1<<l)) > 0) {
            p *= ps[(int)i][l];
            bcnt++;
          }
        }
        if (bcnt%2 == 1) {
          rcnt += n/p;
        } else {
          rcnt -= n/p;
        }
      }
    }
    if (rcnt < k) {
      printf("-1\n");
    } else {
      long long ans[2] = {};
      long long l = 0LL;
      long long r = (1LL<<19LL);
      long long div = 1LL;
      long long lcnt = 0LL;
      int frcnt = 0;
      while (div <= n) {
        long long nxt = l+r;
        long long cnt = 0LL;
        if (nxt%2LL == 1LL) {
          div *= 2LL;
          l *= 2LL;
          r *= 2LL;
        } else {
          nxt /= 2LL;
        }
        for (long long i = 1LL; i <= n; i += 1LL) {
          long long d = (nxt*i)/div;
          long long diff = 0LL;
          if (d > n) {
            d = n;
          }
          for (int j = 1; j < (1<<pcnt[(int)i]); j++) {
            int bcnt = 0;
            long long p = 1LL;
            for (int l = 0; l < pcnt[(int)i]; l++) {
              if ((j&(1<<l)) > 0) {
                p *= ps[(int)i][l];
                bcnt++;
              }
            }
            if (bcnt%2 == 1) {
              diff += d/p;
            } else {
              diff -= d/p;
            }
          }
          cnt += d-diff;
        }
        if (cnt < k) {
          l = nxt;
          lcnt = cnt;
        } else {
          r = nxt;
          rcnt = cnt;
        }
      }
      printf("%lld/%lld\n", ans[0], ans[1]);
    }
    t--;
  }
  
  return 0;
}
0