牛客国庆集训派对Day1-C:Utawarerumono

https://www.nowcoder.com/acm/contest/201/C

题目描述

算术是为数不多的会让Kuon感到棘手的事情。通常她会找Haku帮忙,但是Haku已经被她派去买东西了。于是她向你寻求帮助。
给出一个关于变量x,y的不定方程ax+by=c,显然这个方程可能有多个整数解。Kuon想知道如果有解,使得p2*x2+p1*x+q2*y2+q1*y最小的一组整数解是什么。为了方便,你只需要输出p2*x2+p1*x+q2*y2+q1*y的最小值。

输入描述:

第一行三个空格隔开的整数a,b,c(0 ≤ a,b,c≤ 105)。
第二行两个空格隔开的整数p1,p2(1 ≤ p1,p2 ≤ 105)。
第三行两个空格隔开的整数q1,q2(1 ≤ q1,q2 ≤ 105)。

输出描述:

如果方程无整数解,输出“Kuon”。
如果有整数解,输出p2*x2+p1*x+q2*y2+q1*y的最小值。

示例1

输入

2 2 1
1 1
1 1

输出

Kuon

示例2

输入

1 2 3
1 1
1 1

输出

4

先利用GCD来判断不定方程ax+by=c是否有解,如果有解,将不定方程和带入方程p_{2}x^{2}+p_{1}x+q_{2}y^{2}+q_{1}y进行化简;

过程如下:

y=\dfrac {c-ax}{b}

将y带入方程可得:

p_{2}x^{2}+p_{1}x+q_{2},\dfrac {c^{2}+a^{2}x^{2}-2acx}{b^{2}}+q_{1}\dfrac {c-ax}{b}

=\left( p_{2}+\dfrac {a^{2}q_{2}}{b^{2}}\right) x^{2}+\left( p_{1}-\dfrac {2acq_{2}}{b^{2}}-\dfrac {aq_{1}}{b}\right) x+\dfrac {q_{2}c^{2}}{b^{2}}+\dfrac {q_{1}c}{b}

由抛物线的性质可知:方程p_{2}x^{2}+p_{1}x+q_{2}y^{2}+q_{1}y的最小值在对称轴附近

由上式化简可知,抛物线的对称轴为:\dfrac {\dfrac {q_{1}a}{b}+\dfrac {2acq_{2}}{b^{2}}-p_{1}}{2\left( p_{2}+\dfrac {a^{2}}{b^{2}}q_{2}\right) }

然后在这条抛物线两边来找满足ax+by=c的整数解就可以了(即满足\left( c-ax\right) \% b=0),然后求两边满足要求的最小解就可以了

#include<bits/stdc++.h>
#define ll long long
using namespace std;
int main()
{
    ll a,b,c;
    ll p1,p2,q1,q2,y;
    scanf("%lld%lld%lld%lld%lld%lld%lld",&a,&b,&c,&p1,&p2,&q1,&q2);
    ll gcd=__gcd(a,b);
    if(c%gcd) puts("Kuon");
    else{

        ll aa=(a*a*q2+b*b*p2);
        ll bb=b*b*p1-2*a*c*q2-a*q1*b;
        ll x=(-bb)/(2*aa);
		ll x1=x,x2=x;
        while((c-a*x1)%b)//并不是每个x都满足条件
            x1--;
        y=(c-a*x1)/b;
        ll ans1=p2*x1*x1+p1*x1+q2*y*y+q1*y;
        while((c-a*x2)%b)
            x2++;
        y=(c-a*x2)/b;
        ll ans2=p2*x2*x2+p1*x2+q2*y*y+q1*y;
        printf("%lld\n",min(ans1,ans2));
    }
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值