欧拉函数之HDU2588 GCD

N的所有因子中，求大于等于M的因子a的euler(N/a)之和。

原因，可以参考上一篇文章~HDU4983

#include <stdio.h>
#include <ctype.h>
#include <string.h>
#include <stdlib.h>
#include <limits.h>
#include <math.h>
#include <algorithm>
#include <stack>
#include <queue>
#include <vector>
#include <map>
#include <set>
#include <string>
#include <sstream>
using namespace std;

#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1

typedef long long LL;
const double pi=4.0*atan(1.0);
const int MAXN=105;
const int INF=1<<30;
const LL M=1000000007;
int euler(int n)
{
	LL ans=n;
	for(int i=2;(LL)i*i<=n;i++)
		if(n%i==0)
		{
			ans=ans/i*(i-1);
			while(n%i==0)
				n/=i;
		}
		if(n>1)
			ans=ans/n*(n-1);
		return ans;
}

int main()
{
	int n,m;
	int T;
	int i,t;
	int ans=0;
	while(scanf("%d",&T)!=EOF)
	{
		while(T--)
		{
			scanf("%d%d",&n,&m);
			ans=0;
			for(i=1;i*i<=n;i++)
			{
				if(n%i==0)
				{
					if(i*i==n&&i>=m)
						ans+=euler(i);
					else
					{
						t=n/i;
						if(i>=m)
							ans+=euler(t);
						if(t>=m)
							ans+=euler(i);
					}
				}
			}
			printf("%d\n",ans);
		}
	}
	return 0;
}