#include <bits/stdc++.h>
using namespace std;
class {
long long MultiplyMod(long long int a, long long int b, long long int mod) { //computes a * b % mod
long long int r = 0;
a %= mod, b %= mod;
while (b) {
if (b & 1)
r = (r + a) % mod;
b >>= 1, a = ((long long int) a << 1) % mod;
}
return r;
}
template<typename T>
T PowerMod(T a, T n, T mod) { // computes a^n % mod
T r = 1;
while (n) {
if (n & 1)
r = MultiplyMod(r, a, mod);
n >>= 1, a = MultiplyMod(a, a, mod);
}
return r;
}
template<typename T>
bool isprime(T n) { // Miller-Rabbin
// sprawdza czy n jest liczba pierwsza
const int pn = 9, p[] = {2, 3, 5, 7, 11, 13, 17, 19, 23};
for (int i = 0; i < pn; ++i)
if (n % p[i] == 0)
return n == p[i];
if (n < p[pn - 1])
return 0;
T s = 0, t = n - 1;
while (~t & 1)
t >>= 1, ++s;
for (int i = 0; i < pn; ++i) {
T pt = PowerMod<T>(p[i], t, n);
if (pt == 1)
continue;
bool ok = 0;
for (int j = 0; j < s && !ok; ++j) {
if (pt == n - 1)
ok = 1;
pt = MultiplyMod(pt, pt, n);
}
if (!ok)
return 0;
}
return 1;
}
public:
bool sprawdz() {
cin >> ws;
// bez wiodacych zer
if (cin.peek() == '0')
return false;
cin >> numer;
/////////////////////////////////////////////////////////////
long long int lewaczesc = 0, prawaczesc, lewospraw;
for (long long int pow10 = 1; pow10 < numer; pow10 *= 10) {
lewospraw = lewaczesc;
lewaczesc = numer / pow10;
prawaczesc = numer - (lewaczesc * pow10);
// wykryj i omin wiodace zera
if (lewospraw == (lewospraw / 10) * 10)
continue;
if (isprime<long long int>(prawaczesc) && isprime<long long int>(lewaczesc))
return true;
}
return false;
}
long long int numer;
} k;
int main()
{
cout << (k.sprawdz() ? "TAK" : "NIE");
return 0;
}
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 | #include <bits/stdc++.h> using namespace std; class { long long MultiplyMod(long long int a, long long int b, long long int mod) { //computes a * b % mod long long int r = 0; a %= mod, b %= mod; while (b) { if (b & 1) r = (r + a) % mod; b >>= 1, a = ((long long int) a << 1) % mod; } return r; } template<typename T> T PowerMod(T a, T n, T mod) { // computes a^n % mod T r = 1; while (n) { if (n & 1) r = MultiplyMod(r, a, mod); n >>= 1, a = MultiplyMod(a, a, mod); } return r; } template<typename T> bool isprime(T n) { // Miller-Rabbin // sprawdza czy n jest liczba pierwsza const int pn = 9, p[] = {2, 3, 5, 7, 11, 13, 17, 19, 23}; for (int i = 0; i < pn; ++i) if (n % p[i] == 0) return n == p[i]; if (n < p[pn - 1]) return 0; T s = 0, t = n - 1; while (~t & 1) t >>= 1, ++s; for (int i = 0; i < pn; ++i) { T pt = PowerMod<T>(p[i], t, n); if (pt == 1) continue; bool ok = 0; for (int j = 0; j < s && !ok; ++j) { if (pt == n - 1) ok = 1; pt = MultiplyMod(pt, pt, n); } if (!ok) return 0; } return 1; } public: bool sprawdz() { cin >> ws; // bez wiodacych zer if (cin.peek() == '0') return false; cin >> numer; ///////////////////////////////////////////////////////////// long long int lewaczesc = 0, prawaczesc, lewospraw; for (long long int pow10 = 1; pow10 < numer; pow10 *= 10) { lewospraw = lewaczesc; lewaczesc = numer / pow10; prawaczesc = numer - (lewaczesc * pow10); // wykryj i omin wiodace zera if (lewospraw == (lewospraw / 10) * 10) continue; if (isprime<long long int>(prawaczesc) && isprime<long long int>(lewaczesc)) return true; } return false; } long long int numer; } k; int main() { cout << (k.sprawdz() ? "TAK" : "NIE"); return 0; } |
English