题意:给定一个序列,多次询问

将a数组从小到大排序,下面那个值只有不超过32种,于是预处理f[i][j],表示分母为i时,aj/i的前缀和是多少。

然后对于一个给定的p,一定将分母划分成了一些连续的段落,通过枚举这些分母,二分获得分母变化的位置,将区间和累计进答案。

注意,对于给定的p,一个分母控制的某段区间是a的元素属于p^i+1到p^(i+1)的这段。

复杂度log^2n。

#include<cstdio>
//#include<cmath>
#include<algorithm>
using namespace std;
int f,cc;
void R(int &x){
cc=0;f=1;
for(;cc<'0'||cc>'9';cc=getchar())if(cc=='-')f=-1;
for(x=0;cc>='0'&&cc<='9';cc=getchar())(x*=10)+=(cc-'0');
x*=f;
}
#define MOD 1000000000
typedef long long ll;
int T,n,m,a[35][500005];
ll c[500005];
//int poses[100005][35];
//double Log(const int &p,const int &x){
// return log((double)x)/log((double)p);
//}
int main(){
R(T);
for(;T;--T){
R(n); R(m);
for(int i=1;i<=n;++i){
R(a[1][i]);
}
sort(a[1]+1,a[1]+n+1);
for(int i=1;i<=n;++i){
c[i]=(ll)a[1][i];
}
for(int i=1;i<=n;++i){
for(int j=2;j<=32;++j){
a[j][i]=a[1][i]/j;
}
}
for(int i=1;i<=32;++i){
for(int j=1;j<=n;++j){
a[i][j]=(a[i][j]+a[i][j-1])%MOD;
}
}
// for(int p=2;p*p<=MOD;++p){
// for(int i=1;i<=32;++i){
// if((int)(ceil(Log(p,c[n]))+0.5)<i){
// poses[p][i]=n+1;
// continue;
// }
// int l=1,r=n;
// while(l<r){
// int mid=(l+r>>1);
// if((int)(ceil(Log(p,c[mid]))+0.5)>=i){
// r=mid;
// }
// else{
// l=mid+1;
// }
// }
// poses[p][i]=l;
// }
// }
int p;
ll ans=0;
int popo[35];
for(int zu=1;zu<=m;++zu){
R(p);
// if((ll)p*(ll)p<=(ll)MOD){
// for(int i=1;i<32;++i){
// if(poses[p][i]>n){
// break;
// }
// ans=(ans+((ll)zu*(ll)((a[i][poses[p][i+1]-1]-a[i][poses[p][i]-1]+MOD)%MOD))%(ll)MOD)%(ll)MOD;
// }
// }
// else{
ll pp=1;
for(int i=1;i<32;++i){
ll* be=upper_bound(c+1,c+n+1,pp);
if(be-c>n){
break;
}
pp*=(ll)p;
ans=(ans+((ll)zu*(ll)((a[i][upper_bound(c+1,c+n+1,pp)-c-1]-a[i][be-c-1]+MOD)%MOD))%(ll)MOD)%(ll)MOD;
}
// }
}
printf("%lld\n",ans);
}
return 0;
}

【二分】【预处理】zoj4029 Now Loading!!!的更多相关文章

  1. NOIP2015聪明的质检员[二分 | 预处理]

    背景 NOIP2011 day2 第二题 描述 小T 是一名质量监督员,最近负责检验一批矿产的质量.这批矿产共有 n 个矿石,从 1到n 逐一编号,每个矿石都有自己的重量 wi 以及价值vi .检验矿 ...

  2. Codeforces Round #321 (Div. 2) B 二分+预处理

    B. Kefa and Company time limit per test 2 seconds memory limit per test 256 megabytes input standard ...

  3. HDU 6274 二分+预处理(CCPC K题

    #include"bits/stdc++.h" #define db double #define ll long long #define vec vector<ll> ...

  4. [luoguP3092] [USACO13NOV]没有找零No Change(状压DP + 二分)

    传送门 先通过二分预处理出来,每个硬币在每个商品处最多能往后买多少个商品 直接状压DP即可 f[i]就为,所有比状态i少一个硬币j的状态所能达到的最远距离,在加上硬币j在当前位置所能达到的距离,所有的 ...

  5. EC R 86 D Multiple Testcases 构造 贪心 二分

    LINK:Multiple Testcases 得到很多种做法.其中O(n)的做法值得一提. 容易想到二分答案 check的时候发现不太清楚分配的策略. 需要先考虑如何分配 容易发现大的东西会对小的产 ...

  6. 「Poetize3」Heaven Cow与God Bull

    描述 Description 给定一个整数n,求一个整数m,满足m<=n,并且m/phi(m)的值最大.注:phi(m)代表m的欧拉函数,即不大于m且与m互质的数的个数. 题解:m/phi(m) ...

  7. D. 实验室传染病

    D. 实验室传染病 题意 给出 n 个人的位置,以及每个人的传染范围,当一个人患病时,他的传染范围内(包括边界上)的人全部会被感染并继续向外传播. 求以每个人为传染源最多有多少人被感染. 分析 首先二 ...

  8. 【倍增】Tak and Hotels II @ABC044&ARC060/upcexam6463

    6463: Tak and Hotels II 时间限制: 1 Sec  内存限制: 128 MB 题目描述 N hotels are located on a straight line. The ...

  9. CF #442 div2

    A 判断下5个名字出现了几次.pre数据巨弱,就这么一水题在std测刷掉了非常多的人.. /** @Date : 2017-10-24 16:04:41 * @FileName: A.cpp * @P ...

随机推荐

  1. Servlet笔记6--Servlet程序改进

    第一步改进,GenericServlet: 我们目前所有放入Servlet类直接实现了javax.servlet.Servlet接口,但是这个接口中有很多方法是目前不需要的,我们可能只需要编写serv ...

  2. Apple Notification Center Service--ANCS【转】

    Apple Notification Center Service 转自:http://studentdeng.github.io/blog/2014/03/22/ancs/ MAR 22ND, 20 ...

  3. Codeforces Round #504 E. Down or Right

    Codeforces Round #504 E. Down or Right 题目描述:交互题. 有一个\(n \times n\)的方阵,有一些格子是障碍,从\((1, 1)\)出发,只能向右向下走 ...

  4. Focal Loss笔记

    论文:<Focal Loss for Dense Object Detection> Focal Loss 是何恺明设计的为了解决one-stage目标检测在训练阶段前景类和背景类极度不均 ...

  5. [HBase]mem store flusher 流程

  6. 数据库-mysql安装

    MySQL 安装 所有平台的Mysql下载地址为: MySQL 下载. 挑选你需要的 MySQL Community Server 版本及对应的平台. Linux/UNIX上安装Mysql Linux ...

  7. Python 常用的内建函数

    内建函数 ​ Build-in Function,启动python解释器,输入dir(__builtins__), 可以看到很多python解释器启动后默认加载的属性和函数,这些函数称之为内建函数, ...

  8. java 证书体系及应用,自已做https证书

    原文: https://blog.csdn.net/wjq008/article/details/49071857 接下来我们将域名www.zlex.org定位到本机上.打开C:\Windows\Sy ...

  9. MySQL学习笔记:delete from与truncate table的区别

    在Mysql数据库的使用过程中,删除表数据可以通过以下2种方式: delete from table_name truncate table table_name (1)delete from语句可以 ...

  10. mybatis中多条件判断---choose when的用法

    <select id="getFunctionByPage" resultMap="FunctionRlt"> SELECT K.FUNCTION_ ...