TC SRM 562 div2

250:  (x/k)*3 + (x/3) = p 给出P和K，求满足的x最小的值，直接从1到3*P枚举即可  

500:给定两个正数 N,M 1 <= A <= N, 1 <= B <= M ,  SSR(A, B) = (sqrt(A)+sqrt(B))^2  求满足SSR(A, B)为整数的可能数；

化简后可得 SSR(A, B) = A　＋ B + 2*sqrt(A)*sqrt(B);只要满足sqrt(A)*sqrt(B)为整数即可

总的来说就两种情况：

1：A和B都是可开方的数；

2：都不是可开方的数，但是化简后的根号里面的数相同： sqrt(8) = 2*sqrt（2）   和 sqrt(2)

这里遇到了一个特别操蛋的问题，就是在计算M里面能够达到多少倍的sqrt（x） 时，用了sqrt(a)/sqrt(b)本地测试正确，提交后在系统上就是不对，无语。后来慢慢找出改成sqrt（a/b）就对了。

#line 5 "TheSquareRootDilemma.cpp"

#include <cstdlib>

#include <cctype>

#include <cstring>

#include <cstdio>

#include <cmath>

#include <algorithm>

#include <vector>

#include <string>

#include <iostream>

#include <sstream>

#include <map>

#include <set>

#include <queue>

#include <stack>

#include <fstream>

#include <numeric>

#include <iomanip>

#include <bitset>

#include <list>

#include <stdexcept>

#include <functional>

#include <utility>

#include <ctime>

using namespace std;



#define CL(a,num) memset((a),(num),sizeof(a))

#define iabs(x)  ((x) > 0 ? (x) : -(x))

#define inf 0x7f7f7f7f

#define MOD 1073741824

#define lc l,m,rt<<1

#define rc m + 1,r,rt<<1|1



#define PB push_back

#define MP make_pair

#define REP(i,n) for(i=0;i<(n);++i)

#define FOR(i,l,h) for(i=(l);i<=(h);++i)

#define FORD(i,h,l) for(i=(h);i>=(l);--i)

typedef vector<int> VI;

typedef vector<string> VS;

typedef vector<double> VD;

typedef long long LL;

typedef pair<int,int> PII;



class TheSquareRootDilemma

{

        public:

        int countPairs(int N, int M)

        {



            int i;

            int qnum = (int)sqrt(1.0*M);//计算出M中有多少个可开方数

            int sum = 0;



            for (i = 1; i <= N; ++i)

            {

                double a = sqrt(1.0*i);

                if (a == (int)a)

                {

                    sum += qnum;

                }

                else

                {

                    int tmp = i;

                    double bi = sqrt(1.0*tmp);

                    int biNum = (int)bi;

                    //对该数化简



                    for (int j = 2; j <= biNum; ++j)

                    {

                        while (tmp % (j*j) == 0)

                        {

                            tmp /= (j*j);

                        }

                    }

                    //计算M中能够表示最大的几倍的sqrt(tmp)

                    double b = sqrt(1.0*M)/sqrt(1.0*tmp);

                    int bnum = (int)b;

                    if (bnum >= 1){

                        sum += bnum;

                    }

                }

            }

            return sum;

        }



};