結果
問題 | 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;}