結果

問題 No.463 魔法使いのすごろく🎲
コンテスト
ユーザー char134217728
提出日時 2017-09-03 17:15:01
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 2,875 bytes
コンパイル時間 1,962 ms
コンパイル使用メモリ 161,760 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-11-06 20:05:08
合計ジャッジ時間 2,690 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 8 WA * 28
権限があれば一括ダウンロードができます
コンパイルメッセージ
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