結果

問題 No.428 小数から逃げる夢
ユーザー 37zigen37zigen
提出日時 2017-09-05 05:28:00
言語 Java21
(openjdk 21)
結果
AC  
実行時間 458 ms / 1,000 ms
コード長 9,464 bytes
コンパイル時間 3,316 ms
コンパイル使用メモリ 89,208 KB
実行使用メモリ 58,552 KB
最終ジャッジ日時 2024-11-06 21:34:17
合計ジャッジ時間 48,899 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 366 ms
54,220 KB
testcase_01 AC 363 ms
55,336 KB
testcase_02 AC 433 ms
55,496 KB
testcase_03 AC 368 ms
56,388 KB
testcase_04 AC 398 ms
56,424 KB
testcase_05 AC 362 ms
54,332 KB
testcase_06 AC 403 ms
56,528 KB
testcase_07 AC 367 ms
54,092 KB
testcase_08 AC 407 ms
58,012 KB
testcase_09 AC 458 ms
58,056 KB
testcase_10 AC 424 ms
58,104 KB
testcase_11 AC 400 ms
56,632 KB
testcase_12 AC 399 ms
56,832 KB
testcase_13 AC 398 ms
58,552 KB
testcase_14 AC 403 ms
56,556 KB
testcase_15 AC 416 ms
56,900 KB
testcase_16 AC 402 ms
56,756 KB
testcase_17 AC 396 ms
56,700 KB
testcase_18 AC 400 ms
57,004 KB
testcase_19 AC 423 ms
56,480 KB
testcase_20 AC 399 ms
56,756 KB
testcase_21 AC 418 ms
58,060 KB
testcase_22 AC 410 ms
58,128 KB
testcase_23 AC 402 ms
56,544 KB
testcase_24 AC 403 ms
56,824 KB
testcase_25 AC 436 ms
58,156 KB
testcase_26 AC 400 ms
56,636 KB
testcase_27 AC 404 ms
58,156 KB
testcase_28 AC 404 ms
56,892 KB
testcase_29 AC 424 ms
56,756 KB
testcase_30 AC 421 ms
56,684 KB
testcase_31 AC 410 ms
57,904 KB
testcase_32 AC 410 ms
57,120 KB
testcase_33 AC 400 ms
57,948 KB
testcase_34 AC 418 ms
56,792 KB
testcase_35 AC 438 ms
56,492 KB
testcase_36 AC 399 ms
57,364 KB
testcase_37 AC 403 ms
58,332 KB
testcase_38 AC 403 ms
56,504 KB
testcase_39 AC 409 ms
56,476 KB
testcase_40 AC 425 ms
56,516 KB
testcase_41 AC 418 ms
56,992 KB
testcase_42 AC 419 ms
56,664 KB
testcase_43 AC 422 ms
57,216 KB
testcase_44 AC 433 ms
56,792 KB
testcase_45 AC 413 ms
56,636 KB
testcase_46 AC 431 ms
56,364 KB
testcase_47 AC 415 ms
56,772 KB
testcase_48 AC 424 ms
57,920 KB
testcase_49 AC 406 ms
56,600 KB
testcase_50 AC 399 ms
57,876 KB
testcase_51 AC 420 ms
56,700 KB
testcase_52 AC 402 ms
56,636 KB
testcase_53 AC 399 ms
56,852 KB
testcase_54 AC 402 ms
56,632 KB
testcase_55 AC 402 ms
58,252 KB
testcase_56 AC 429 ms
57,164 KB
testcase_57 AC 413 ms
56,556 KB
testcase_58 AC 400 ms
57,080 KB
testcase_59 AC 410 ms
56,792 KB
testcase_60 AC 420 ms
57,388 KB
testcase_61 AC 431 ms
58,164 KB
testcase_62 AC 404 ms
58,272 KB
testcase_63 AC 427 ms
56,836 KB
testcase_64 AC 407 ms
56,912 KB
testcase_65 AC 404 ms
56,868 KB
testcase_66 AC 435 ms
56,748 KB
testcase_67 AC 407 ms
57,012 KB
testcase_68 AC 427 ms
57,028 KB
testcase_69 AC 410 ms
56,560 KB
testcase_70 AC 397 ms
58,012 KB
testcase_71 AC 453 ms
56,812 KB
testcase_72 AC 394 ms
56,592 KB
testcase_73 AC 398 ms
57,916 KB
testcase_74 AC 398 ms
57,948 KB
testcase_75 AC 399 ms
55,832 KB
testcase_76 AC 397 ms
58,000 KB
testcase_77 AC 423 ms
58,216 KB
testcase_78 AC 446 ms
58,168 KB
testcase_79 AC 433 ms
56,828 KB
testcase_80 AC 422 ms
56,684 KB
testcase_81 AC 391 ms
56,940 KB
testcase_82 AC 392 ms
56,752 KB
testcase_83 AC 398 ms
56,656 KB
testcase_84 AC 413 ms
56,620 KB
testcase_85 AC 410 ms
58,028 KB
testcase_86 AC 396 ms
55,768 KB
testcase_87 AC 404 ms
57,052 KB
testcase_88 AC 403 ms
57,024 KB
testcase_89 AC 412 ms
57,780 KB
testcase_90 AC 401 ms
56,652 KB
testcase_91 AC 397 ms
56,796 KB
testcase_92 AC 402 ms
58,192 KB
testcase_93 AC 406 ms
56,652 KB
testcase_94 AC 394 ms
56,592 KB
testcase_95 AC 411 ms
56,580 KB
testcase_96 AC 404 ms
56,788 KB
testcase_97 AC 397 ms
55,932 KB
testcase_98 AC 399 ms
55,356 KB
testcase_99 AC 406 ms
56,744 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import java.io.FileNotFoundException;
import java.util.ArrayList;
import java.util.Arrays;
import java.util.Scanner;

