#include <algorithm>
#include <cstdio>
using namespace std;
using ull = unsigned long long;
const int MAX = 1000000;
int cnt;
unsigned int P[25];
int n;
struct uint128 {
ull b;
uint128(int x) {
b = x;
}
bool operator<=(ull x) const {
return b <= x;
}
bool operator>(ull x) const {
return b > x;
}
uint128 operator*(unsigned long long x) const {
uint128 res = 0;
ull u = b >> 32;
ull l = b & 0xffffffff;
ull r = u * x + ((l * x) >> 32);
if (r >> 32) {
res.b = 1ll<<63;
} else {
res.b = b * x;
}
return res;
}
operator ull() const {
return b;
}
};
unsigned long long H[MAX];
int hsize;
void go() {
H[hsize++] = -1;
while (hsize) {
auto t = -H[0];
H[MAX-1-cnt++] = t;
pop_heap(H, H + hsize);
hsize--;
for (int i = 1; i < n; i++){
int r = P[i];
auto y = (ull)t * r;
if (y >= 1000000000ll) {
break;
}
bool vis = false;
for (int j = i + 1; j < n; j++) {
if (y % P[j] == 0) {
vis = true;
break;
}
}
if (!vis) {
H[hsize++] = -y;
push_heap(H, H + hsize);
}
}
}
copy_n(H+MAX-cnt, cnt, H);
reverse(H, H+cnt);
}
int main() {
ull k;
scanf("%d %llu\n", &n, &k);
for (int i=0; i<n; i++) {
scanf("%u", P+i);
}
sort(P, P+n);
go();
ull res = 1;
for (uint128 p1 = 1; p1 <= k; p1 = p1 * (ull)P[0]) {
for (int p = 0; p < n; p++) {
auto a = p1 * (ull)(p ? P[p] : 1);
int j = cnt - 1;
while (a * H[j/2] > k && j > 0) j /= 2;
for (int i = 0; i <= j; i++) {
while (i <= j && a * H[i] * H[j] > k) {
j--;
}
if (i <= j) {
res = max(res, ull(a) * H[i] * H[j]);
}
}
}
}
printf("%lld\n", res);
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 98 99 100 101 102 103 104 105 106 | #include <algorithm> #include <cstdio> using namespace std; using ull = unsigned long long; const int MAX = 1000000; int cnt; unsigned int P[25]; int n; struct uint128 { ull b; uint128(int x) { b = x; } bool operator<=(ull x) const { return b <= x; } bool operator>(ull x) const { return b > x; } uint128 operator*(unsigned long long x) const { uint128 res = 0; ull u = b >> 32; ull l = b & 0xffffffff; ull r = u * x + ((l * x) >> 32); if (r >> 32) { res.b = 1ll<<63; } else { res.b = b * x; } return res; } operator ull() const { return b; } }; unsigned long long H[MAX]; int hsize; void go() { H[hsize++] = -1; while (hsize) { auto t = -H[0]; H[MAX-1-cnt++] = t; pop_heap(H, H + hsize); hsize--; for (int i = 1; i < n; i++){ int r = P[i]; auto y = (ull)t * r; if (y >= 1000000000ll) { break; } bool vis = false; for (int j = i + 1; j < n; j++) { if (y % P[j] == 0) { vis = true; break; } } if (!vis) { H[hsize++] = -y; push_heap(H, H + hsize); } } } copy_n(H+MAX-cnt, cnt, H); reverse(H, H+cnt); } int main() { ull k; scanf("%d %llu\n", &n, &k); for (int i=0; i<n; i++) { scanf("%u", P+i); } sort(P, P+n); go(); ull res = 1; for (uint128 p1 = 1; p1 <= k; p1 = p1 * (ull)P[0]) { for (int p = 0; p < n; p++) { auto a = p1 * (ull)(p ? P[p] : 1); int j = cnt - 1; while (a * H[j/2] > k && j > 0) j /= 2; for (int i = 0; i <= j; i++) { while (i <= j && a * H[i] * H[j] > k) { j--; } if (i <= j) { res = max(res, ull(a) * H[i] * H[j]); } } } } printf("%lld\n", res); return 0; } |
English