bzoj 2820

 

收获：

当一个东西的取值有限时，我们可以枚举它，然后统计它被计算了多少次。

 

#include <cstdio>
#include <iostream>
using namespace std;

typedef long long dnt;

int prm[10000010], isnot[10000010], mu[10000010], f[10000010], ptot;

void init( int n ) {
    mu[1] = 1;
    for( int i=2; i<=n; i++ ) {
        if( !isnot[i] ) {
            prm[++ptot] = i;
            mu[i] = -1;
        }
        for( int j=1; j<=ptot && i*prm[j]<=n; j++ ) {
            isnot[i*prm[j]]=true;
            if( i%prm[j]==0 ) {
                mu[i*prm[j]] = 0;
                break;
            }
            mu[i*prm[j]] = -mu[i];
        }
    }
    for( int i=1; i<=ptot; i++ ) {
        int p = prm[i];
        for( int j=p; j<=n; j+=p )
            f[j] += mu[j/p];
    }
    for( int i=1; i<=n; i++ )
        f[i] += f[i-1];
}
dnt calc( dnt n, dnt m ) {
    if( n>m ) swap(n,m);
    dnt rt = 0;
    for( int i=1; i<=n; i++ ) {
        int ii = min( n/(n/i), m/(m/i) );
        rt += (f[ii]-f[i-1])*(n/i)*(m/i);
        i = ii;
    }
    return rt;
}
int main() {
    int T;
    init( 10000000 );
    scanf( "%d", &T );
    while( T-- ) {
        int n, m;
        scanf( "%d%d", &n, &m );
        printf( "%lld\n", calc(n,m) );
    }
}