题意:
有n*n的格子,问你从(0,0)能看到多少个点
从x点-y点上没有其他点就代表从x能看到y点
思路:
就是求gcd(x,y)==1的对数
遍历x,然后通过寻找与gcd(i,y)!=1的对数来得到答案
gcd(i,y)!=1 y<=n
对i进行素因子分解
i=p1^e1 * p2^e2 * p3^e3 * … pk^ek
Ans+=m/p1+m/p2+…m/pk-m/(p1*p2)-m/(p1*p3)….+m/(p1*p2*p3)
奇加偶减
pi必须是互异的质因子,如果不是互异的假设2和2*2,因为事件2完全包含事件2*2
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<map>
#include<set>
using namespace std;
#define lowbit(x) (x&(-x))
typedef long long LL;
const int maxn = 1000005;
const int inf=(1<<28)-1;
int Fac[50];
LL dfs(int pos,LL now,int step,int m)
{
LL Ans=0;
if(pos==Fac[0]+1)
{
if(now==1) return 0;
if(step&1)
Ans+=m/now;
else
Ans-=m/now;
return Ans;
}
Ans+=dfs(pos+1,now,step,m);
Ans+=dfs(pos+1,now*Fac[pos],step+1,m);
return Ans;
}
LL Solve(int x,int m)
{
Fac[0]=0;
for(int i=2;i*i<=x;++i)
if(x%i==0)
{
Fac[++Fac[0]]=i;
while(x%i==0)
{
x/=i;
}
}
if(x>1) Fac[++Fac[0]]=x;
return dfs(1,1,0,m);
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
int n,m;
scanf("%d%d",&n,&m);
if(n>m) swap(n,m);
LL Ans=m;
for(int i=2;i<=n;++i)
Ans+=m-Solve(i,m);
printf("%lld\n",Ans);
}
return 0;
}