P1072 Hankson 的趣味题

由于质因数只有log个，所以从这方面下手，

(1)由lcm(b0,x)=b1知

①x|b1，用b1直接限定x的质因数及其最大次数

②对于b0的每个质因数，若次数与b1相等，则x的次数任意取，若次数小于b1，则x的次数必为b1的次数。

(2)由gcd(a0,x)=a1知

①a1|x，以此限定最小次数Min2

②gcd(a0,x)=a1，a0的所有质因数的次数与x的次数取较小值等于a1的次数。即当前者大于Min2，后者Max须等于Min2。

然后就没了。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
#define rep(i,j,k) for(i=j;i<=k;++i)
#define per(i,j,k) for(i=j;i>=k;--i)
#define ll long long
#define db double
#define mkp(x,y) make_pair(x,y)
#define pii pair<int,int>
#define X first
#define Y second
const int S=50000;
struct YINSHU{
    int val,Max,Min,Min2;
}ys[30];
int a0,b0,a1,b1,cnt,pri[S],ln,ans;
bool vis[S];
int main(){
//	freopen("r.in","r",stdin);
//	freopen("w.out","w",stdout);
    int _,i,j,x;bool flg;
    for(i=2;i<S;++i){
        if(!vis[i])pri[++cnt]=i;
        for(j=1;j<=cnt&&i*pri[j]<S;++j){
            vis[i*pri[j]]=1;
            if(i%pri[j]==0)break;
        }
    }
    for(scanf("%d",&_);_--;){
        scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
        ln=0;
        rep(i,1,cnt){
            if(b1%pri[i]==0){
                for(x=0;b1%pri[i]==0;b1/=pri[i])++x;
                ys[++ln]=(YINSHU){pri[i],x,x,0};
            }
            if(b1==1)break;
        }
        if(b1>1)ys[++ln]=(YINSHU){b1,1,1,0};
        j=1;
        rep(i,1,cnt){
            if(b0%pri[i]==0){
                for(x=0;b0%pri[i]==0;b0/=pri[i])++x;
                while(ys[j].val!=pri[i])++j;
                if(ys[j].Max==x)ys[j].Min=0;
            }
            if(b0==1)break;
        }
        if(b0>1){
            while(ys[j].val!=b0)++j;
            if(ys[j].Max==1)ys[j].Min=0;
        }
        flg=0;j=1;
        rep(i,1,cnt){
            if(a1%pri[i]==0){
                for(x=0;a1%pri[i]==0;a1/=pri[i])++x;
                while(ys[j].val!=pri[i]&&j<=ln)++j;
                if(j>ln||ys[j].Max<x){
                    flg=1;break;
                }
                ys[j].Min2=x;
            }
            if(a1==1)break;
        }
        if(a1>1){
            while(j<=ln&&ys[j].val!=a1)++j;
            if(j>ln)flg=1;
            else ys[j].Min2=1;
        }
        if(flg){
            puts("0");continue;
        }
        j=1; 
        rep(i,1,cnt){
            if(a0%pri[i]==0){
                for(x=0;a0%pri[i]==0;a0/=pri[i])++x;
                while(j<=ln&&ys[j].val<pri[i])++j;
                if(j>ln)break;
                if(ys[j].val==pri[i]){
                    if(x>ys[j].Min2){
                        if(ys[j].Min>ys[j].Min2){
                            flg=1;break;
                        }
                        ys[j].Max=ys[j].Min2;
                    }
                }
            }
            if(a0==1)break;
        }
        if(a0>1){
            while(j<=ln&&ys[j].val<a0)++j;
            if(j<=ln&&ys[j].val==a0){
                if(!ys[j].Min2){
                    if(ys[j].Min)flg=1;
                    else ys[j].Max=0;
                }
            }
        }
        if(flg){
            puts("0");continue;
        }
        ans=1;
        rep(i,1,ln)ans*=ys[i].Max-max(ys[i].Min,ys[i].Min2)+1;
        printf("%d\n",ans);
    }
    return 0;
}

这道题爆了十几发OJ，我还纳闷零几年的NOIPT2怎么会要写一百多行，看完题解才发现是自己傻逼了：其实在线性筛的时候直接对每个质数处理就好了。

鈤狗，我觉得应该把这篇博客删了。