#362710 (D) No.856 増える演算

提出ソース

結果

問題	No.856 増える演算
ユーザー	👑 hos.lyric
提出日時	2019-07-26 22:30:58
言語	D (dmd 2.106.1)
結果	WA
実行時間	-
コード長	4,420 bytes
コンパイル時間	2,438 ms
コンパイル使用メモリ	176,340 KB
実行使用メモリ	34,336 KB
最終ジャッジ日時	2024-06-22 02:02:46
合計ジャッジ時間	21,596 ms
ジャッジサーバーID （参考情報）	judge2 / judge4

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	190 ms 28,260 KB
testcase_01	AC	186 ms 27,920 KB
testcase_02	AC	186 ms 30,652 KB
testcase_03	AC	185 ms 27,848 KB
testcase_04	AC	185 ms 28,896 KB
testcase_05	AC	192 ms 29,552 KB
testcase_06	AC	189 ms 30,248 KB
testcase_07	AC	187 ms 29,568 KB
testcase_08	AC	192 ms 28,384 KB
testcase_09	AC	189 ms 30,692 KB
testcase_10	AC	189 ms 28,124 KB
testcase_11	AC	190 ms 30,612 KB
testcase_12	AC	187 ms 29,812 KB
testcase_13	AC	191 ms 30,504 KB
testcase_14	AC	185 ms 30,320 KB
testcase_15	AC	194 ms 30,576 KB
testcase_16	AC	189 ms 28,104 KB
testcase_17	AC	188 ms 28,584 KB
testcase_18	AC	191 ms 29,224 KB
testcase_19	AC	187 ms 28,104 KB
testcase_20	AC	189 ms 31,192 KB
testcase_21	AC	198 ms 27,688 KB
testcase_22	AC	192 ms 30,612 KB
testcase_23	AC	193 ms 32,916 KB
testcase_24	AC	187 ms 29,012 KB
testcase_25	AC	187 ms 29,748 KB
testcase_26	AC	188 ms 31,952 KB
testcase_27	AC	188 ms 29,252 KB
testcase_28	AC	191 ms 30,000 KB
testcase_29	AC	182 ms 27,948 KB
testcase_30	AC	185 ms 30,144 KB
testcase_31	AC	187 ms 29,820 KB
testcase_32	AC	184 ms 29,840 KB
testcase_33	AC	189 ms 31,692 KB
testcase_34	AC	199 ms 28,420 KB
testcase_35	AC	211 ms 31,496 KB
testcase_36	AC	203 ms 29,004 KB
testcase_37	AC	204 ms 28,608 KB
testcase_38	AC	200 ms 29,084 KB
testcase_39	AC	205 ms 29,172 KB
testcase_40	AC	204 ms 30,888 KB
testcase_41	AC	209 ms 31,596 KB
testcase_42	AC	214 ms 30,368 KB
testcase_43	AC	201 ms 29,276 KB
testcase_44	AC	201 ms 29,152 KB
testcase_45	AC	182 ms 30,548 KB
testcase_46	AC	185 ms 28,272 KB
testcase_47	AC	191 ms 28,920 KB
testcase_48	AC	197 ms 28,976 KB
testcase_49	AC	194 ms 30,448 KB
testcase_50	AC	193 ms 31,864 KB
testcase_51	AC	193 ms 29,620 KB
testcase_52	AC	197 ms 30,360 KB
testcase_53	AC	227 ms 28,492 KB
testcase_54	AC	214 ms 29,720 KB
testcase_55	AC	229 ms 29,988 KB
testcase_56	AC	213 ms 29,340 KB
testcase_57	AC	235 ms 33,208 KB
testcase_58	AC	225 ms 30,084 KB
testcase_59	AC	259 ms 31,544 KB
testcase_60	AC	220 ms 30,484 KB
testcase_61	AC	241 ms 31,456 KB
testcase_62	AC	238 ms 32,508 KB
testcase_63	AC	187 ms 29,504 KB
testcase_64	AC	233 ms 28,872 KB
testcase_65	AC	208 ms 28,616 KB
testcase_66	AC	220 ms 28,616 KB
testcase_67	AC	225 ms 30,152 KB
testcase_68	AC	249 ms 30,300 KB
testcase_69	AC	244 ms 30,908 KB
testcase_70	AC	251 ms 31,688 KB
testcase_71	AC	250 ms 32,744 KB
testcase_72	AC	237 ms 29,840 KB
testcase_73	AC	255 ms 32,028 KB
testcase_74	AC	267 ms 34,336 KB
testcase_75	AC	254 ms 33,116 KB
testcase_76	AC	277 ms 31,316 KB
testcase_77	AC	260 ms 30,856 KB
testcase_78	AC	258 ms 33,732 KB
testcase_79	AC	274 ms 32,072 KB
testcase_80	AC	274 ms 33,196 KB
testcase_81	AC	258 ms 31,716 KB
testcase_82	WA	-

権限があれば一括ダウンロードができます

コンパイルメッセージ

Complex!real[]
Main.d(145): Deprecation: function `std.math.exponential.log` is deprecated - `std.math.exponential.log` called with argument types `(long)` matches both `log(real)`, `log(double)`, and `log(float)`. Cast argument to floating point type instead.
Main.d(145): Deprecation: function `std.math.exponential.log` is deprecated - `std.math.exponential.log` called with argument types `(long)` matches both `log(real)`, `log(double)`, and `log(float)`. Cast argument to floating point type instead.

ソースコード

raw source code

