結果

問題 No.800 四平方定理
ユーザー iicafiaxusiicafiaxus
提出日時 2019-03-17 21:44:10
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 1,072 ms / 2,000 ms
コード長 1,169 bytes
コンパイル時間 995 ms
コンパイル使用メモリ 131,212 KB
実行使用メモリ 66,300 KB
最終ジャッジ日時 2024-06-13 04:19:43
合計ジャッジ時間 15,787 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.stdio, std.conv, std.string;
import std.array, std.range, std.algorithm, std.container;
import std.math, std.random, std.bigint, std.datetime, std.format;
string read(){ static string[] ss; while(!ss.length) ss = readln.chomp.split; string res = ss[0]; ss.popFront; return res; }

int DEBUG_LEVEL = 0;
void print()(){ writeln(""); }
void print(T, A ...)(T t, lazy A a){ write(t), print(a); }
void print(int level, T, A ...)(T t, lazy A a){ if(level <= DEBUG_LEVEL) print(t, a); }

void main(string[] args){
	if(args.length > 1 && args[1] == "-debug"){
		if(args.length > 2 && args[2].isNumeric) DEBUG_LEVEL = args[2].to!int;
		else DEBUG_LEVEL = 1;
	}
	
	/*
	
	xx + yy をカウント
	ww - zz + D をカウント ＊0以上のもののみ
	
	*/
	
	
	long n = read.to!long;
	long d = read.to!long;
	
	long[long] acnt, bcnt;
	
	long xx, yy, k;
	foreach(x; 1 .. n + 1){
		xx = x * x;
		foreach(y; 1 .. n + 1){
			k = xx + y * y;
			if(k !in acnt) acnt[k] = 1;
			else acnt[k] += 1;
		}
	}
	
	long ans;
	long zz;
	foreach(z; 1 .. n + 1){
		zz = z * z;
		foreach(w; 1 .. n + 1){
			k = zz - w * w + d;
			if(k in acnt) ans += acnt[k];
		}
	}
	
	ans.writeln;
	
}
0