結果

問題 No.800 四平方定理
ユーザー bal4ubal4u
提出日時 2019-04-13 07:22:45
言語 C
(gcc 12.3.0)
結果
AC  
実行時間 143 ms / 2,000 ms
コード長 1,481 bytes
コンパイル時間 274 ms
コンパイル使用メモリ 32,640 KB
実行使用メモリ 40,832 KB
最終ジャッジ日時 2024-09-15 06:49:34
合計ジャッジ時間 3,040 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 3 ms
5,376 KB
testcase_02 AC 3 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 4 ms
5,376 KB
testcase_08 AC 4 ms
5,376 KB
testcase_09 AC 3 ms
5,376 KB
testcase_10 AC 58 ms
20,864 KB
testcase_11 AC 63 ms
20,864 KB
testcase_12 AC 66 ms
23,168 KB
testcase_13 AC 51 ms
17,920 KB
testcase_14 AC 50 ms
19,072 KB
testcase_15 AC 66 ms
22,912 KB
testcase_16 AC 53 ms
19,456 KB
testcase_17 AC 51 ms
19,328 KB
testcase_18 AC 80 ms
27,136 KB
testcase_19 AC 77 ms
27,008 KB
testcase_20 AC 1 ms
5,376 KB
testcase_21 AC 1 ms
5,376 KB
testcase_22 AC 90 ms
27,136 KB
testcase_23 AC 141 ms
40,832 KB
testcase_24 AC 106 ms
33,152 KB
testcase_25 AC 143 ms
40,704 KB
testcase_26 AC 1 ms
5,376 KB
testcase_27 AC 1 ms
5,376 KB
testcase_28 AC 102 ms
32,768 KB
testcase_29 AC 121 ms
36,864 KB
testcase_30 AC 104 ms
32,768 KB
testcase_31 AC 95 ms
30,720 KB
testcase_32 AC 111 ms
33,920 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c: In function 'in':
main.c:10:14: warning: implicit declaration of function 'getchar_unlocked' [-Wimplicit-function-declaration]
   10 | #define gc() getchar_unlocked()
      |              ^~~~~~~~~~~~~~~~
main.c:18:24: note: in expansion of macro 'gc'
   18 |         int n = 0, c = gc();
      |                        ^~

ソースコード

diff # 

// yukicoder: No.800 四平方定理
// 2019.4.13 bal4u
// ハッシュテーブル x^2 + y^2 = w^2 + D - z^2

#include <stdio.h>
#include <math.h>

//// 高速入力
#if 1
#define gc() getchar_unlocked()
#define pc(c) putchar_unlocked(c)
#else
#define gc() getchar()
#define pc(c) putchar(c)
#endif
int in()    // 非負整数の入力
{
	int n = 0, c = gc();
	do n = 10 * n + (c & 0xf), c = gc(); while (c >= '0');
	return n;
}

//// ハッシュテーブル
#define HASHSIZ 19999999 //9999991
typedef struct { int s, cnt; } HASH;
HASH hash[HASHSIZ+5], *hashend = hash + HASHSIZ;

int lookup(int s)
{
	HASH *p = hash + s % HASHSIZ;
	while (p->s) {
		if (p->s == s) return p->cnt;
		if (++p == hashend) p = hash;
	}
	return 0;
}

void insert(int s)
{
	HASH *p = hash + s % HASHSIZ;
	while (p->s) {
		if (p->s == s) { p->cnt++; return; }
		if (++p == hashend) p = hash;
	}
	p->s = s, p->cnt = 1;
}

#define MAX 2000
int n2[MAX+5];

int main()
{
	int i, k, N, D, x, y, z, zmax, w, ans;

	for (i = 0; i < MAX; i++) {
		n2[i+1] = n2[i] + (i<<1) + 1;
	}
	
	N = in(), D = in();

	// w^2 + D - z^2
	for (w = 1; w <= N; w++) {
		zmax = (int)sqrt(n2[w] + D - 2);
		if (zmax > N) zmax = N;
		k = n2[w] + D;
		for (z = 1; z <= zmax; z++) if (k-n2[z] > 1) insert(k - n2[z]);
	}

	// x ^ 2 + y ^ 2
	ans = 0;
	for (x = 1; x <= N; x++) for (y = x + 1; y <= N; y++) {
		ans += lookup(n2[x] + n2[y]) << 1;
	}
	for (x = 1; x <= N; x++) ans += lookup(n2[x] << 1);

	printf("%d\n", ans);
	return 0;
}
0