QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#127463#6634. Central SubsetoreoioiwyWA 10ms17664kbC++174.6kb2023-07-19 18:25:432023-07-19 18:25:44

Judging History

你现在查看的是最新测评结果

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-07-19 18:25:44]
  • 评测
  • 测评结果:WA
  • 用时:10ms
  • 内存:17664kb
  • [2023-07-19 18:25:43]
  • 提交

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)