# |
제출 시각 |
아이디 |
문제 |
언어 |
결과 |
실행 시간 |
메모리 |
23222 |
2017-05-05T06:56:39 Z |
ainta |
hi (KRIII5_HI) |
C++14 |
|
500 ms |
7128 KB |
#include<cstdio>
#include<algorithm>
using namespace std;
int P[110], Mod = 1000000007, n, w[60], F[5110], InvF[5110], C[110][110];
int D[128][5010], Comb[5110][110], Num[60], TP[60][5010], CC[110][110];
int Pow(int a, int b){
int r = 1;
while(b){
if(b&1)r=1ll*r*a%Mod;
a=1ll*a*a%Mod;b>>=1;
}
return r;
}
void Do(int nd, int b, int e, int sum){
if(b==e){
Num[b] = nd;
return;
}
int m = (b+e)>>1, i, j, k, ss = sum;
for(i=0;i<=sum;i++)TP[0][i] = D[nd][i];
for(i=1;i<=m-b+1;i++)for(j=0;j<=5000;j++)TP[i][j] = 0;
for(i=0;i<=m-b;i++){
for(j=0;j<=ss;j++){
for(k=1;k<=w[i+b];k++){
TP[i+1][j+k] = (TP[i+1][j+k] + 1ll * Comb[j+k][k] * C[w[i+b]][k] % Mod * TP[i][j])%Mod;
}
}
ss += w[b+i];
}
for(i=0;i<=ss;i++)D[nd+nd+1][i] = TP[m-b+1][i];
Do(nd+nd+1, m+1, e, ss);
for(i=0;i<=sum;i++)TP[0][i] = D[nd][i];
for(i=1;i<=e-m;i++)for(j=0;j<=5000;j++)TP[i][j] = 0;
ss = sum;
for(i=1;i<=e-m;i++){
for(j=0;j<=ss;j++){
for(k=1;k<=w[i+m];k++){
TP[i][j+k] = (TP[i][j+k] + 1ll * Comb[j+k][k] * C[w[i+m]][k] % Mod * TP[i-1][j])%Mod;
}
}
ss += w[m+i];
}
for(i=0;i<=ss;i++)D[nd+nd][i] = TP[e-m][i];
Do(nd+nd, b, m, ss);
}
int main(){
int i, j, k;
for(i=1;i<=100;i++){
int s = 0;
for(j=1;j<i;j++){
s = (s + 1ll * P[j]*(i-j))%Mod;
}
P[i] = (i - s + Mod)%Mod;
}
F[0]=1;
for(i=1;i<=5000;i++){
F[i]=1ll*F[i-1]*i%Mod;
}
InvF[5000] = Pow(F[5000], Mod-2);
for(i=5000;i>=1;i--)InvF[i-1] = 1ll*InvF[i]*i%Mod;
for(i=0;i<=5000;i++){
for(j=0;j<=i&&j<=100;j++)Comb[i][j] = 1ll*F[i]*InvF[j]%Mod*InvF[i-j]%Mod;
}
C[0][0] = 1;
for(i=0;i<100;i++){
for(j=0;j<=100;j++){
for(k=0;k+j<=100;k++){
C[j+k][i+1] = (C[j+k][i+1] + 1ll * C[j][i] * P[k])%Mod;
}
}
}
scanf("%d",&n);
int sum = 0;
for(i=1;i<=n;i++){
scanf("%d",&w[i]);
sum += w[i];
}
D[1][0] = 1;
Do(1, 1, n, 0);
if(n==1){
printf("%d\n",sum);
return 0;
}
for(i=1;i<=100;i++){
for(j=1;j<i;j++){
for(k=1;k+j<=i;k++){
CC[i][j] = (CC[i][j] + 1ll*k*k*C[i-k][j])%Mod;
}
}
}
int res = 0;
for(i=1;i<=n;i++){
for(int ii = 1; ii <= sum;ii++){
if(!D[Num[i]][ii])continue;
int ss = w[i] * w[i];
for(j=1;j<w[i];j++){
ss = (ss + 1ll * Comb[ii + j - 2][j] * CC[w[i]][j])%Mod;
}
res = (res + 1ll * ss * D[Num[i]][ii])%Mod;
}
}
printf("%d\n",res);
}
Compilation message
HI.cpp: In function 'int main()':
HI.cpp:72:19: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
scanf("%d",&n);
^
HI.cpp:75:26: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
scanf("%d",&w[i]);
^
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
26 ms |
7128 KB |
Output is correct |
2 |
Correct |
26 ms |
7128 KB |
Output is correct |
3 |
Correct |
26 ms |
7128 KB |
Output is correct |
4 |
Correct |
26 ms |
7128 KB |
Output is correct |
5 |
Correct |
23 ms |
7128 KB |
Output is correct |
6 |
Correct |
23 ms |
7128 KB |
Output is correct |
7 |
Correct |
23 ms |
7128 KB |
Output is correct |
8 |
Correct |
23 ms |
7128 KB |
Output is correct |
9 |
Correct |
23 ms |
7128 KB |
Output is correct |
10 |
Correct |
23 ms |
7128 KB |
Output is correct |
11 |
Correct |
23 ms |
7128 KB |
Output is correct |
12 |
Correct |
16 ms |
7128 KB |
Output is correct |
13 |
Correct |
23 ms |
7128 KB |
Output is correct |
14 |
Correct |
23 ms |
7128 KB |
Output is correct |
15 |
Correct |
19 ms |
7128 KB |
Output is correct |
16 |
Correct |
19 ms |
7128 KB |
Output is correct |
17 |
Correct |
23 ms |
7128 KB |
Output is correct |
18 |
Correct |
26 ms |
7128 KB |
Output is correct |
19 |
Correct |
23 ms |
7128 KB |
Output is correct |
20 |
Correct |
23 ms |
7128 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
26 ms |
7128 KB |
Output is correct |
2 |
Correct |
33 ms |
7128 KB |
Output is correct |
3 |
Correct |
59 ms |
7128 KB |
Output is correct |
4 |
Correct |
143 ms |
7128 KB |
Output is correct |
5 |
Correct |
203 ms |
7128 KB |
Output is correct |
6 |
Correct |
326 ms |
7128 KB |
Output is correct |
7 |
Correct |
466 ms |
7128 KB |
Output is correct |
8 |
Execution timed out |
500 ms |
7128 KB |
Execution timed out |
9 |
Halted |
0 ms |
0 KB |
- |