#include #include #include using namespace std; #define N 1000010 long long h[N],g[N],ans; int n,m,prime[N],miu[N],idx; bool notprime[N]; int main() { miu[1]=1; for(int i=2;im) swap(n,m); memset(g,0,sizeof g),memset(h,0,sizeof h); for(int i=1;i<=n;i++) for(int j=i;j<=n;j+=i) g[i]+=miu[j]; for(int i=1;i<=m;i++) for(int j=i;j<=m;j+=i) h[i]+=miu[j]; for(int i=1;i<=n;i++) for(int j=i;j<=n;j+=i) ans+=miu[i]*miu[j/i]*g[j]*h[j]; printf("%I64d\n",ans); } }