結果

問題 No.463 魔法使いのすごろく🎲
ユーザー char134217728char134217728
提出日時 2017-09-03 17:15:01
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 2,875 bytes
コンパイル時間 1,668 ms
コンパイル使用メモリ 162,780 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-04-24 12:15:12
合計ジャッジ時間 2,719 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 WA -
testcase_03 AC 2 ms
5,376 KB
testcase_04 WA -
testcase_05 WA -
testcase_06 AC 1 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 AC 1 ms
5,376 KB
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 AC 1 ms
5,376 KB
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 AC 2 ms
5,376 KB
testcase_23 WA -
testcase_24 AC 2 ms
5,376 KB
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
testcase_32 WA -
testcase_33 WA -
testcase_34 WA -
testcase_35 AC 1 ms
5,376 KB
testcase_36 AC 2 ms
5,376 KB
testcase_37 AC 1 ms
5,376 KB
testcase_38 AC 1 ms
5,376 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:81:1: warning: ISO C++ forbids declaration of ‘main’ with no type [-Wreturn-type]
   81 | main(){
      | ^~~~

ソースコード

diff # 

#include <bits/stdc++.h>
#define FOR(i,a,b) for (int i=(a);i<(b);i++)
#define FORR(i,a,b) for (int i=(a);i>=(b);i--)
#define pb push_back
#define pcnt __builtin_popcount
#define show(x) cout<<#x<<" = "<<x<<endl;
#define maxs(x,y) x = max(x,y)
#define mins(x,y) x = min(x,y)
#define fi first
#define se second
#define rng(a) a.begin(),a.end()
#define each(it,c) for(__typeof((c).begin()) it=(c).begin();it!=(c).end();it++)
#define sz(x) (int)(x).size()

using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<pii> vpii;
typedef set<int> si;
typedef pair<ll,ll> pll;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<pll> vpll;
typedef set<ll> sl;
template<typename T>string join(vector<T>&v)
{stringstream s;FOR(i,0,sz(v))s<<' '<<v[i];return s.str().substr(1);}
ll gcd(ll a,ll b){if(a>b)swap(a,b);for(;a>0;b%=a,swap(a,b));return b;}
int modpow(ll a,ll n,int m){if(a==0)return a;ll p=1;for(;n>0;n/=2,a=a*a%m)if(n&1)p=p*a%m;return(int)p;}
void dout(double d){printf("%.12f\n",d);}

const int iinf = 1e9;
const ll linf = 1e18;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-10;


const int N=103;
const int GN = N;

struct Gauss{
  int solve(int n, double o[GN][GN], double u[GN], double r[GN]) {
    int q;
    double m[GN][GN],v[GN], p;
    memmove(m,o,sizeof(m));
    memmove(v,u,sizeof(v));
    FOR(i, 0, n) {
      p = 0;
      q = 0;
      FOR(j, i, n) if(p < abs(m[i][j])){
        p = abs(m[i][j]);
        q = j;
      }
      if(p < eps) return -1;
      swap(v[i],v[q]);
      FOR(k, i, n) swap(m[i][k], m[q][k]);
      v[i]/=m[i][i];
      FORR(k, n-1, i) m[i][k]/=m[i][i];
      FOR(j, i+1, n){
        v[j]-=v[i]*m[j][i];
        FORR(k, n-1, i) m[j][k]-=m[i][k]*m[j][i];
      }
    }
    FORR(i, n-1, 0) {
      FORR(j, n-1, i+1){
        v[i]-=m[i][j]*v[j];
        m[i][j]=0;
      }
      r[i]=v[i];
    }
    return 0;
  }
};

double M[GN][GN], V[GN], R[GN], DP[2][GN];
int n, m, c[N];
Gauss g;

main(){
  cin.tie(0);
  ios::sync_with_stdio(false);
  cin >> n >> m;
  FOR(i, 2, n) cin >> c[i];
  double rm = 1./m;
  int o = (m+1)/2, _j;
  FOR(i, 1, o+1){
    FOR(j, 1, i){
      M[i-1][j-1] += rm;
      V[i-1] += rm * c[n-j];
    }
    FOR(j, 1, m-i+1){
      _j = j > o ? m+1-j : j;
      M[i-1][_j-1] += rm;
      V[i-1] += rm * c[n-j];
    }
  }
  g.solve(o, M, V, R);
  FOR(i, 0, o) DP[0][n-1-i] = R[i];


  FORR(i, n-o-1, 1){
    FOR(j, 1, m+1){
      _j = i + j;
      if(_j > n) _j = 2 * n - _j;
      DP[0][i] += (DP[0][_j] + c[_j]) * rm;
    }
  }


  FORR(i, n-m-1, 1){
    FOR(j, 1, m+1){
      _j = i + j;
      if(_j > n) _j = 2 * n - _j;
      DP[1][i] += (DP[1][_j] + c[_j]) * rm;
    }
    FOR(j, 1, m){
      _j = i + j;
      if(_j > n) _j = 2 * n - _j;
      mins(DP[1][i], DP[0][_j] + c[_j]);
    }
  }
  dout(DP[1][1]);
}
0