import std.conv, std.functional, std.range, std.stdio, std.string;
import std.algorithm, std.array, std.bigint, std.complex, std.container, std.math, std.numeric, std.regex, std.typecons;
import core.bitop;

class EOFException : Throwable { this() { super("EOF"); } }
string[] tokens;
string readToken() { for (; tokens.empty; ) { if (stdin.eof) { throw new EOFException; } tokens = readln.split; } auto token = tokens.front; tokens.popFront; return token; }
int readInt() { return readToken.to!int; }
long readLong() { return readToken.to!long; }
real readReal() { return readToken.to!real; }

bool chmin(T)(ref T t, in T f) { if (t > f) { t = f; return true; } else { return false; } }
bool chmax(T)(ref T t, in T f) { if (t < f) { t = f; return true; } else { return false; } }

int binarySearch(alias pred, T)(in T[] as) { int lo = -1, hi = cast(int)(as.length); for (; lo + 1 < hi; ) { const mid = (lo + hi) >> 1; (unaryFun!pred(as[mid]) ? hi : lo) = mid; } return hi; }
int lowerBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a >= val)); }
int upperBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a > val)); }


struct ModInt(long M) {
  long x;
  this(in ModInt a) { x = a.x; }
  this(in long a) { x = a % M; if (x < 0) x += M; }
  ref ModInt opAssign(in long a) { return this = ModInt(a); }
  ref ModInt opOpAssign(string op)(in ModInt a) {
    static if (op == "+") { x += a.x; if (x >= M) x -= M; }
    else static if (op == "-") { x -= a.x; if (x < 0) x += M; }
    else static if (op == "*") { x *= a.x; x %= M; }
    else static assert(false);
    return this;
  }
  ref ModInt opOpAssign(string op)(in long a) { return mixin("this " ~ op ~ "= ModInt(a)"); }
  ModInt opUnary(string op)() const if (op == "-") { return ModInt(-x); }
  ModInt opBinary(string op, T)(in T a) const { return mixin("ModInt(this) " ~ op ~ "= a"); }
  ModInt opBinaryRight(string op)(in long a) const { return mixin("ModInt(a) " ~ op ~ "= this"); }
  string toString() const { return x.to!string; }
}


// a^-1 (mod m)
long modInv(long a, long m)
in {
  assert(m > 0, "modInv: m > 0 must hold");
}
do {
  long b = m, x = 1, y = 0, t;
  for (; ; ) {
    t = a / b;
    a -= t * b;
    if (a == 0) {
      assert(b == 1 || b == -1, "modInv: gcd(a, m) != 1");
      if (b == -1) {
        y = -y;
      }
      return (y < 0) ? (y + m) : y;
    }
    x -= t * y;
    t = b / a;
    b -= t * a;
    if (b == 0) {
      assert(a == 1 || a == -1, "modInv: gcd(a, m) != 1");
      if (a == -1) {
        x = -x;
      }
      return (x < 0) ? (x + m) : x;
    }
    y -= t * x;
  }
}


enum long MO = 1000000007;
alias Mint = ModInt!MO;

Mint power(Mint a, long e) {
  Mint x = a, y = 1;
  for (; e; e >>= 1) {
    if (e & 1) y *= x;
    x *= x;
  }
  return y;
}

enum L = 2^^18;

int N;
long[] A;

void main() {
  auto fft = new Fft(L);
  
  try {
    for (; ; ) {
      N = readInt();
      A = new long[N];
      foreach (i; 0 .. N) {
        A[i] = readInt();
      }
foreach(i;0..N)if(!(1<=A[i]&&A[i]<=10^^5))for(;;){}
      
      auto f = new real[L];
      f[] = 0.0;
      foreach (i; 0 .. N) {
        f[cast(int)(A[i])] += 1;
      }
      auto ff = fft.fft!real(f);
      foreach (l; 0 .. L) {
        ff[l] = ff[l] * ff[l];
      }
      auto g = fft.inverseFft!real(ff);
pragma(msg,typeof(g));
      
      auto cnt = new int[L];
      foreach (l; 0 .. L) {
        cnt[l] = cast(int)(round(g[l].re));
      }
      foreach (i; 0 .. N) {
        --cnt[cast(int)(2 * A[i])];
      }
// foreach(l;0..L){
 // if(!(cnt[l]>=0))for(;;){}
 // if(!(cnt[l]%2==0))for(;;){}
// }
      cnt[] /= 2;
      debug {
        writeln("cnt = ", cnt[0 .. 20]);
      }
      
      Mint ans = 1;
      foreach (l; 0 .. L) {
        if (cnt[l] > 0) {
          ans *= power(Mint(l), cnt[l]);
        }
      }
      long ASum;
      foreach_reverse (i; 0 .. N) {
        ans *= power(Mint(A[i]), ASum);
        ASum += A[i];
      }
      
      long AMin = long.max;
      auto opt = tuple(real.infinity, 0L, 0L);
      foreach_reverse (i; 0 .. N) {
        if (i < N - 1) {
          chmin(opt, tuple(log(A[i] + AMin) + log(A[i]) * AMin, A[i], AMin));
        }
        chmin(AMin, A[i]);
      }
      debug {
        writeln("opt = ", opt);
      }
      Mint dnm = (opt[1] + opt[2]) * power(Mint(opt[1]), opt[2]);
      ans *= modInv(dnm.x, MO);
      
      writeln(ans);
    }
  } catch (EOFException e) {
  }
}

yukicoder

結果

テストケース

コンパイルメッセージ

ソースコード