public class Main {
	class Int {
		long[] mag;
		int sig = 0;
		long long_mask = -1L;
		long int_mask = long_mask >>> 32;
		long radix = int_mask + 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);
		}

		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;
			mag = new long[] { int_mask & v, v >>> 32, };
		}

		// public Int(String str) {
		// if (str.equals("0")) {
		// mag = new long[] { 0 };
		// return;
		// }
		// if (str.charAt(str.length() - 1) == '-') {
		// sig = -1;
		// }

		// 1~999は2^10> より10ビット
		// 1000も10ビット
		// 2013まで10ビット
		// 2^10>10^3
		// 2^(10/3)>10
		// 1文字 2^(10/3)
		// 2文字 2^(10/3)*2
		// long numBits = (str.length() * 10 + 2) / 3;
		// mag = new long[(int) numBits];
		// for (int i = 0; i < str.length(); ++i) {
		//
		// add(new Int(str.charAt(i) - '0'));
		// }
		// }

		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;
				res = v >>> 32;
			}
			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;
				res = v >>> 32;
				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;
			long res = prd >>> 32;
			if (res > 0) {
				++intLen;
				mag = Arrays.copyOf(mag, 2);
				mag[1] = res;
			}
			sig *= a.sig;
		}

		// n^1.58
		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);
			// (a*r+b)(c*r+d)
			// (a+b)(c+d)r-(ac+bd)r+acr^2+db
			// XYr -(Z+W)r+Zr^2+W;
			// a=0かつb=1 より
			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;
		}

		void div(Int b) {
			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(tmp);
				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);
			mag = X.mag;
			sig = nextsig;
			intLen = X.intLen;
		}

		// num<<k;
		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;
			// tr("shift", mag, intLen);
		}

		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());
			// rs.get(i) = b^{2^i}
			// rs.back=b^{2^cnt} s.t b^{2^(cnt+1)}>=this
			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;
		}

		String toString(int radix, int cnt) {
			if (intLen <= 1) {
				if (radix == 10) {
					return String.valueOf(mag[0]);
				} else if (radix == 2) {
					return Long.toBinaryString(mag[0]);
				} else {
					throw new AssertionError();
				}
			}
			// radix==10として計算する
			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();
		}

		String toStringRadix2() {
			String ret = "";
			if (sig == 0)
				return "0";
			if (mag[intLen - 1] == 0)
				throw new AssertionError();
			for (int i = 0; i < intLen; ++i) {
				String s = Long.toBinaryString(mag[i]);
				while (s.length() < 32) {
					s = "0" + s;
				}
				ret = s + ret;
			}
			return ret;
		}
	}

	public void run() {
		Scanner sc = new Scanner(System.in);
		String str = "0.1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889909192939495969798991";
		// 1e190
		Int a = new Int(0);
		for (int i = 2; i < str.length(); ++i) {
			a.karatsuba(new Int(10));
			a.add(new Int(str.charAt(i) - '0'));
		}
		Int n = new Int(sc.nextInt());
		a.karatsuba(n);
		String ans = a.toString(10);
		if (ans.length() > 190) {
			System.out.println(
					ans.substring(0, ans.length() - 190) + "." + ans.substring(+ans.length() - 190, ans.length()));
		} else {
			System.out.println(
					0 + "." + ans.substring(+ans.length() - 190, ans.length()));
		}
	}

	public static void main(String[] args) throws FileNotFoundException {
		new Main().run();
	}

	static void tr(Object... objects) {
		System.out.println(Arrays.deepToString(objects));
	}
}
0