QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#127463 | #6634. Central Subset | oreoioiwy | WA | 10ms | 17664kb | C++17 | 4.6kb | 2023-07-19 18:25:43 | 2023-07-19 18:25:44 |
Judging History
answer
#include <bits/stdc++.h>
#define int long long
using namespace std;
typedef long long ll;
// #define tpyeinput int
// inline char nc() {static char buf[1000000],*p1=buf,*p2=buf;return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;}
// inline void read(tpyeinput &sum) {char ch=nc();sum=0;while(!(ch>='0'&&ch<='9')) ch=nc();while(ch>='0'&&ch<='9') sum=(sum<<3)+(sum<<1)+(ch-48),ch=nc();}
template<typename T>
void read(T &x) {
int f = 1;
x = 0;
char ch = getchar();
while (ch < '0' || ch > '9') {
if (ch == '-')f = -1;
ch = getchar();
}
while (ch >= '0' && ch <= '9') {
x = x * 10 + (ch ^ 48);
ch = getchar();
}
x *= f;
}
const int N = 1e6+10;
int n, m, K;
int a[N], f[N][30];
int M[N],MM;
bool prime[N];
vector<int> primes, fac;
void ST_pre() {
for(int i=1; i<=m; i++) f[i][0] = M[i];
int t = log(m)/ log(2) + 1;
for(int j=1; j<t; j++) {
for(int i=1; i<=m-(1<<j)+1; i++) {
f[i][j] = min(f[i][j-1], f[i+(1<<(j-1))][j-1]);
}
}
}
int query(int l, int r) {
int k = log(r-l+1)/log(2);
return min(f[l][k], f[r-(1<<k)+1][k]);
}
int ksm(int x,int y)
{
int res = 1;
while(y)
{
if(y & 1)
{
res *= x;
}
x *= x;
y >>= 1;
}
return res;
}
void preprim()
{
for(int i = 2; i <= N; i++)
{
if(!prime[i])
{
primes.push_back(i);
for(int j = i + i; j <= N; j += i)
{
prime[j] = true;
}
}
}
for(int i = 2; i <= N; i++)
{
if(!prime[i])
{
fac.push_back(i);
}else
{
primes.push_back(i);
}
}
}
void precost()
{
for(int i = 0; fac[i] <= m; i++)
{
int mn = M[i];
vector<pair<int,int>> prims;
int ii = i;
for(int j = 0; primes[j] <= sqrt(i); j++)
{
if(ii % j == 0)
{
prims.push_back({j,0});
while(ii % j == 0)
{
prims.back().second++;
ii /= j;
}
}
}
if(ii > 1)
{
prims.push_back({ii,1});
}
if(prims.size() == 1)
{
auto [x,y] = prims[0];
for(int j = 1; j <= (y >> 1); j++)
{
mn = min(mn, M[ksm(x,j)] + M[ksm(x,y-j)]);
}
}else
{
int sum = 0;
for(auto [x,y] : prims)
{
sum += M[ksm(x,y)];
}
mn = min(mn, sum);
}
M[i] = mn;
}
for(int i=2; i<=m;i++)
{
for(int j=(m+1)/i - 1; j <= m; j++)
{
if(i * j <= m) continue;
MM = min(MM, M[i] + M[j]);
}
}
}
void pre_M() {
for(int i=1; i<=m; i++) {
for(int j=1; j<i; j++) {
if(i % j == 0) {
M[i] = min(M[i], M[j] + M[i/j]);
}
}
}
for(int i=2; i<=m;i++)
{
for(int j=(m+1)/i - 1; j <= m; j++)
{
if(i * j <= m) continue;
MM = min(MM, M[i] + M[j]);
}
}
}
ll check(int mid) {
if(mid == a[n/2+1]) return 0;
if(mid > a[n/2+1]) return -1;
ll cost = 0;
int l = lower_bound(a+1, a+1+n, mid+1) - a;
// int ll = lower_bound(a+1, a+1+n, mid) - a;
int cnt = l-1;
// if(cntt > n/2) return -1;
while(cnt < n/2+1) {
int tmp = mid ? (a[l]-1)/mid+1 : a[l]+1;
if(tmp > m)
{
cost += MM;
cnt++; l++;
continue;
}
while(a[l]/tmp > mid) tmp++;
while(a[l]/(tmp-1) <= mid) tmp--;
cost += min(MM,query(tmp, m));
cnt++; l++;
}
return cost;
}
void solve() {
read(n); read(m); read(K);
MM = 1e18;
for(int i=1; i<=n; i++) read(a[i]);
for(int i=1; i<=m; i++) read(M[i]);
// M_pre();
precost();
// pre_M();
ST_pre();
sort(a+1, a+1+n);
int l = 0, r = a[n/2+1];
while(l < r) {
int mid = l + r >> 1;
ll cost = check(mid);
if(cost >= 0 && cost <= K) r = mid;
else l = mid + 1;
}
printf("%lld\n", l);
}
signed main() {
#ifdef LOCAL_TEST
freopen("test.in", "r", stdin);
#endif
int times = 1;
read(times);
preprim();
while(times--) {
solve();
}
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 0
Wrong Answer
time: 10ms
memory: 17664kb
input:
2 4 3 1 2 2 3 3 4 6 7 1 2 2 3 3 1 1 4 4 5 5 6 6 4
output:
3 0
result:
wrong answer Integer 3 violates the range [1, 2] (test case 1)