HDU1299--Diophantus of Alexandria HDU（109）

看到这题的时候想到了刘汝佳里面的原题，就直接暴力枚举一下。

#pragma warning(disable:4786)
#include<iostream>
#include<algorithm>
#include<math.h>
#include<stdio.h>
#include<time.h>
#include<stdlib.h>
#include<queue>
#include<set>
#include<vector>
#include<string>
#include<ctime>
#include<string.h>
using namespace std;
#define pi acos(-1.0)
//#define LL __int64
typedef long long LL;
#define INF 0x7fffffffffffffff
#define bug puts("hear!")
#define inf 0x7fffffff
#define eps 1e-10
#define FRE freopen("in.txt","r",stdin)
#define E exp(1.0)
#define mod 1000000007
int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}
int main(){
    int n,l;
    l=1;
    scanf("%d",&n);
    while(n--)
    {
        int a,q;
        q=0;
        scanf("%d",&a);
        for(int i=a;i<=2*a;i++)
        {
            int z,f,g;
            z=i/gcd(i,a)*a;
            f=z/a-z/i;
            g=z%f;
            if(f==1)
                q++;
        }
        printf("Scenario #%d:\n",l);
        printf("%d\n\n",q);
        l++;
    }
    return 0;
}

完全的TLE了。好吧，其实测试了一下99999999，发现果然过不了。虽然只是遍历了a~2*a但是好像超了范围。

然后试了试定义i为long long或者__int64，居然直接蹦了，你妹了。

索性看了看题解，发现大家代码基本上都一模一样（天下文章一大抄）。

看了看思路，发现原来不用这么做，直接用数论的素因子定理就好。于是明白了什么..............直接开工。

但是做到最后发现为什么q*3没明白，索性去查查

分解素因子时，如果n没法完全分解（n包含大于sqrt（n）素因子）sum*=3最后在sum=（sum+1）/2

在一个大牛的博客发现了这句话。

这就明白了，所以就A了...............

还是数论好玩。

#pragma warning(disable:4786)
#include<iostream>
#include<algorithm>
#include<cmath>
#include<stdio.h>
#include<time.h>
#include<stdlib.h>
#include<queue>
#include<set>
#include<vector>
#include<string>
#include<ctime>
#include<string.h>
using namespace std;
#define pi acos(-1.0)
//#define LL __int64
typedef long long LL;
#define INF 0x7fffffffffffffff
#define bug puts("hear!")
#define inf 0x7fffffff
#define eps 1e-10
#define FRE freopen("in.txt","r",stdin)
#define E exp(1.0)
#define mod 1000000007
const int N=40000;
int s[N],m[N],l;
void is_prime()
{
    int i,j;
    for (i=2;i<=N;i++)
    {
        for (j=2;i*j<N;j++)
        {
            s[i*j]=0;
        }
    }
    l=0;
    for (i=2;i<N;i++)
    {
        if(s[i])
            m[l++]=i;
    }
}
int main()
{
    int t,n,l,num;
    memset(s,1,sizeof(s));
    is_prime();
    while (scanf("%d",&t)!=EOF)
    {
        int z=1;
        while (t--)
        {
            scanf("%d",&n);
            int q=1;
            for(int i=0;i<l;i++)
            {
                int o=(int)sqrt(n*1.0)+1;
                if(m[i]>o)
                break;
                num=0;
                while (n%m[i]==0)
                {
                    num++;
                    n/=m[i];
                }
                q*=(2*num+1);
            }
            if(n>1)
                q*=3;
            printf("Scenario #%d:\n",z++);
            printf("%d\n\n",(q+1)/2);

        }
    }
    return 0;
}