結果

問題 No.2510 Six Cube Sum Counting
ユーザー 沙耶花沙耶花
提出日時 2023-10-20 21:50:09
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 1,542 bytes
コンパイル時間 4,558 ms
コンパイル使用メモリ 261,284 KB
実行使用メモリ 647,168 KB
最終ジャッジ日時 2023-10-20 21:51:28
合計ジャッジ時間 20,619 ms
ジャッジサーバーID
（参考情報） 		judge10 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 MLE -
testcase_07 MLE -
testcase_08 MLE -
testcase_09 MLE -
testcase_10 MLE -
testcase_11 MLE -
testcase_12 MLE -
testcase_13 MLE -
testcase_14 MLE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 MLE -
testcase_19 MLE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 MLE -
testcase_26 MLE -
testcase_27 RE -
testcase_28 MLE -
testcase_29 MLE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

int c0[82000000],c1[82000000];

int main(){
	//?
	long long X;
	cin>>X;
	/*
	for(int i=0;i<10;i++){
		for(int j=i;j<10;j++){
			for(int k=j;k<10;k++){
				for(int l=k;l<10;l++){
					for(int m=l;m<10;m++){
						for(int n=m;n<10;n++){
							if(i*i*i+j*j*j+k*k*k+l*l*l+m*m*m+n*n*n==444){
								cout<<i<<j<<k<<l<<m<<n<<endl;
							}
						}
					}
				}
			}
		}
	}
	*/
	
	
	long long ans = 0;
	c0[0]++;
	rep(i,301){
		for(int j=i;j<=300;j++){
			for(int k=j;k<=300;k++){
				c1[i*i*i+j*j*j+k*k*k]++;
			}
		}
	}
	
	long long cur = 0;
	rep(i,82000000){
		if(X-i>=82000000)continue;
		cur += (long long)c0[i] * c1[X-i];
	}
	
	ans += cur;
	for(int i=1;i<=300;i++){
		for(int j=i;j>=i;j--){
			for(int k=j;k>=0;k--){
				for(int l=k;l>=0;l--){
					int nn = j*j*j+k*k*k+l*l*l;
					if(X-nn<82000000&&X-nn>=0)cur -= (long long)c0[nn] * c1[X-nn];
					c0[nn]++;
					if(X-nn<82000000&&X-nn>=0)cur += (long long)c0[nn] * c1[X-nn];
				}
			}
		}
		
		for(int j=i-1;j<=i-1;j++){
			for(int k=j;k<=300;k++){
				for(int l=k;l<=300;l++){
					int nn = j*j*j+k*k*k+l*l*l;
					if(X-nn<82000000&&X-nn>=0)cur -= (long long)c1[nn] * c0[X-nn];
					c1[nn]--;
					if(X-nn<82000000&&X-nn>=0)cur += (long long)c1[nn] * c0[X-nn];
					
				}
			}
		}
		
		ans += cur;
		
	}
	cout<<ans<<endl;
	
	return 0;
}
0