Mon, 04 Nov 2019 17:21:22 GMT

tmp/s01.cpp

@@ -0,0 +1,61 @@
+#define _CRT_SECURE_NO_WARNINGS
+#define _SILENCE_CXX17_C_HEADER_DEPRECATION_WARNING
+#include <bits/stdc++.h>
+using namespace std;
+#define CRP(t, x) const t &x
+#define OPX(op, t, x) operator op(CRP(t, x))
+#define OPL(t, x) bool OPX(<, t, x) const
+#define FIL(x, v) memset(x, v, sizeof(x))
+#define CLR(x) FIL(x, 0)
+#define NE1(x) FIL(x, -1)
+#define INF(x) FIL(x, 0x3f)
+#ifndef _DEBUG
+#define _DEBUG 0
+#endif // !_DEBUG
+#define IFD if (_DEBUG)
+typedef int64_t ll, i64;
+typedef uint64_t ull, u64;
+template <typename T>
+using enable_if_arithmetic = typename enable_if<is_arithmetic<T>::value>::type;
+template <typename T>
+using enable_if_integral = typename enable_if<is_integral<T>::value>::type;
+inline char getchar(int)
+{
+    static char buf[64 << 20], *S = buf, *T = buf;
+    if (S == T) T = fread(S = buf, 1, 64 << 20, stdin) + S;
+    return S == T ? EOF : *S++;
+}
+template <typename T, typename = enable_if_integral<T>>
+inline void read(T &x)
+{
+    int ch = x = 0, f = 1;
+    while (!isdigit(ch = getchar()))
+        if (ch == '-') f = -1;
+    for (; isdigit(ch); ch = getchar()) x = x * 10 + ch - '0';
+    x *= f;
+}
+template <typename T, typename... Args, typename = enable_if_integral<T>>
+inline void read(T &x, Args &... args) { read(x), read(args...); }
+ll gcd(ll a, ll b) { return b == 0 ? a : gcd(b, a % b); }
+ll fpow(ll a, ll b, ll m)
+{
+    ll r = 1;
+    for (; b; b >>= 1, a = a * a % m)
+        if (b & 1) r = r * a % m;
+    return r;
+}
+constexpr double eps = 1e-8;
+inline int sgn(double x) { return x > eps ? 1 : x < -eps ? -1 : 0; }
+const int N = 1e4 + 50;
+int minfac[N], primes[N], pcnt;
+int main()
+{
+    for (ll i = 2; i < N; i++)
+    {
+        if (minfac[i] == 0)
+            minfac[i] = i, primes[pcnt++] = i;
+        for (int j = 0; j < pcnt && primes[j] <= minfac[i] && i * primes[j] < N; j++)
+            minfac[i * primes[j]] = primes[j];
+    }
+    return 0;
+}