結果

問題 No.174 カードゲーム（Hard）
ユーザー te-sh
提出日時 2018-07-13 16:27:08
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 353 ms / 2,000 ms
コード長 2,068 bytes
コンパイル時間 2,586 ms
コンパイル使用メモリ 159,976 KB
実行使用メモリ 35,200 KB
最終ジャッジ日時 2024-06-13 01:30:28
合計ジャッジ時間 6,148 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.algorithm, std.container, std.conv, std.math, std.range, std.typecons, std.stdio, std.string;

auto rdsp(){return readln.splitter;}
void pick(R,T)(ref R r,ref T t){t=r.front.to!T;r.popFront;}
void readV(T...)(ref T t){auto r=rdsp;foreach(ref v;t)pick(r,v);}
void readA(T)(size_t n,ref T[]t){t=new T[](n);auto r=rdsp;foreach(ref v;t)pick(r,v);}

void main()
{
  int n; real pa, pb; readV(n, pa, pb);
  int[] a; readA(n, a);
  int[] b; readA(n, b);

  a.sort();
  b.sort();

  auto calc(int[] a, real p)
  {
    auto q = new real[][](n, n);
    foreach (i; 0..n) q[i][] = 0;

    auto r = new real[](1<<n);
    r[] = 0; r[$-1] = 1;

    foreach_reverse (i; 1..1<<n) {
      auto c = i.popcnt, f = i.bsf;
      if (c == 1) {
        q[c-1][f] += r[i];
      } else {
        auto d = r[i]*p;
        q[c-1][f] += d;
        r[i.bitReset(f)] += d;
        foreach (j; f+1..n)
          if (i.bitTest(j)) {
            d = r[i]*(1-p)/(c-1);
            q[c-1][j] += d;
            r[i.bitReset(j)] += d;
          }
      }
    }

    return q;
  }

  auto r1 = calc(a, pa), r2 = calc(b, pb), e = 0.0L;

  foreach (i; 0..n)
    foreach (j; 0..n)
      foreach (k; 0..n)
        if (a[j] > b[k])
          e += r1[i][j]*r2[i][k]*(a[j]+b[k]);

  writefln("%.10f", e);
}

pragma(inline) {
  pure bool bitTest(T)(T n, size_t i) { return (n & (T(1) << i)) != 0; }
  pure T bitSet(T)(T n, size_t i) { return n | (T(1) << i); }
  pure T bitReset(T)(T n, size_t i) { return n & ~(T(1) << i); }
  pure T bitComp(T)(T n, size_t i) { return n ^ (T(1) << i); }

  pure T bitSet(T)(T n, size_t s, size_t e) { return n | ((T(1) << e) - 1) & ~((T(1) << s) - 1); }
  pure T bitReset(T)(T n, size_t s, size_t e) { return n & (~((T(1) << e) - 1) | ((T(1) << s) - 1)); }
  pure T bitComp(T)(T n, size_t s, size_t e) { return n ^ ((T(1) << e) - 1) & ~((T(1) << s) - 1); }

  import core.bitop;
  pure int bsf(T)(T n) { return core.bitop.bsf(ulong(n)); }
  pure int bsr(T)(T n) { return core.bitop.bsr(ulong(n)); }
  pure int popcnt(T)(T n) { return core.bitop.popcnt(ulong(n)); }
}
0