UVA 10277 Boastin' Red Socks

最新推荐文章于 2014-09-14 14:04:13 发布

原创最新推荐文章于 2014-09-14 14:04:13 发布 · 540 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

UVA

当前文章被收录于：

UVA

54 篇文章 0 人学习

订阅专栏查看详情

当前文章被以下社区和专栏收录：

本文介绍了一种高效算法，用于解决特定数学问题。通过引入gcd函数和数学推导，实现快速计算并找到最优解。

#include<iostream>
#include<cmath>
#include<cstdio>
#define LL long long
using namespace std;
LL gcd(LL a,LL b) {return b==0?a:gcd(b,a%b);}
int main()
{
    LL p,q;
    while (scanf("%lld%lld",&p,&q)!=EOF)
    {
        if (p==0 && q==0) break;
        if (p==0) {puts("0 2");continue;}
        LL t=gcd(p,q);p/=t;q/=t;LL N;
        //printf("%lld %lld\n",p,q);
        for (N=2;N<=50000;N++)
        {
            if (N*(N-1)%q==0)
          {
              LL b=N*(N-1)/q*p;
              double thta=sqrt((double)(4*b+1));
              double ans=(1+thta)/2.0;
          //    printf("%lf %lf %d\n",thta,ans,N);
              if (ans-floor(ans)<1e-7)
              {
                  printf("%lld %lld\n",(long long)ans,(long long)(N-ans));
                  break;
              }
          }
        }
        if (N>50000) puts("impossible");
    }
    return 0;
}

标签

#uva #数论