oj.uz

답안 #19657

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
19657 2016-02-25T04:15:22 Z kaTkaHr 카드 (kriii4_Z) C++14
100 / 100
 535 ms 142448 KB 
#include <stdio.h>
#include<vector>
#include <algorithm>
#include <map>

using namespace std;

typedef long long ll;

const int MX = 3005, MM = 1000000007;

template<typename T>
T pw(T A, ll B){
  T R = 1;
  while(B){
    if( B&1 ) R = R * A;
    A = A * A; B /= 2;
  }
  return R;
}

ll rv(ll A){ 
  ll R = 1, B = MM-2;
  while(B){
    if( B&1 ) R = R * A % MM;
    A = A * A % MM; B /= 2;
  }
  return R;
}
/*
struct frac{
  ll A, B;
		frac(ll A):A(A), B(1){}
  frac(ll a, ll b){
    A = (a%MM+MM) % MM;
    B = (b%MM+MM) % MM;
  }
  frac(){A = 0, B = 1;}
  frac operator+ (const frac &l)const{
    return frac((A * l.B + B * l.A) % MM, B * l.B % MM);
  }
  frac operator*(const frac &l)const{
    return frac(A*l.A % MM, B*l.B % MM);
  }
  frac operator/(const frac &l)const{
    return frac(A*l.B % MM, B*l.A % MM);
  }
  frac operator- (const frac &l)const{
    return frac((A*l.B - B*l.A%MM + MM) % MM, B*l.B % MM);
  }
  ll v(){ return A * rv(B) % MM; }
};// */
//*
struct frac{
  ll A;
		frac(ll A):A(A){}
  frac(ll a, ll b){
			a = (a%MM+MM)%MM;
			b = (b%MM+MM)%MM;
			A = rv(b) * a % MM;
  }
  frac(){A = 0;}
  frac operator+ (const frac &l)const{
			return l.A + A >= MM? l.A + A - MM: l.A + A;
  }
  frac operator*(const frac &l)const{
			return l.A * A % MM;
  }
  frac operator/(const frac &l)const{
			return A * rv(l.A) % MM;
  }
  frac operator- (const frac &l)const{
			return A >= l.A? A-l.A: A-l.A + MM;
  }
  ll v(){ return A; }
};// */

int D[MX];
frac T[MX][MX];
frac C[MX][MX];
frac A[11][MX];

int main()
{
	int N, L;
	scanf("%d%d", &N, &L);
	for(int i = 1; i <= N; i++){
		int A;
		scanf("%d", &A);
		D[A]++;
	}
	C[0][0] = 1;
	for(int i = 1; i <= 3000; i++){
		C[i][0] = 1;
		for(int j = 1; j <= i; j++){
			C[i][j] = C[i-1][j] + C[i-1][j-1];
		}
	}
	A[0][0] = 1;
	for(int i = 1; i <= L; i++) A[0][i] = A[0][i-1] * D[0];
	for(int k = 1; k <= 10; k++){
		T[0][0] = 1;
		for(int i = 1; i <= D[k]; i++){
			for(int j = 1; j <= L; j++){
				T[i][j] = T[i][j-1] * i;
				if( j >= k) T[i][j] = T[i][j] + T[i-1][j-k] * C[j-1][k-1];
			}
		}
		frac mul = 1;
		for(int i = 1; i <= D[k]; i++) mul = mul * i;
		for(int i = k*D[k]; i <= L; i++){
			for(int j = k*D[k]; j <= i; j++){
				A[k][i] = A[k][i] + A[k-1][i-j] * T[D[k]][j] * C[i][j];
			}
		}
		for(int i = 0; i <= L; i++) A[k][i] = A[k][i] * mul;
	}

	frac ans = A[10][L];
	for(int i = 1; i <= L; i++) ans = ans / N;
	printf("%lld\n", ans.v());
}

Subtask #1
13.0 / 13.0

# 결과 실행 시간 메모리 Grader output
1 Correct 431 ms 142448 KB Output is correct
2 Correct 455 ms 142448 KB Output is correct
3 Correct 385 ms 142448 KB Output is correct
4 Correct 464 ms 142448 KB Output is correct
5 Correct 389 ms 142448 KB Output is correct
6 Correct 339 ms 142448 KB Output is correct
7 Correct 414 ms 142448 KB Output is correct
8 Correct 364 ms 142448 KB Output is correct
9 Correct 486 ms 142448 KB Output is correct
10 Correct 368 ms 142448 KB Output is correct
11 Correct 426 ms 142448 KB Output is correct
12 Correct 379 ms 142448 KB Output is correct
13 Correct 373 ms 142448 KB Output is correct
14 Correct 442 ms 142448 KB Output is correct
15 Correct 446 ms 142448 KB Output is correct
16 Correct 330 ms 142448 KB Output is correct
17 Correct 455 ms 142448 KB Output is correct
18 Correct 362 ms 142448 KB Output is correct
19 Correct 435 ms 142448 KB Output is correct
20 Correct 390 ms 142448 KB Output is correct
21 Correct 529 ms 142448 KB Output is correct
22 Correct 535 ms 142448 KB Output is correct
23 Correct 524 ms 142448 KB Output is correct
24 Correct 532 ms 142448 KB Output is correct
25 Correct 534 ms 142448 KB Output is correct

Subtask #2
87.0 / 87.0

# 결과 실행 시간 메모리 Grader output
1 Correct 394 ms 142448 KB Output is correct
2 Correct 423 ms 142448 KB Output is correct
3 Correct 400 ms 142448 KB Output is correct
4 Correct 385 ms 142448 KB Output is correct
5 Correct 427 ms 142448 KB Output is correct
6 Correct 390 ms 142448 KB Output is correct
7 Correct 366 ms 142448 KB Output is correct
8 Correct 430 ms 142448 KB Output is correct
9 Correct 384 ms 142448 KB Output is correct
10 Correct 375 ms 142448 KB Output is correct
11 Correct 386 ms 142448 KB Output is correct
12 Correct 409 ms 142448 KB Output is correct
13 Correct 365 ms 142448 KB Output is correct
14 Correct 419 ms 142448 KB Output is correct
15 Correct 384 ms 142448 KB Output is correct
16 Correct 443 ms 142448 KB Output is correct
17 Correct 448 ms 142448 KB Output is correct
18 Correct 398 ms 142448 KB Output is correct
19 Correct 382 ms 142448 KB Output is correct
20 Correct 405 ms 142448 KB Output is correct