結果
問題 | No.152 貯金箱の消失 |
ユーザー |
👑 |
提出日時 | 2015-02-15 23:28:53 |
言語 | D (dmd 2.109.1) |
結果 |
AC
|
実行時間 | 51 ms / 5,000 ms |
コード長 | 2,101 bytes |
コンパイル時間 | 700 ms |
コンパイル使用メモリ | 124,868 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-12 02:19:30 |
合計ジャッジ時間 | 1,411 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 12 |
ソースコード
import core.thread;import std.conv, std.stdio, std.string;import std.algorithm, std.array, std.bigint, std.container, std.math, std.range, std.regex;// Inputclass EOFException : Throwable { this() { super("EOF"); } }string[] tokens;string readToken() { for (; tokens.empty; ) { tokens = readln.split; if (tokens.empty && stdin.eof) throw new EOFException; } auto token = tokens[0];tokens.popFront; return token; }int readInt() { return to!int(readToken); }long readLong() { return to!long(readToken); }real readReal() { return to!real(readToken); }// chmin/chmaxvoid chmin(T)(ref T t, in T f) { if (t > f) t = f; }void chmax(T)(ref T t, in T f) { if (t < f) t = f; }// Pairstruct Pair(S, T) {S x; T y;int opCmp( const Pair p) const { return (x < p.x) ? -1 : (x > p.x) ? +1 : (y < p.y) ? -1 : (y > p.y) ? +1 : 0; }int opCmp(ref const Pair p) const { return (x < p.x) ? -1 : (x > p.x) ? +1 : (y < p.y) ? -1 : (y > p.y) ? +1 : 0; }string toString() const { return "(" ~ to!string(x) ~ ", " ~ to!string(y) ~ ")"; }}auto pair(S, T)(inout(S) x, inout(T) y) { return Pair!(S, T)(x, y); }// Arrayint binarySearch(T)(in T[] as, in bool delegate(T) test) { int low = -1, upp = as.length; for (; low + 1 < upp; ) { int mid = (low + upp) >> 1; (test(as[mid]) ? low : upp) = mid; } return upp; }int lowerBound(T)(in T[] as, in T val) { return as.binarySearch((T a) { return (a < val); }); }int upperBound(T)(in T[] as, in T val) { return as.binarySearch((T a) { return (a <= val); }); }T[] unique(T)(in T[] as) { T[] bs; foreach (a; as) if (bs.empty || bs[$ - 1] != a) bs ~= a; return bs; }immutable MO = 1000003;int gcd(int a, int b) {return b ? gcd(b, a % b) : a;}int L;void main(string[] args) {try {for (; ; ) {L = readInt;int ans;for (int m = 1; 4 * m^^2 <= L; ++m) for (int n = 1; n < m; ++n) {if ((m - n) % 2 != 0 && gcd(m, n) == 1) {const a = m^^2 - n^^2;const b = 2 * m * n;const c = m^^2 + n^^2;if ((a + b + c) * 4 <= L) {++ans;}}}ans %= MO;writeln(ans);}} catch (EOFException) {}}