2015 Multi-University Training Contest 3 1002 RGCDQ

本文介绍了一种算法,用于解决求解给定区间[L,R]内任意两数的素因子个数的最大公约数的最大值的问题。通过预处理得到每个数的素因子个数,并利用前缀和的方法加速区间查询。

RGCDQ 

Problem's Link:  http://acm.hdu.edu.cn/showproblem.php?pid=5317


 

Mean: 

定义函数f(x)表示:x的不同素因子个数。

如:f(2) = 1, f(6) = 2;

给定L和R(L<=i<j<=R),求区间内任意不相等的两个数f(x)的最大公约数的最大值。

analyse:

因为2*3*5*7*11*13*17 >1e6,所以f(x)的值最大为7;

我们先打表求出每个数的f(x)值;

f[i][j]表示2~i中质因数个数为j的个数。

然后再利用前缀和f[r][i] - f[l-1][i],求出区间[l, r]的值。

Time complexity: O(NlogN)

 

Source code: 

/*
* this code is made by crazyacking
* Verdict: Accepted
* Submission Date: 2015-07-28-19.01
* Time: 0MS
* Memory: 137KB
*/
#include <queue>
#include <cstdio>
#include <set>
#include <string>
#include <stack>
#include <cmath>
#include <climits>
#include <map>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
#define  LL long long
#define  ULL unsigned long long
using namespace std;

int T , a , b;
int f [ 1000009 ][ 10 ];
const int NN = 1000005;
int p [NN ];
bool v [NN ];
void cntPrimeNum()
{
      for( int i = 2; i < NN; ++ i )
      {
            if( ! v [ i ] )
                  for( int j = i; j < NN; j += i )
                  {
                        v [ j ] = true;
                        ++p [ j ];
                  }
      }
}

int main()
{
      cntPrimeNum();
      for( int i = 1; i < NN; ++ i )
            for( int j = 7; j > 0; -- j )
                  if( p [ i ] == j ) f [ i ][ j ] = 1;
      for( int i = 1; i <= 1000000; ++ i )
            for( int j = 1; j <= 7; ++ j )
                  f [ i ][ j ] += f [ i - 1 ][ j ];
      scanf( "%d" , & T );
      while( T -- )
      {
            int ans = 1;
            scanf( "%d %d" , & a , &b );
            for( int i = 7; i > 0; -- i )
            {
                  if( f [b ][ i ] - f [ a - 1 ][ i ] > 1 )
                  {
                        ans = i;
                        break;
                  }
            }
            printf( "%d \n " , ans );
      }
      return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值