Submission #897716

#	Submission time	Handle	Problem	Language	Result	Execution time	Memory
897716	2024-01-03T15:32:43 Z	denniskim	NoM (RMI21_nom)	C++17	100 / 100	155 ms	31992 KB

Main

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef __int128 lll;
typedef long double ld;
typedef pair<ll, ll> pll;
typedef pair<ld, ld> pld;
#define MAX 9223372036854775807LL
#define MIN -9223372036854775807LL
#define INF 0x3f3f3f3f3f3f3f3f
#define fi first
#define se second
#define fastio ios_base::sync_with_stdio(false); cin.tie(NULL); cout.tie(NULL); cout << fixed; cout.precision(10);
#define sp << " "
#define en << "\n"
#define compress(v) sort(v.begin(), v.end()), v.erase(unique(v.begin(), v.end()), v.end())

ll n, m;
ll ss = 1000000007;
ll ans;
ll cou[2010];
ll dp[2010][2010];

struct NCK
{
	ll fac[4010] = {0, }, inv[4010] = {0, };
	
	ll power(ll xx, ll yy)
	{
	    ll ret = 1;
	     
	    while(yy > 0)
	    {
	        if(yy % 2)
	            ret = (ret * xx) % ss;
	        xx = (xx * xx) % ss;
	        yy /= 2;
	    }
	     
	    return ret;
	}
	 
	ll C(ll xx, ll yy)
	{
		if(xx < yy)
			return 0;
		
	    if(xx == yy || !yy)
	        return 1;
	    
	    long long r = (fac[xx] * inv[xx - yy]) % ss;
	    r = (r * inv[yy]) % ss;
	    return r;
	}
	
	void init(void)
	{
		fac[0] = 1;
	    fac[1] = 1;
	     
	    for(ll i = 2 ; i <= 4000 ; i++)
	        fac[i] = (fac[i - 1] * i) % ss;
	     
	    inv[4000] = power(fac[4000], ss - 2);
	    
	    for(ll i = 3999 ; i >= 0 ; i--)
	        inv[i] = (inv[i + 1] * (i + 1)) % ss;
	}
}nck;

int main(void)
{
	fastio
	
	cin >> n >> m;
	
	nck.init();
	
	ans = nck.fac[n * 2];
	
	for(ll i = 1 ; i <= n * 2 ; i++)
	{
		ll gap = i % m;
		
		if(gap == 0)
			gap = m;
		
		cou[gap]++;
	}
	
	for(ll i = 0 ; i <= m ; i++)
		dp[0][i] = 1;
	
	for(ll i = 1 ; i <= n ; i++)
	{
		for(ll j = 1 ; j <= m ; j++)
		{
			ll X = cou[j];
			
			for(ll k = 1 ; k * 2 <= X && k <= i ; k++)
				dp[i][j] = (dp[i][j] + dp[i - k][j - 1] * nck.C(X, k * 2) % ss * nck.fac[k * 2] % ss * nck.C(i, k) % ss) % ss;
		}
		
		for(ll j = 1 ; j <= m ; j++)
			dp[i][j] = (dp[i][j] + dp[i][j - 1]) % ss;
	}
	
	for(ll i = 1 ; i <= n ; i++)
	{
		ll gap = dp[i][m] * nck.C(n, i) % ss * nck.fac[n * 2 - i * 2] % ss;
		
		if(i & 1)
			ans = (ans - gap) % ss;
		else
			ans = (ans + gap) % ss;
	}
	
	ans = (ans + ss) % ss;
	
	cout << ans;
	return 0;
}

Subtask #1
9.0 / 9.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	1 ms	344 KB	Output is correct
2	Correct	0 ms	348 KB	Output is correct
3	Correct	0 ms	348 KB	Output is correct
4	Correct	0 ms	348 KB	Output is correct
5	Correct	0 ms	348 KB	Output is correct
6	Correct	0 ms	424 KB	Output is correct
7	Correct	0 ms	600 KB	Output is correct
8	Correct	0 ms	348 KB	Output is correct
9	Correct	0 ms	348 KB	Output is correct
10	Correct	0 ms	348 KB	Output is correct

Subtask #2
12.0 / 12.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	1 ms	344 KB	Output is correct
2	Correct	0 ms	348 KB	Output is correct
3	Correct	0 ms	348 KB	Output is correct
4	Correct	0 ms	348 KB	Output is correct
5	Correct	0 ms	348 KB	Output is correct
6	Correct	0 ms	424 KB	Output is correct
7	Correct	0 ms	600 KB	Output is correct
8	Correct	0 ms	348 KB	Output is correct
9	Correct	0 ms	348 KB	Output is correct
10	Correct	0 ms	348 KB	Output is correct
11	Correct	1 ms	860 KB	Output is correct
12	Correct	1 ms	604 KB	Output is correct
13	Correct	1 ms	604 KB	Output is correct
14	Correct	1 ms	604 KB	Output is correct
15	Correct	1 ms	860 KB	Output is correct
16	Correct	1 ms	860 KB	Output is correct
17	Correct	1 ms	604 KB	Output is correct
18	Correct	1 ms	604 KB	Output is correct
19	Correct	1 ms	860 KB	Output is correct

