（素数暴力筛）HDU - 6069 Counting Divisors

分析：

现场的时候，学长A的很快，队友写了一个瞎几把优化的素数筛，C++T了。

然后本地试了15组最差数据，跑了6s，然后就用java重写一波，交了之后7sA了。。

又学了一招，，C++勉强超时java说不定可过。。

 

不过正解应该在《挑战程序设计》的120-121页，一模一样的数据规模，，学长就是改这个A的。。

据说标程也差不多。

 

代码（java，和C++一模一样的）：

import java.io.*;  
import java.math.*;  
import java.util.*;  
import java.text.*;  

public class Main {
    static int p[] = new int[100010];
    static long m[] = new long[1000010];
    static long m2[] = new long[1000010];
    static int num=0;
    static long mod = 998244353;
    
    static void init(){
        for(long i=2;i<=1000000;i++){
            if(m[(int)i]==0){
                p[num++]=(int)i;
                for(long j=i*i;j<=1000000;j+=i){
//                    System.out.println(j);
                    m[(int)j]=1;
                }
            }
        }
    }
    
    static public void main(String args[]){
        init();
        int t;
        Scanner cin = new Scanner (new BufferedInputStream(System.in));
        t=cin.nextInt();
        long l,r,k;
        while(t-->0){
            long sum=0;
            l=cin.nextLong();
            r=cin.nextLong();
            k=cin.nextLong();
            for(int i=0;i<=(int)(r-l);i++){
                m2[i]=m[i]=1;
            }
            for(int i=0;i<num;i++){
                long x=l/p[i];
                if(l%p[i]!=0)x++;
                long a;
                x*=p[i];
                while(x<=r){
                    a=x;
                    long c=0;
                    while(a%p[i]==0){
                        c++;
                        a/=p[i];
                    }
                    m[(int)(x-l)]=m[(int)(x-l)]*(c*k+1)%mod;
                    m2[(int)(x-l)]*=(x/a);
                    x+=p[i];
                }
            }
            for(long i=l;i<=r;i++){
                if(m2[(int)(i-l)]!=i)m[(int)(i-l)]=m[(int)(i-l)]*(k+1)%mod;
            }
            for(long i=l;i<=r;i++){
                sum=(sum+m[(int)(i-l)])%mod;
            }
            System.out.println(sum);
        }
        
        
    }
}