結果
問題 | No.472 平均順位 |
ユーザー |
![]() |
提出日時 | 2018-03-28 12:12:51 |
言語 | C++17(clang) (17.0.6 + boost 1.87.0) |
結果 |
MLE
|
実行時間 | - |
コード長 | 1,290 bytes |
コンパイル時間 | 2,049 ms |
コンパイル使用メモリ | 163,412 KB |
実行使用メモリ | 449,536 KB |
最終ジャッジ日時 | 2024-11-30 12:15:03 |
合計ジャッジ時間 | 14,100 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 MLE * 2 |
other | AC * 12 TLE * 3 MLE * 1 |
ソースコード
#include<bits/stdc++.h>typedef long long ll;typedef unsigned long long ull;using namespace std;#define pb push_backint dy[]={0, 0, 1, -1, 1, 1, -1, -1};int dx[]={1, -1, 0, 0, 1, -1, -1, 1};#define FOR(i,a,b) for (int i=(a);i<(b);i++)#define RFOR(i,a,b) for (int i=(b)-1;i>=(a);i--)#define REP(i,n) for (int i=0;i<(n);i++)#define RREP(i,n) for (int i=(n)-1;i>=0;i--)#define mp make_pair#define fi first#define sc second#define INF (1e9)int N,P;int a[4][5000];bool done[5001][15001];double memo[5001][15001];double solve(int i,int j) {if(done[i][j]) {return memo[i][j];}if(i == N - 1) {double ret;if(j < 4) {ret = a[j][i];}elseret = INF;done[i][j] = true;return memo[i][j] = ret;}double x = INF,y = INF,z = INF,w = INF;x = (solve(i + 1,j) * (N - i - 1) + a[0][i]) / (N - i);if(j > 0) {y = (solve(i + 1,j - 1) * (N - i - 1) + a[1][i]) / (N - i);}if(j > 1) {z = (solve(i + 1,j - 2) * (N - i - 1) + a[2][i]) / (N - i);}if(j > 2) {w = (solve(i + 1,j - 3) * (N - i - 1) + a[3][i]) / (N - i);}done[i][j] = true;return memo[i][j] = min(min(x,y),min(z,w));}int main(){cin >> N >> P;REP(i,N) {cin >> a[0][i] >> a[1][i] >> a[2][i];a[3][i] = 1;}printf("%0.10f\n",solve(0,P));return 0;}