[bzoj4454]C Language Practice 解题报告

这题标算据说是O(n)-O(1)的gcd，然后一篇来自东欧的论文里说n必有拆分$n=x_1x_2x_3$,$x_i$要么是质数要么小于等于$\sqrt n$，这个我就完全没懂为什么。。然后打表出$(i,j) \forall i,j \in [0,\lfloor \sqrt n \rfloor ]$若求(n,m),就$x_1,x_2,x_3$依次与m求gcd。
但是这个拆分。。我感觉好像非常不科学的样子。
不过倒是看到另一个看起来有点靠谱的结论。。就是n的质因子个数为$O({\log n \over \log \log n})$ .所以这题我们可以分解质因子求gcd，设$S=10^6$时间复杂度$O(S+T((n+m)\log S+nm{\log S \over \log \log S}))≈2.4*10^9$。。实际跑起来的话，打表发现10^6以内一个数的平均质因子个数只有2…也就是说如果是随机数据的话，比上面那个看不懂的常数3的标算快到不知道哪里去了。。
代码：

#include<cstdio>
#include<iostream>
using namespace std;
#include<cmath>
#include<cstring>
int prime[100005][21];
int ind[1000005];
int smp[1000005];
void in(int &x){
    char c=getchar();
    while(c<'0'||c>'9')c=getchar();
    for(x=0;c>='0'&&c<='9';c=getchar())x=x*10+(c^'0');
}

struct NS{
    int data;
    int tot;
    pair<int,int> dvs[7];
}a[2005],b[2005];
void in(NS & x){
    in(x.data);
    x.tot=0;
    if(!x.data)return;
    //printf("%d:",x);
    for(int tmp=x.data;tmp!=1;++x.tot){
        x.dvs[x.tot]=make_pair(smp[tmp],0);
        for(;smp[tmp]==x.dvs[x.tot].first;tmp/=prime[smp[tmp]][1])++x.dvs[x.tot].second;
        //printf("(%d,%d) ",prime[x.dvs[x.tot].first][1],x.dvs[x.tot].second);
    }
    //puts("");
}
int main(){
    freopen("bzoj_4454.in","r",stdin);
    int T,n,m,i,j,ptot=0;
    //init
    memset(smp,-1,sizeof(smp));
    for(i=2;i<=1000000;++i){
        if(smp[i]==-1){
            prime[ptot][0]=ind[i]=1;
            prime[ptot][1]=i;
            smp[i]=ptot++;
        }
        if(ind[i])
            for(j=0;j<ptot&&i*prime[j][1]<=1000000;++j){
                smp[i*prime[j][1]]=j;
                if(i%prime[j][1]==0){
                    prime[j][ind[i*prime[j][1]]=ind[i]+1]=i*prime[j][1];
                    break;
                }
            }
        else
            for(j=0;j<ptot&&i*prime[j][1]<=1000000;++j){
                smp[i*prime[j][1]]=j;
                if(i%prime[j][1]==0)break;
            }
    }
    //work
    in(T);
    unsigned ans;
    int k,o,gcd;
    while(T--){
        in(n),in(m);
        for(i=0;i<n;++i)in(a[i]);
        for(j=0;j<m;++j)in(b[j]);
        ans=0;
        for(i=n;i--;)
            for(j=m;j--;)
                if(a[i].data==0)ans+=b[j].data^i^j;
                else if(b[j].data==0)ans+=a[i].data^i^j;
                else{
                    //printf("query(%d,%d):",a[i].data,b[j].data);
                    gcd=1;
                    for(k=o=0;k<a[i].tot&&o<b[j].tot;)
                        if(a[i].dvs[k].first<b[j].dvs[o].first)++k;
                        else if(a[i].dvs[k].first>b[j].dvs[o].first)++o;
                        else{
                            gcd*=prime[a[i].dvs[k].first][min(a[i].dvs[k].second,b[j].dvs[o].second)];
                            //printf("(%d,min(%d,%d)=%d)=%d ",a[i].dvs[k].first,(a[i].dvs[o].second,b[j].dvs[o].second),prime[a[i].dvs[k].first][min(a[i].dvs[o].second,b[j].dvs[o].second)]);
                            ++k,++o;
                        }
                    //printf("%d\n",gcd);
                    ans+=gcd^i^j;
                }
        printf("%u\n",ans);
    }
}

总结：
①n的最多质因子数是很少的，只有$\log n \over \log \log n$.
②1~n的平均质因子数可以视为常数。（$10^6$以内只有2，$10^7$以内只有3）