#include #include #include using namespace atcoder; using mint = modint; using namespace std; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf 1000000001 int N,L; //long long dp[51][18002][51]; long long rdp[51][18002][51]; /* void get(vector t){ dp[0][0][0] = 1LL; rep(i,N){ rep(j,L){ rep(k,N+1){ if(dp[i][j][k]==0)continue; dp[i+1][j][k] += dp[i][j][k]; if(j+t[i] t){ rdp[0][0][0] = 1LL; rep(i,N){ rep(j,L){ rep(k,N+1){ if(rdp[i][j][k]==0)continue; rdp[i+1][j][k] += rdp[i][j][k]; if(j+t[i]>N>>L; L *= 60; vector t(N); long long sum = 0; rep(i,N){ int m,s; scanf("%d:%d",&m,&s); t[i] = m*60+s; sum += t[i]; } vector f(N+5,1.0); for(int i=1;i(N+1,0)); dp[0][0] = 1; rep(i,N+1){ rep(j,L-1){ rep(k,N+1){ rdp[i][j+1][k] += rdp[i][j][k]; } } } rep(i,N){ rep(j,L){ rep(k,N+1){ if(dp[j][k]==0)continue; rep(kk,N+1){ double temp = dp[j][k] * rdp[N-1-i][L-1-j][kk]; temp *= f[k+kk]; temp *= f[(N-1)-k-kk]; ans += temp; } } } vector ndp(L,vector(N+1,0)); rep(j,L){ rep(k,N+1){ if(dp[j][k]==0)continue; ndp[j][k] += dp[j][k]; if(j+t[i]