結果
問題 | No.381 名声値を稼ごう Extra |
ユーザー | 37zigen |
提出日時 | 2017-09-05 18:58:46 |
言語 | Java21 (openjdk 21) |
結果 |
TLE
|
実行時間 | - |
コード長 | 10,907 bytes |
コンパイル時間 | 2,937 ms |
コンパイル使用メモリ | 90,576 KB |
実行使用メモリ | 85,624 KB |
最終ジャッジ日時 | 2024-11-06 22:28:23 |
合計ジャッジ時間 | 12,476 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 55 ms
37,184 KB |
testcase_01 | TLE | - |
ソースコード
import java.io.BufferedReader; import java.io.FileNotFoundException; import java.io.IOException; import java.io.InputStream; import java.io.InputStreamReader; import java.io.PrintStream; import java.io.PrintWriter; import java.util.ArrayList; import java.util.Arrays; import java.util.Iterator; public class Main { class Int { long[] mag; int sig = 0; // 2進数 // long long_mask = -1L; // long int_mask = long_mask >>> 32; // long radix = int_mask + 1; // 10進数 long radix = (long) 1e9; // long radix = -1; int intLen = 0; // public Int(int v) { // intLen = v == 0 ? 0 : 1; // mag = new long[] { (long) Math.abs(v) }; // sig = (int) Math.signum(v); // } // // radix == 10で取得を想定 public Int(String s) { if (s.length() == 1 && s.charAt(0) == 0) { sig = 0; intLen = 0; mag = new long[] { 0 }; return; } if (s.indexOf(0) == '-') sig = -1; else sig = 1; int pos = 0; mag = new long[s.length() / 9 + 1]; for (int i = s.length() - 1; i >= (s.indexOf(0) == '-' ? 1 : 0); i -= 9) { long v = 0; for (int j = Math.max(i - 8, 0); j <= i; ++j) { v *= 10; v += s.charAt(j) - '0'; } mag[pos] = v; ++pos; ++intLen; } } public Int(long v) { if (v == 0) { intLen = 0; sig = 0; mag = new long[] { 0 }; return; } sig = (int) Math.signum(v); v = (long) (Math.abs(v)); // intLen = (v >= (1L << 32)) ? 2 : 1; intLen = (v >= radix) ? 2 : 1; // mag = new long[] { int_mask & v, v >>> 32, }; mag = new long[] { v % radix, v / radix }; } void add(Int a) { if (a.sig == 0) { return; } if (sig * a.sig >= 0) { unsignedAdd(a); if (sig == 0) { sig = a.sig; } } else { Int u = a.copy(); u.sig *= -1; sub(u); } } void sub(Int a) { if (a.sig == 0) return; if (sig != a.sig) { unsignedAdd(a); if (sig == 0) { sig = (-1) * a.sig; } } else { if (sig * compare(a) > 0) { unsignedSub(a); } else { Int t = a.copy(); t.unsignedSub(this); sig *= -t.sig; mag = t.mag; intLen = t.intLen; } } } void unsignedAdd(Int a) { long res = 0; if (mag.length < a.mag.length) { mag = Arrays.copyOf(mag, a.mag.length); } intLen = Math.max(intLen, a.intLen); for (int i = 0; i < a.intLen; ++i) { long v = mag[i] + a.mag[i] + res; // mag[i] = int_mask & v; mag[i] = v % radix; // res = v >>> 32; res = v / radix; } for (int i = a.intLen; res != 0; ++i) { if (i >= mag.length) { mag = Arrays.copyOf(mag, 1 + mag.length); } long v = mag[i] + res; // mag[i] = int_mask & v; mag[i] = v % radix; // res = v >>> 32; res = v / radix; intLen = Math.max(intLen, i + 1); } } void unsignedSub(Int a) { if (sig == 0) { mag = Arrays.copyOf(a.mag, a.mag.length); sig = -1 * a.sig; } long res = 0; for (int i = 0; i < a.intLen; ++i) { long t = mag[i] - a.mag[i] - res; if (t < 0) { res = 1; t += radix; } else { res = 0; } mag[i] = t; } for (int i = a.intLen; res != 0; ++i) { long t = mag[i] - res; if (t < 0) { res = 1; t += radix; } else { res = 0; } mag[i] = t; } boolean f = false; for (int i = intLen - 1; i >= 0; --i) { if (mag[i] != 0) { intLen = i + 1; f = true; break; } } if (!f) { sig = 0; intLen = 0; } } int compare(Int a) { if (sig != a.sig) { return Integer.compare(sig, a.sig); } if (intLen != a.intLen) { return sig * Integer.compare(intLen, a.intLen); } for (int i = intLen - 1; i >= 0; --i) { if (mag[i] == a.mag[i]) continue; return sig * Long.compare(mag[i], a.mag[i]); } return 0; } void singleMul(Int a) { if (intLen > 1 || a.intLen > 1) { throw new AssertionError(); } if (a.sig == 0) { intLen = 0; sig = 0; mag = new long[] { 0 }; return; } long prd = mag[0] * a.mag[0]; // mag[0] = int_mask & prd; mag[0] = prd % radix; // long res = prd >>> 32; long res = prd / radix; if (res > 0) { ++intLen; mag = Arrays.copyOf(mag, 2); mag[1] = res; } sig *= a.sig; } // n^1.58 // ( a* r+ b )( c*r + d ) // = ( a + b )( c + d )r-( ac + bd )r + acr^2 + db // = XYr - ( Z + W )r + Zr^2 + W; void karatsuba(Int o) { int len_ = Math.max(intLen, o.intLen); if (len_ <= 1) { singleMul(o); return; } int Len = 1; while (Len < len_) { Len *= 2; } Int ret = new Int(0); Int a = copy(Len / 2, Len); Int b = copy(0, Len / 2); Int c = o.copy(Len / 2, Len); Int d = o.copy(0, Len / 2); Int X = a.copy(); X.add(b); Int Y = c.copy(); Y.add(d); X.karatsuba(Y); X.shift(Len / 2); ret.add(X); Int Z = a.copy(); Z.karatsuba(c); Int W = b.copy(); W.karatsuba(d); ret.add(W); Int tmp = Z.copy(); tmp.shift(Len); ret.add(tmp); Z.add(W); Z.shift(Len / 2); ret.sub(Z); this.sig = ret.sig; this.mag = Arrays.copyOf(ret.mag, ret.mag.length); this.intLen = ret.intLen; } // Newton method // n^1.58 void div(Int b) { Int a = copy(); if (b.sig == 0) throw new ArithmeticException(); int nextsig = sig * b.sig; sig = (int) (Math.abs(sig)); b.sig = (int) (Math.abs(b.sig)); if (compare(b) < 0) { sig = 0; mag = new long[] { 0 }; intLen = 0; return; } Int M = new Int(1); Int X = new Int(1); M.shift(intLen + 1); Int pre = X.copy(); while (true) { Int tmp = X.copy(); tmp.karatsuba(X); tmp.karatsuba(b); tmp.shift(-M.intLen); X.add(X); X.sub(tmp); if (X.compare(pre) == 0) { break; } pre = X.copy(); } X.karatsuba(this); X.shift(-M.intLen); { Int d = X.copy(); d.karatsuba(b); if (a.compare(d) < 0) { X.sub(new Int(1)); } } mag = X.mag; sig = nextsig; intLen = X.intLen; } void shift(int k) { if (sig == 0) return; if (intLen + k <= 0) { mag = new long[] { 0 }; sig = 0; intLen = 0; return; } Int u = copy(); mag = new long[intLen + k]; if (k >= 0) System.arraycopy(u.mag, 0, mag, k, intLen); else System.arraycopy(u.mag, -k, mag, 0, intLen + k); intLen += k; } Int copy(int l, int r) { Int a = copy(); a.mag = new long[r - l]; r = Math.min(r, a.intLen); if (l >= r) { return new Int(0); } System.arraycopy(mag, l, a.mag, 0, r - l); a.intLen = r - l; return a; } Int copy() { Int ret = new Int(0); ret.sig = sig; ret.mag = Arrays.copyOf(mag, mag.length); ret.intLen = intLen; return ret; } void check() { // v=0->sig==0を満たすか? boolean f = true; for (int i = 0; i < mag.length; ++i) { f &= mag[i] == 0; } if (f && sig != 0) { throw new AssertionError(); } if (sig == 0) return; // intLenは正しいか? f = mag[intLen - 1] > 0; for (int i = mag.length - 1; i > intLen - 1; --i) { f &= mag[i] == 0; } if (!f) { throw new AssertionError(); } } void mod(Int mod) { Int tmp = copy(); tmp.div(mod); tmp.karatsuba(mod); sub(mod); } ArrayList<Int> rs = new ArrayList<>(); // rsサイズの変更には対応していない。 String toString(int radix) { if (sig == 0) return "0"; rs.clear(); Int r = new Int(radix); Int b = r.copy(); int cnt = 0; rs.add(b.copy()); while (true) { ++cnt; b.karatsuba(b); rs.add(b.copy()); Int tmp = b.copy(); tmp.karatsuba(b); if (tmp.compare(this) >= 0) { break; } } String str = toString(radix, cnt); while (str.charAt(0) == '0') { str = str.substring(1, str.length()); } if (sig == -1) str = "-" + str; return str; } long popCount(long MODULO) { if (sig == 0) return 0; rs.clear(); Int r = new Int(2); Int b = r.copy(); int cnt = 0; rs.add(b.copy()); while (true) { ++cnt; b.karatsuba(b); rs.add(b.copy()); Int tmp = b.copy(); tmp.karatsuba(b); if (tmp.compare(this) >= 0) { break; } } return popCount(MODULO, cnt); } // radix==10を想定 long popCount(long MODULO, int cnt) { if (intLen <= 2) { return Long.bitCount(mag[0] + mag[1] * radix); } Int b = rs.get(cnt); Int u = copy(); u.div(b); Int l = copy(); u.rs = rs; long us = u.popCount(MODULO, cnt - 1); u.karatsuba(b); l.sub(u); l.rs = rs; long ls = l.popCount(MODULO, cnt - 1); long ret = us + ls; if (ret >= MODULO) ret -= MODULO; return ret; } String toString(int radix, int cnt) { if (intLen <= 1) { if (radix == 10) { String ret = String.valueOf(mag[0]); while (ret.length() < 1 << (cnt + 1)) { ret = "0" + ret; } return ret; } else if (radix == 2) { return Long.toBinaryString(mag[0]); } else { throw new AssertionError(); } } Int b = rs.get(cnt); Int u = copy(); u.div(b); Int l = copy(); u.rs = rs; String us = u.toString(radix, cnt - 1); while (us.length() < 1 << cnt) { us = "0" + us; } u.karatsuba(b); l.sub(u); l.rs = rs; String ls = l.toString(radix, cnt - 1); while (ls.length() < 1 << cnt) { ls = "0" + ls; } return us + ls; } long toLong() { if (sig == 0) { return 0; } if (intLen == 1) { return sig * mag[0]; } if (intLen == 2) { return sig * (mag[0] + (mag[1] << 32)); } throw new AssertionError(); } } public void run() { Scanner sc = new Scanner(System.in); Int v = new Int(sc.next()); long MODULO = 1004535809L; System.out.println(v.popCount(MODULO)); } public static void main(String[] args) throws FileNotFoundException { new Main().run(); } static void tr(Object... objects) { System.out.println(Arrays.deepToString(objects)); } private static class Scanner { BufferedReader br; Iterator<String> it; Scanner(InputStream in) { br = new BufferedReader(new InputStreamReader(in)); } String next() throws RuntimeException { try { if (it == null || !it.hasNext()) it = Arrays.asList(br.readLine().split(" ")).iterator(); return it.next(); } catch (IOException e) { throw new IllegalStateException(); } } int nextInt() throws RuntimeException { return Integer.parseInt(next()); } long nextLong() throws RuntimeException { return Long.parseLong(next()); } double nextDouble() throws RuntimeException { return Double.parseDouble(next()); } void close() { try { br.close(); } catch (IOException e) { throw new IllegalStateException(); } } } private static class Printer extends PrintWriter { Printer(PrintStream out) { super(out); } } }