洛谷 1072 Hankson 的趣味题——质因数界限讨论

题目：https://www.luogu.org/problemnew/show/P1072

思路是把每个数质因数分解，答案对于每个质因数的次数有选择的区间，通过这个计算。

指数的限制就是上限是b1，下限是a1；a0-a1后有剩余的自己不能有；b1-b0有剩余的自己不能剩（即必须满上限）。

分解质因数用了那个好像是 O( n^(1/4) ) 的方法。其实如果给的都是大质数是不是会被卡？

然后写了无比冗长的代码。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=;
int T,a0,a1,b0,b1,p[N],tot,hi[N],lo[N],hi2[N],lo2[N],lm[N],ans;
bool flag=;
int rdn()
{
  int ret=;bool fx=;char ch=getchar();
  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}
  while(ch>=''&&ch<='') ret=(ret<<)+(ret<<)+ch-'',ch=getchar();
  return fx?ret:-ret;
}
int main()
{
  T=rdn();
  while(T--)
    {
      a0=rdn(); a1=rdn(); b0=rdn(); b1=rdn();
      int n=b1; tot=;//limit
      for(int i=;i*i<=n;i++)
    {
      if(n%i==)
        {
          p[++tot]=i;lo[tot]=hi[tot]=lo2[tot]=hi2[tot]=;
          while(n%i==)n/=i,hi[tot]++;
          hi2[tot]=hi[tot];
        }
    }
      if(n>)
    {
      p[++tot]=n;lo[tot]=hi[tot]=lo2[tot]=hi2[tot]=;
      hi[tot]=hi2[tot]=;
    }

      flag=;
      n=a1; int p0=;//bottom
      for(int i=;i*i<=n;i++)
    if(n%i==)
      {
        int d=;
        while(n%i==)n/=i,d++;
        while(p0<tot&&p[p0]<i)p0++;
        if(p[p0]!=i){flag=;break;}
        lo[p0]=d;
      }
      if(flag){puts("");continue;}
      if(n>)
    {
      while(p0<tot&&p[p0]<n)p0++;
      if(p[p0]==n)
        lo[p0]=;
    }

      flag=;
      for(int i=;i<=tot;i++)if(hi[i]<lo[i]){flag=;break;}
      if(flag){puts("");continue;}

      n=b0; p0=;
      for(int i=;i*i<=n;i++)
    {
      if(n%i==)
        {
          int d=;
          while(n%i==)n/=i,d++;
          while(p0<tot&&p[p0]<i)p0++;
          if(p[p0]!=i)continue;
          lo2[p0]=d;//can't leave
        }
    }
      if(n>)
    {
      while(p0<tot&&p[p0]<n)p0++;
      if(p[p0]==n) lo2[p0]=;
    }
      for(int i=;i<=tot;i++)
    {
      int d=hi[i]-lo2[i];
      if(d) lo2[i]=hi[i]; else lo2[i]=;
    }

      n=a0; p0=;
      for(int i=;i*i<=n;i++)
    {
      if(n%i==)
        {
          int d=;
          while(n%i==)n/=i,d++;
          while(p0<tot&&p[p0]<i)p0++;
          if(p[p0]!=i)continue;
          if(d>lo[p0])hi2[p0]=lo[p0];//can't exist
        }
    }
      if(n>)
    {
      while(p0<tot&&p[p0]<n)p0++;
      if(p[p0]==n)
        if(>lo[p0])hi2[p0]=lo[p0];
    }

      for(int i=;i<=tot;i++)
    hi[i]=min(hi[i],hi2[i]),lo[i]=max(lo[i],lo2[i]);
      ans=; flag=;
      for(int i=;i<=tot;i++)
    {
      if(lo[i]>hi[i]){flag=;break;}
      ans*=hi[i]-lo[i]+;
    }
      printf("%d\n",flag?ans:);
    }
  return ;
}