#921593 (C++17) No.2510 Six Cube Sum Counting

• 問題: No.2510 Six Cube Sum Counting
• ユーザー: 👑 binap
• 提出日時: 2023-10-21 00:51:55
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: TLE
• 実行時間: -
• コード長: 1,811 bytes
• コンパイル時間: 10,638 ms
• コンパイル使用メモリ: 442,296 KB
• 実行使用メモリ: 362,356 KB
• 最終ジャッジ日時: 2023-10-21 00:53:35
• 合計ジャッジ時間: 92,048 ms
• ジャッジサーバーID （参考情報）: judge14 / judge12

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	3 ms; 4,328 KB
testcase_01	AC	2 ms; 4,332 KB
testcase_02	TLE	-
testcase_03	AC	2 ms; 4,332 KB
testcase_04	AC	2 ms; 4,332 KB
testcase_05	AC	3 ms; 4,332 KB
testcase_06	TLE	-
testcase_07	TLE	-
testcase_08	TLE	-
testcase_09	TLE	-
testcase_10	TLE	-
testcase_11	TLE	-
testcase_12	AC	3,953 ms; 362,356 KB
testcase_13	TLE	-
testcase_14	AC	3,930 ms; 362,356 KB
testcase_15	AC	2 ms; 4,332 KB
testcase_16	AC	2 ms; 4,332 KB
testcase_17	AC	3 ms; 4,332 KB
testcase_18	AC	3,896 ms; 362,356 KB
testcase_19	AC	3,917 ms; 362,356 KB
testcase_20	AC	3,011 ms; 221,360 KB
testcase_21	TLE	-
testcase_22	AC	3,755 ms; 321,436 KB
testcase_23	AC	3,272 ms; 252,776 KB
testcase_24	TLE	-
testcase_25	AC	3,947 ms; 362,356 KB
testcase_26	AC	3,963 ms; 362,356 KB
testcase_27	AC	3,982 ms; 359,188 KB
testcase_28	AC	3,937 ms; 362,356 KB
testcase_29	AC	3,844 ms; 339,916 KB

ソースコード

#include<bits/stdc++.h>
#include<atcoder/all>
#include <boost/multiprecision/cpp_int.hpp>
#define rep(i,n) for(int i=0;i<n;i++)
using namespace std;
using namespace atcoder;
typedef long long ll;
typedef vector<int> vi;
typedef vector<long long> vl;
typedef vector<vector<int>> vvi;
typedef vector<vector<long long>> vvl;
typedef long double ld;

template <int m>
std::ostream& operator<<(std::ostream& os, const atcoder::static_modint<m>& a) {os << a.val(); return os;}
template<typename T>
istream& operator>>(istream& is, vector<T>& v){int n = v.size(); assert(n > 0); rep(i, n) is >> v[i]; return is;}
template<typename T>
ostream& operator<<(ostream& os, const vector<T>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : " "); return os;}
template <typename T>
ostream& operator<<(ostream& os, const vector<vector<T>>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : ""); return os;}

namespace mp = boost::multiprecision;

int main(){
	int n;
	cin >> n;
	auto cu = [&](int i){
		return i * i * i;
	};
	unordered_map<int, mp::cpp_int> ma;
	for(int a = 0; a <= 300; a++){
		for(int b = a; b <= 300; b++){
			for(int c = b; c <= 300; c++){
				int x = cu(a) + cu(b) + cu(c);
				if(x > n) break;
				ma[x] *= 302;
				ma[x] += (c + 1);
			}
		}
	}
	ll ans = 0;
	for(int d = 0; d <= 300; d++){
		if(3 * cu(d) > n) break;
		if(4 * cu(d) + 54000000 < n) continue;
		for(int e = d; e <= 300; e++){
			if(cu(d) + 2 * cu(e) > n) break;
			if(4 * cu(d) + cu(e) + 27000000 < n) continue;
			for(int f = e; f <= 300; f++){
				int x = cu(d) + cu(e) + cu(f);
				if(x > n) break;
				int y = n - x;
				if(ma.find(y) != ma.end()){
					mp::cpp_int z = ma[y];
					while(z > 0){
						if(((z % 302) - 1) <= d) ans++;
						z /= 302;
					}
				}
			}
		}
	}
	cout << ans;
	return 0;
}