結果

問題 No.2510 Six Cube Sum Counting
ユーザー 👑 binap
提出日時 2023-10-21 00:51:55
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3,856 ms / 4,000 ms
コード長 1,811 bytes
コンパイル時間 9,402 ms
コンパイル使用メモリ 426,688 KB
最終ジャッジ日時 2025-02-17 11:55:35
ジャッジサーバーID
（参考情報） 		judge3 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0