数学+dp HDOJ 5317 RGCDQ

题目传送门

 /*
     题意：给一个区间，问任意两个数的素数因子的GCD最大
     数学+dp：预处理出f[i]，发现f[i] <= 7，那么用dp[i][j] 记录前i个f[]个数为j的数有几个，
             dp[r][j] - dp[l-1][j]表示区间内j的个数，情况不多，分类讨论一下
 */
 #include <cstdio>
 #include <algorithm>
 #include <cstring>
 #include <vector>
 #include <cmath>
 using namespace std;

 const int MAXN = 1e6 + ;
 const int INF = 0x3f3f3f3f;
 bool is_prime[MAXN];
 int f[MAXN];
 int dp[MAXN][];
 int cnt[];
 int n;

 int GCD(int a, int b)    {
     return (b == ) ? a : GCD (b, a % b);
 }

 void solve(void)    {
     memset (is_prime, true, sizeof (is_prime));
     for (int i=; i<=n; ++i)    {
         if (is_prime[i])    {
             for (int j=i; j<=n; j+=i)    {
                 f[j]++;    is_prime[j] = false;
             }
         }
     }
 }

 void work(void)    {
     memset (dp, , sizeof (dp));
     for (int i=; i<=n; ++i)    {
         for (int j=; j<=; ++j)    {
             dp[i][j] = dp[i-][j];
         }
         dp[i][f[i]]++;
     }
 }

 int main(void)    {       //HDOJ 5317 RGCDQ
     //freopen ("B.in", "r", stdin);
     n = ;
     solve ();    work();

     int T;
     scanf ("%d", &T);
     while (T--)    {
         int l, r;
         scanf ("%d%d", &l, &r);
         bool ok = false;
         for (int i=; i>=; --i)    {
             if (dp[r][i] - dp[l-][i] >= )    {
                 printf ("%d\n", i);    ok = true;    break;
             }
         }
         if (!ok)    {
             if(dp[r][] - dp[l-][] ==  && dp[r][] - dp[l-][] == ) {
                 printf("2\n");
             } else if(dp[r][] - dp[l-][] ==  && dp[r][] - dp[l-][] == ) {
                 printf("2\n");
             } else if(dp[r][] - dp[l-][] ==  && dp[r][] - dp[l-][] == ) {
                 printf("2\n");
             } else {
                 printf("1\n");
             }
         }
     }

     return ;
 }