Subtask #3
13.0 / 13.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	1 ms	344 KB	Output is correct
2	Correct	0 ms	348 KB	Output is correct
3	Correct	0 ms	348 KB	Output is correct
4	Correct	0 ms	348 KB	Output is correct
5	Correct	0 ms	348 KB	Output is correct
6	Correct	0 ms	424 KB	Output is correct
7	Correct	0 ms	600 KB	Output is correct
8	Correct	0 ms	348 KB	Output is correct
9	Correct	0 ms	348 KB	Output is correct
10	Correct	0 ms	348 KB	Output is correct
11	Correct	1 ms	860 KB	Output is correct
12	Correct	1 ms	604 KB	Output is correct
13	Correct	1 ms	604 KB	Output is correct
14	Correct	1 ms	604 KB	Output is correct
15	Correct	1 ms	860 KB	Output is correct
16	Correct	1 ms	860 KB	Output is correct
17	Correct	1 ms	604 KB	Output is correct
18	Correct	1 ms	604 KB	Output is correct
19	Correct	1 ms	860 KB	Output is correct
20	Correct	2 ms	1112 KB	Output is correct
21	Correct	2 ms	1124 KB	Output is correct
22	Correct	4 ms	2140 KB	Output is correct
23	Correct	4 ms	1528 KB	Output is correct
24	Correct	4 ms	2136 KB	Output is correct
25	Correct	3 ms	1880 KB	Output is correct
26	Correct	4 ms	1368 KB	Output is correct
27	Correct	4 ms	2396 KB	Output is correct

Subtask #4
18.0 / 18.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	1 ms	344 KB	Output is correct
2	Correct	0 ms	348 KB	Output is correct
3	Correct	0 ms	348 KB	Output is correct
4	Correct	0 ms	348 KB	Output is correct
5	Correct	0 ms	348 KB	Output is correct
6	Correct	0 ms	424 KB	Output is correct
7	Correct	0 ms	600 KB	Output is correct
8	Correct	0 ms	348 KB	Output is correct
9	Correct	0 ms	348 KB	Output is correct
10	Correct	0 ms	348 KB	Output is correct
11	Correct	1 ms	860 KB	Output is correct
12	Correct	1 ms	604 KB	Output is correct
13	Correct	1 ms	604 KB	Output is correct
14	Correct	1 ms	604 KB	Output is correct
15	Correct	1 ms	860 KB	Output is correct
16	Correct	1 ms	860 KB	Output is correct
17	Correct	1 ms	604 KB	Output is correct
18	Correct	1 ms	604 KB	Output is correct
19	Correct	1 ms	860 KB	Output is correct
20	Correct	2 ms	1112 KB	Output is correct
21	Correct	2 ms	1124 KB	Output is correct
22	Correct	4 ms	2140 KB	Output is correct
23	Correct	4 ms	1528 KB	Output is correct
24	Correct	4 ms	2136 KB	Output is correct
25	Correct	3 ms	1880 KB	Output is correct
26	Correct	4 ms	1368 KB	Output is correct
27	Correct	4 ms	2396 KB	Output is correct
28	Correct	23 ms	2904 KB	Output is correct
29	Correct	20 ms	4700 KB	Output is correct
30	Correct	17 ms	5888 KB	Output is correct
31	Correct	23 ms	3136 KB	Output is correct
32	Correct	22 ms	4960 KB	Output is correct
33	Correct	20 ms	6932 KB	Output is correct
34	Correct	33 ms	10252 KB	Output is correct

Subtask #5
48.0 / 48.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	1 ms	344 KB	Output is correct
2	Correct	0 ms	348 KB	Output is correct
3	Correct	0 ms	348 KB	Output is correct
4	Correct	0 ms	348 KB	Output is correct
5	Correct	0 ms	348 KB	Output is correct
6	Correct	0 ms	424 KB	Output is correct
7	Correct	0 ms	600 KB	Output is correct
8	Correct	0 ms	348 KB	Output is correct
9	Correct	0 ms	348 KB	Output is correct
10	Correct	0 ms	348 KB	Output is correct
11	Correct	1 ms	860 KB	Output is correct
12	Correct	1 ms	604 KB	Output is correct
13	Correct	1 ms	604 KB	Output is correct
14	Correct	1 ms	604 KB	Output is correct
15	Correct	1 ms	860 KB	Output is correct
16	Correct	1 ms	860 KB	Output is correct
17	Correct	1 ms	604 KB	Output is correct
18	Correct	1 ms	604 KB	Output is correct
19	Correct	1 ms	860 KB	Output is correct
20	Correct	2 ms	1112 KB	Output is correct
21	Correct	2 ms	1124 KB	Output is correct
22	Correct	4 ms	2140 KB	Output is correct
23	Correct	4 ms	1528 KB	Output is correct
24	Correct	4 ms	2136 KB	Output is correct
25	Correct	3 ms	1880 KB	Output is correct
26	Correct	4 ms	1368 KB	Output is correct
27	Correct	4 ms	2396 KB	Output is correct
28	Correct	23 ms	2904 KB	Output is correct
29	Correct	20 ms	4700 KB	Output is correct
30	Correct	17 ms	5888 KB	Output is correct
31	Correct	23 ms	3136 KB	Output is correct
32	Correct	22 ms	4960 KB	Output is correct
33	Correct	20 ms	6932 KB	Output is correct
34	Correct	33 ms	10252 KB	Output is correct
35	Correct	147 ms	7252 KB	Output is correct
36	Correct	121 ms	28532 KB	Output is correct
37	Correct	145 ms	7128 KB	Output is correct
38	Correct	114 ms	17492 KB	Output is correct
39	Correct	98 ms	25272 KB	Output is correct
40	Correct	155 ms	31992 KB	Output is correct

Submission #897716

Compilation message

Subtask #1 9.0 / 9.0

Subtask #2 12.0 / 12.0

Subtask #3 13.0 / 13.0

Subtask #4 18.0 / 18.0

Subtask #5 48.0 / 48.0

Subtask #1
9.0 / 9.0

Subtask #2
12.0 / 12.0

Subtask #3
13.0 / 13.0

Subtask #4
18.0 / 18.0

Subtask #5
48.0 / 48.0