適切な栄養

Vasyaにはn個のバールがあります。Ber-Colaのボトル1本はバール、バーバー1本はbバールです。彼は、負でない整数のBer-Colaボトルと、負でない整数のBarsバーを購入できます。

Ber-ColaとBarsのバーのボトルをいくらか購入して、正確にnバールを使うことが可能かどうかを調べてください。

言い換えると、VasyaがBer-Colaとy Barsバーのxボトルを購入でき、x・a + y・b = nであるか、不可能であると言うことができるように、2つの非負の整数xとyを見つける必要があります。

入力
最初の行には、単一の整数n(1≤n≤10000000)—Vasyaが持っている金額が含まれています。

2行目には、単一の整数a(1≤a≤10000000)が含まれています—Ber-Colaのボトル1本のコスト。

3行目には単一の整数b(1≤b≤10000000)が含まれています—1バーバーのコスト。

出力
Vasyaが正確にnバールを費やすような方法でBarsとBer-Colaを購入できない場合は、«NO»(引用符なし)を印刷します。

それ以外の場合は、最初の行に«YES»(引用符なし)を出力します。2行目に、2つの非負の整数xとyを印刷します—正確にnバールを使うために、Vasyaが購入する必要のあるBer-Colaのボトルの数とBarsバーの数、つまりx・a + y・b = n。複数の回答がある場合は、それらのいずれかを印刷してください。

数値xとyのいずれも0に等しくすることができます。

例としては、
入力
7
2
3
出力
YES
2 1
入力
100
25
10
出力
YES
0 10
入力
15
4
8
出力
NOの
入力
9960594
2551
2557
出力
YES
1951 1949
注意
最初の例でVasyaは、BER-コーラ二ボトルと一つのバーのバーを購入することができます。彼は正確に2・2 + 1・3 = 7バールを費やします。

2番目の例では、Vasyaは正確にnバールを複数の方法で使用できます。

Ber-Colaのボトル2本とBarsバー5本を購入します。
Ber-Colaのボトルを4本購入し、Barsバーは購入しないでください。
Ber-Colaを購入したり、10Barsバーを購入したりしないでください。
3番目の例では、正確にnバールを費やすために、Ber-ColaバーとBarsバーを使用することは不可能です。

ax + by = cax + by = c
set d = gcd(a、b)d = gcd(a、b)の場合、d | cd | cが
次のように解決される場合にのみ、解決策があります。

拡張ユークリッドを使用して、ax0 + by0 = dax0 + by0 = dの解を見つけます

この場合、a(x0 ∗ c / d)+ b(y0 ∗ c / d)= ca(x0 ∗ c / d)+ b(y0 ∗ c / d)= cである
ため、特別な解はx '= x0 ∗ c / dになります。 、y '= y0 ∗ c / dx' = x0 ∗ c / d、y '= y0 ∗ c / d
および一般解=特殊解+均一解

均一解は、方程式ax + by = 0 ax + by = 0の解です。

したがって、一般的な解はx = x '+ k ∗ b / d、y = y'−k ∗ a /dk∈zです。

xの場合、最初にその倍数を乗算し、次にxが残りをb / dとして取得し、xが残りを取得してxを取得します。次に、この残りの場合、yを取得でき、現在のyが正であるかどうかを判断します。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <set>
#include <string>
#include <queue>
#include <map>
#include <stack>
#include <map>
#include <unordered_map>
#include <vector>
#include <cmath>
#include <ext/rope>
#include <bits/stdc++.h> 

using namespace std;

#define gt(x) x = read()
#define int long long
#define ios ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)

typedef pair<int, int> PII;
typedef unsigned long long ULL;

inline int read(int out = 0)
{
    
    
    char c;
    while((c=getchar()) < 48 || c > 57);
    while(c >= 48 && c <= 57) out=out*10+c-48,c=getchar();
    return out; 
}

const int N = 1010;
const int M = 35;
const int mod = 1e9 + 7;
const int PP = 131;
const double eps = 1e-10;

int n;
int a, b;
int x, y;
int exgcd(int a, int b, int &x, int &y)
{
    
    
    if (!b)
    {
    
    
        x = 1, y = 0;
        return a;
    }
    int d = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}

signed main(){
    
    
	gt(n), gt(a), gt(b);
	
	int d = exgcd(a,b, x, y);
	
//	cout << d << endl;
	if (n % d)    cout << "NO" << endl;
	else{
    
    
		x *= (n / d);
		int mod = (b / d);
		x= (x % mod + mod) % mod;
		y = (n - x * a) / b;
		if (x < 0 || y < 0){
    
    
			cout << "NO" << endl;
		}
		else{
    
    
			cout << "YES" << endl;
			cout << x << " " << y << endl;
		}
	}
	
	return 0;
}

おすすめ

転載: blog.csdn.net/qq_45772483/article/details/112790262