hdu5255 魔法因子

　　题目地址:http://acm.hdu.edu.cn/showproblem.php?pid=5255

　　首先先预处理出一个p,使得p*因子X等于一个整数，且p最小，设q=p*X。

　　则题目则可以看成存在一个r，使得p*r和q*r，数字的首位和末位交换位置，而其他位置上的数字恰好不变。

　　那么怎么求这个r呢，考虑到abs(q*r-p*r)=abs(q-p)*r,且abs(q*r-p*r)=X99...99Y，X+Y=9，且中间数字全是9，由于只要统计10位数以内的答案，因此去枚举这种形式的数字，且判断q*r和p*r是否可行。

　　读入的时候没考虑到精度问题直接*10^6导致wa了10+次，因此最好用字符串读入。

　　代码:

#include<cstdio>
#include<algorithm>
#include<map>
#include<cstring>
#include<vector>
#include<cmath>
#include<string>
#define N 2000010
#define M 1000
using namespace std;
double n;
long long a,b,z,pp,qq,ans[N],P,p,q,o,k,tot,i,j,l;
long long gcd(long long a,long long b)
{
    if (b==) return a;
    return gcd(b,a%b);
}
int check(long long p,long long q)
{
    if (p==) return ;
    if (p>=P) return ;
    long long a,b,c=,pp;
    a=p%;
    pp=p;
    while (pp)
    {
        b=pp%;
        pp=pp/;
        c=c*;
    }
    c=c/;
    p=p-a+b-b*c+a*c;
    if (p==q) return ;return ;
}
int main()
{
    int test,ii;
    P=;
    P=P*;
    scanf("%d",&test);
    for (ii=;ii<=test;ii++)
    {
    tot=;
    scanf("%lf",&n);
    p=(n+0.0000005)*;
    q=p/gcd(p,);
    p=/gcd(p,);
    o=abs(p-q);

    b=;a=;
    for (k=;k<=;k++)
    {
    for (i=;i<=;i++)
    for (j=;j<=;j++)
    {
        z=b*i+a*+j;
        if (z%o==)
        {
            z=z/o;
            pp=p*z;qq=q*z;
            if (check(pp,qq))
            {
                tot++;
                ans[tot]=pp;
            }
        }
    }
    a=a*+;
    b=b*;
    }
    sort(ans+,ans++tot);
    j=;ans[]=-;
    for (i=;i<=tot;i++)
    if (ans[i]!=ans[i-])
    {
        j++;ans[j]=ans[i];
    }
    tot=j;
    printf("Case #%d:\n",ii);
    printf("%I64d\n",tot);
    for (i=;i<=tot-;i++)
    printf("%I64d ",ans[i]);
    if (tot)
    printf("%I64d\n",ans[tot]);
    }
}