#include <bits/stdc++.h>
using namespace std;
#define rep(i,n) for (int (i)=(0);(i)<(int)(n);++(i))
using ll = long long;
//using P = pair<int, int>;
using namespace std;

template<class T> void vin(vector<T>& v, int n) {
    v.resize(n);
    for (int i = 0; i < n; ++i) {
        cin >> v[i];
    }
}

const int MAXN = 5010;
const int INF = 1e9;

int dp[2][3*MAXN+1];
int N, P;
// int a[MAXN], b[MAXN], c[MAXN];
//
// int solve(int idx, int p) {
//     int& res = dp[p];
//     if (res != INF) return res;
//     if (idx >= N) {
//         if (p > 0) res = 1e7;
//         else if (p == 0) res = 0;
//     }
//     else {
//         if (p >= 3) res = min(res, solve(idx+1, p-3)+1);
//         if (p >= 2) res = min(res, solve(idx+1, p-2)+c[idx]);
//         if (p >= 1) res = min(res, solve(idx+1, p-1)+b[idx]);
//         res = min(res, solve(idx+1, p) + a[idx]);
//     }
//     return res;
// }

int abc[MAXN][5];

int main() {
    cin >> N >> P;
    fill(dp[0], dp[2], INF);
    //rep(i, N) cin >> a[i] >> b[i] >> c[i];
    rep(i, N) cin >> abc[i][0] >> abc[i][1] >> abc[i][2], abc[i][3] = 1;
    dp[0][0] = 0;
    rep(i, N) {
        rep(j, P+1) {
            rep(k, 3+1) {
                if (j + k <= P) {
                    dp[1][j+k] = min(dp[0][j] + abc[i][k], dp[1][j+k]);
                }
            }
        }
        rep(j, P+1) {
            dp[0][j] = dp[1][j];
            dp[1][j] = INF;
        }
    }
    printf("%.10f\n", (double)dp[0][P] / (double)N);
    return 0;
}