牛客xb月赛39 G 区间询问优化小技巧 欧拉筛 求区间最小质因子 树状数组求区间总数

题目

题解思路

1—n 有多少数可以表示为大于等于 K 的质数的乘积
也就是 求1-n 有多少个数，这个数的最小质因子大于等于K。
欧拉筛就是用数的最小质因子筛这个数为合数的，这样我们可以利用欧拉筛来筛出区间的每个数的最小质因子。
这样我们将1-n的每一个数的最小质因子当成树状数组的指针。
每多一个就将这个指针下的值++。
利用树状数组logn求前缀和的特性，求出比指针K大的区间总和数。

但是每次询问的区间大小不同。如果我们每次都清空再加入数，肯定会浪费很多时间复杂度。

可以得到一个性质区间变大会影响小区间答案，这样我们将区间从小到大排序记录编号，从小区间开始慢慢变大，记录小区间的答案并加入编号。

输出的时候从编号小的开始输出即可。

AC代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <vector>
#include <algorithm>
#include <map>
#include <string>
using namespace std;

const  int  INF =  0x3f3f3f3f;
const int N=3000010;

int primes[N],res;
int tre[N] ;
bool  st[N];
int fw[N] ;

struct node
{
    int x,y,p ;
}qu[N];

bool cmp ( node A , node B )
{
    return A.x < B.x ;
}
int lbit(int x )
{
    return x & ( -x );
}

void xg(int x , int y  )
{
    for (int i = x ; i <= 3000000 ; i += lbit(i) )
        tre[i] += y;
}
long long sum(int x  )
{
    long long ans = 0 ;
    for (int i = x ; i ; i -= lbit(i) )
        ans += tre[i];
    return ans;
}

void get_primes(int n)
{
    // LOOP1
    for(int i=2;i<=n;++i)
    {
        if(!st[i])
        {
            primes[res++]=i;
            fw[i] = i ;
        }
        // LOOP2
        for(int j=0;primes[j]<=n/i;++j)
        {
            st[primes[j]*i] = 1;
            fw[primes[j]*i] =  primes[j] ;
            if(i%primes[j]==0) break;
        }
    }
}


int main ()
{
    ios::sync_with_stdio(false);cin.tie(0);
    get_primes(3e6) ;
    int q ;
    cin >> q ;
    for ( int i = 1 ; i <= q ; i++ )
    {
        int t1 , t2 ;
        cin >> t1 >>  t2 ;
        qu[i] = {t1,t2,i} ;
    }
    sort(qu+1,qu+1+q,cmp) ;
    vector < pair <int , int > > ans ;
    qu[0] = {0,0,0} ;
    for ( int i = 1 ; i <= q ; i++ )
    {
        cout << qu[i].x << "\n" ;
        cout << fw[qu[i].x] << "\n" ;

        if ( qu[i].x > qu[i-1].x )
        {
            int t1 = qu[i-1].x  , t2 = qu[i].x ;
            for ( int j = t1 + 1 ; j <= t2 ; j++ )
                xg(fw[j],1) ;
            int t3 = sum(t2) - sum(qu[i].y-1);
            ans.push_back({i,t3}) ;
        }else
        {
            int t3 = sum(qu[i].x) - sum(qu[i].y-1) ;
            ans.push_back({i,t3}) ;
        }
    }
    sort(ans.begin(),ans.end()) ;
    for ( int i = 0 ; i < ans.size() ; i++ )
    {
        cout << ans[i].second << "\n" ;
    }
    return 0 ;
}