【0521模拟赛】小Z爱数学

题目描述

小Z想求F(n,k)，F(n,k)表示n的所有因数pi中，满足n/pi <= k 的和。

小Z发现还是很水，所以他决定加大难度。

求

小Z还准备了很多个询问。现在你来解决一下吧。

输入输出格式

输入格式：

第一行两个整数m 表示询问的个数

接下来每行两个数ni,ki，表示这个询问的n和k

输出格式：

，

也就是把刚好等于k的答案加进去  然后k变大

对于每个询问输出一行一个整数，表示对应的答案。

题目理解：F（n，k）表示n的所有因数qi中，满足n/qi<=k的和，

m组，每组给定n和k，求        ；

m<=500000,ni<=100000

解析：k比较大的的答案  包含了k比较小的答案，所以要先把每次询问按k来排序，

然后我们每次回答询问的时候  如果k不同  就把k变大成现在要问的k，

         也就是把刚好等于k的答案插入线段树， 然后k变大

复杂度 O (nlog^2n+mlogn)

代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#define N 100000
#define ll long long
using namespace std;
inline int read(){
    int x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}
ll t[N*],ans[N*+];
int m;
struct que{int n,k,id;}e[N*+];
bool cmp(que i,que j){return i.k<j.k;}
ll query(int l,int r){
    ll sum=;
    l+=N-,r+=N-;
    sum+=t[l]+t[r];
    if(l==r) return sum-=t[l];
    for(;l^r^;l>>=,r>>=){
        if(~l&) sum+=t[l^];
        if( r&) sum+=t[r^];
    }
    return sum;
}
void Change(int x,int v){
    t[x+=N-]+=v;
    while(x) t[x>>=]=t[x<<]+t[x<<|];
}
void sol(int x){
    for(int i=;i*x<=N;i++) Change(i*x,i);
}
int main(){
    m=read();
    for(int i=;i<=m;i++){e[i].n=read();e[i].k=read();e[i].id=i;}
    sort(e+,e+m+,cmp);
    int j=;
    for(int i=;i<=m;i++){
        while(j<e[i].k) sol(++j);
        ans[e[i].id]=query(,e[i].n);
    }
    for(int i=;i<=m;i++) printf("%lld\n",ans[i]);
    return ;
}