牛客网练习赛15 C、吉姆的奇思妙想（数学单调性 + 二分法）

最新推荐文章于 2020-07-23 11:20:53 发布

隆1

最新推荐文章于 2020-07-23 11:20:53 发布

阅读量245

点赞数

分类专栏：数学文章标签：牛客练习赛15 数学公式题

本文链接：https://blog.csdn.net/qq_38786088/article/details/79938883

版权

数学专栏收录该内容

13 篇文章 0 订阅

订阅专栏

题目链接

这里写图片描述

题意： 给你正整数a、b，求出s取某个整数时，Ei的最小值。（各参数的范围相当复杂）

E i = a i \times \sum 1 \leq j \leq L, d e g j \leq s (d e g 2 j \times f r e q j) + b i \times \sum 1 \leq j \leq L, d e g j > s (M \times f r e q j) = \sum 1 \leq j \leq L, d e g j \leq s (a i \times d e g 2 j \times f r e q j - b i \times M \times f r e q j) + b i \times M \times \sum 1 \leq j \leq L f r e q j = \sum 1 \leq j \leq L, d e g j \leq s (a i \times d e g 2 j - b i \times M) \times f r e q j + b i \times M \times \sum 1 \leq j \leq L f r e q j b i \times M \times \sum 1 \leq j \leq L f r e q j, 这 是 一 个 常 量 ， 无 论 s 取 什 么 值 ， 它 都 不 变 。 因 为 a 为 正 整 数 ， 那 么 (a i \times d e g 2 j - b i \times M) 随 j 增 大 ， d e g 2 j 增 大 ， 所 以 我 们 仅 需 要 将 所 有 (a i \times d e g 2 j - b i \times M) 小 于 0 的 项 找 出 ， （ 为 正 值 就 开 始 加 了 ） 即 找 出 最 大 的 j 满 足 即 可 (54) (55) (56) (57) (58) (59) (60)

$\begin{align} Ei &=a_i\times \displaystyle\sum_{1\leq{j}\leq{L},deg_j\leq{s}} (deg_j^2 \times freq_j) + b_i\times \displaystyle\sum_{1\leq{j}\leq{L},deg_j>s} (M \times freq_j)\\ &=\displaystyle\sum_{1\leq{j}\leq{L},deg_j\leq{s}} ( a_i\times deg_j^2\times freq_j-b_i\times M\times freq_j)+ b_i\times M \times\displaystyle\sum_{1\leq{j}\leq{L}} freq_j\\ &= \displaystyle\sum_{1\leq{j}\leq{L},deg_j\leq{s}} (a_i\times deg_j^2-b_i\times M)\times freq_j+ b_i\times M \times\displaystyle\sum_{1\leq{j}\leq{L}} freq_j\\ &\\ &b_i\times M \times\displaystyle\sum_{1\leq{j}\leq{L}} freq_j ,这是一个常量，无论s取什么值，它都不变。\\ &因为a为正整数，那么(a_i\times deg_j^2-b_i\times M)随j增大，deg_j^2增大，所以我们仅需要将\\ &所有(a_i\times deg_j^2-b_i\times M)小于0的项找出，（为正值就开始加了）即找出最大的j满足即可\\ \end{align}$

解法：已经很明显了，预处理求出所有前i个deg²*freq的和、前i个freq的和。
     二分找出最大j的取值t，使得（a*deg[j]²<b*M）,然后答案就能以O(1)立刻求出。
     记住不能按推导公式直接求解： 因为b*M*sum_freq[L]超了LTM，只能求代入原公式
     求出b*M*（sum_freq[L]-sum_freq[t]）.

     二分也可用upper_bound()代替。

#include <iostream>
#include <fstream>
#include <stdio.h>
#include <string.h>
#include <stack>
#include <queue>
#include <vector>
#include <map>
#include <algorithm>
#define llt unsigned long long int
using namespace std;

const int Size=2*1e5+7;
llt d[Size],f[201000],Sd[201000],Sf[201000];

llt M,L;
int binfen(llt x,llt a){
    if(a*d[1]>x) return -1;
    int l=1,r=L;
    while(l<r){
        int mid=(l+r+1)>>1;
        if(a*d[mid]<x) l=mid;
        else r=mid-1;
    }

    return l;
}
int main(){

    scanf("%llu%llu",&M,&L);
    llt a,b;
    Sd[0]=Sf[0]=0;
    for(int i=1;i<=L;++i){
        scanf("%llu%llu",&a,&f[i]);
        d[i]=a*a;
        Sd[i]=Sd[i-1]+d[i]*f[i];
        Sf[i]=Sf[i-1]+f[i];
    }

    int Q;
    scanf("%d",&Q);
    while(Q--){
        llt a,b;
        scanf("%llu%llu",&a,&b);
        //cout<<ans<<endl;
        int t=binfen(b*M,a);
        if(t==-1)printf("%llu\n",b*M*(Sf[L]));
        else printf("%llu\n",a*(Sd[t])+b*M*(Sf[L]-Sf[t]));
    }

    return 0;
}