結果

問題 No.472 平均順位
コンテスト
ユーザー kyawashell
提出日時 2018-03-28 12:10:55
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
MLE  
実行時間 -
コード長 1,290 bytes
コンパイル時間 1,943 ms
コンパイル使用メモリ 192,552 KB
最終ジャッジ日時 2025-01-05 09:39:02
ジャッジサーバーID
（参考情報） 		judge5 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2 MLE * 2
other AC * 10 TLE * 3 MLE * 3
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
typedef long long ll;
typedef unsigned long long ull;
using namespace std;
#define pb push_back
int 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][6000];
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];
		}
		else 
			ret = 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;
}
0