結果

問題 No.2846 Birthday Cake
ユーザー 沙耶花沙耶花
提出日時 2024-08-23 21:59:03
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 1,380 bytes
コンパイル時間 4,766 ms
コンパイル使用メモリ 272,752 KB
実行使用メモリ 132,224 KB
最終ジャッジ日時 2024-08-23 21:59:38
合計ジャッジ時間 29,075 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 1,255 ms
103,424 KB
testcase_03 AC 1 ms
6,940 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 1,467 ms
115,968 KB
testcase_06 AC 1,583 ms
124,288 KB
testcase_07 AC 1,743 ms
129,024 KB
testcase_08 AC 1,747 ms
131,200 KB
testcase_09 AC 1,819 ms
131,968 KB
testcase_10 AC 1,929 ms
132,096 KB
testcase_11 TLE -
testcase_12 AC 1,046 ms
102,272 KB
testcase_13 AC 1,370 ms
126,976 KB
testcase_14 AC 3 ms
6,940 KB
testcase_15 AC 20 ms
6,944 KB
testcase_16 AC 86 ms
14,720 KB
testcase_17 AC 1,582 ms
129,536 KB
testcase_18 TLE -
testcase_19 AC 2 ms
6,940 KB
testcase_20 AC 3 ms
6,944 KB
testcase_21 AC 2 ms
6,944 KB
testcase_22 AC 12 ms
6,940 KB
testcase_23 AC 9 ms
6,940 KB
testcase_24 AC 1,419 ms
129,408 KB
testcase_25 AC 2 ms
6,944 KB
testcase_26 AC 6 ms
6,940 KB
testcase_27 AC 8 ms
6,940 KB
testcase_28 AC 2 ms
6,940 KB
testcase_29 AC 20 ms
6,940 KB
testcase_30 AC 2 ms
6,944 KB
testcase_31 AC 4 ms
6,940 KB
testcase_32 AC 283 ms
34,688 KB
testcase_33 AC 1,483 ms
129,408 KB
testcase_34 AC 74 ms
14,464 KB
testcase_35 AC 81 ms
14,848 KB
testcase_36 AC 399 ms
40,960 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using mint = modint998244353;
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf32 1000000001
#define Inf64 1000000000000000001

long long N,K;

long long En = 1;
vector<long long> la,lb;
long long L;
map<long long,__int128> mp[25];
map<long long,__int128> mp2[25];
void dfs(int ci,int ck,long long cx,__int128 cv){
	if(cx > L)return;

	if(ck>=la.size()){
		__int128 t = 1;
		rep(i,K){
			t *= i+1;
		}
		mp[ci][cx] += t/cv;
		return;
	}
	__int128 t = 1LL;
	rep(i,K-ci+1){
		if(i==0){
			dfs(ci,ck+1,cx,cv);
		}
		else{
			long long ncx = cx;
			ncx += (L / la[ck])*i;
			dfs(ci+i,ck+1,ncx,cv*t);
		}
		t *= i+1;
	}
}

int main(){
	
	cin>>K>>N;
	rep(i,K)En *= i+1;
	L = 1;
	rep(i,N)L = lcm(L,i+1);
	rep(i,N){
		if(i%2)la.push_back(i+1);
		else lb.push_back(i+1);
	}
	dfs(0,0,0,1);
	swap(la,lb);
	swap(mp,mp2);
	dfs(0,0,0,1);
	__int128 ans  =0;
	rep(i,K+1){
		for(auto x:mp[i]){
			__int128 A = x.first,B = x.second;
			__int128 X = B;
			__int128 Y = mp2[K-i][L-A];
			__int128 rem = 1LL;
			rep(j,K+1){
				if(j==0)continue;
				int cc = j;
				if(X%j==0){
					X /= j;
					continue;
				}
				if(Y%j==0){
					Y /= j;
					continue;
				}
				rem *= j;
			}
			X *= Y;
			X /= rem;
			ans += X;
		}
	}
	long long aa = ans;
	cout<<aa<<endl;
    return 0;
}
0