結果

問題 No.28 末尾最適化
ユーザー LayCurseLayCurse
提出日時 2014-09-30 00:01:11
言語 C++11
(gcc 4.8.5)
結果
AC  
実行時間 239 ms / 5,000 ms
コード長 915 Byte
コンパイル時間 1,543 ms
使用メモリ 7,032 KB
最終ジャッジ日時 2021-01-23 10:48:20
合計ジャッジ時間 2,336 ms
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 2 ms
4,980 KB
testcase_01 AC 239 ms
4,980 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>
using namespace std;

#define REP(i,a,b) for(i=a;i<b;i++)
#define rep(i,n) REP(i,0,n)

#define ll long long
#define ull unsigned ll

#define INF 1000000000

int Q;
int seed, N, K, B;
int X[100000];
int Y[100000];

int p[20], pn[20], ps;

int main(){
  int i, j, k;
  int res;
  
  scanf("%d",&Q);
  while(Q--){
    scanf("%d%d%d%d",&seed,&N,&K,&B);
    N++;
    
    X[0] = seed;
    REP(i,1,N) X[i] = 1 + (((ll)X[i-1]*X[i-1] + (ll)X[i-1]*12345) % 100000009);

    ps = 0;
    for(i=2;i<=B;i++){
      if(B%i==0){
        p[ps] = i;
        pn[ps] = 0;
        while(B%i==0) pn[ps]++, B/=i;
        ps++;
      }
    }

    res = INF;
    rep(k,ps){
      rep(i,N){
        Y[i] = 0;
        j = X[i];
        while(j%p[k]==0) j/=p[k], Y[i]++;
      }
      sort(Y, Y+N);
      j = 0;
      rep(i,K) j += Y[i];
      res = min(res, j/pn[k]);
    }

    printf("%d\n",res);
  }

  return 0;
}
0