結果
問題 | No.800 四平方定理 |
ユーザー |
|
提出日時 | 2019-03-17 21:27:21 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 73 ms / 2,000 ms |
コード長 | 1,085 bytes |
コンパイル時間 | 534 ms |
コンパイル使用メモリ | 74,740 KB |
実行使用メモリ | 57,700 KB |
最終ジャッジ日時 | 2024-07-07 20:25:25 |
合計ジャッジ時間 | 2,439 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#include <iostream> #include <algorithm> #include <vector> #include <map> using namespace std; #define REP(i, n) for (int i = 0; i < (n); i++) const int MOD = 1e9 + 7; struct mint { int n; mint(int n_ = 0) : n(n_) {} }; mint operator+(mint a, mint b) { a.n += b.n; if (a.n >= MOD) a.n -= MOD; return a; } mint operator-(mint a, mint b) { a.n -= b.n; if (a.n < 0) a.n += MOD; return a; } mint operator*(mint a, mint b) { return (long long)a.n * b.n % MOD; } mint &operator+=(mint &a, mint b) { return a = a + b; } mint &operator-=(mint &a, mint b) { return a = a - b; } mint &operator*=(mint &a, mint b) { return a = a * b; } int C1[2000*2000 + 2000*2000 + 1]; int C2[2000*2000 + 2000*2000 + 1]; int main() { int N, D; cin >> N >> D; for (int i = 1; i <= N; i++) { for (int j = 1; j <= N; j++) { C1[i*i + j*j]++; int k = i*i - j*j + D; if (0 <= k && k <= 2000*2000 + 2000*2000) { C2[k]++; } } } long long ans = 0; for (int i = 0; i <= 2000*2000 + 2000*2000; i++) { ans += C1[i] * C2[i]; } cout << ans << endl; }