USACO3.2.4 Feed Ratios (ratios)

本文介绍了一种使用行列式求解最简配比的方法,包括如何通过给定的参数计算a, b, c, t,并判断是否有解。详细步骤和示例代码帮助读者理解和实现这一数学问题的应用。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

其实就是求一个最简配比a,b,c,t使a(x1,y1,z1)+b(x2,y2,z2)+c(x3,y3,z3)=t(k1,k2,k3)
用行列式求出a,b,c,t。如果t=0或abct不同号则说明无解。
否则分别除以gcd输出就行了。
/*
ID:shijiey1
PROG:ratios
LANG:C++
*/
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

int p[4], x[4], y[4], z[4];

int gcd(int a, int b) {
	while (b ^= a ^= b ^= a %= b);
	return a;
}

void scan() {
	for (int i = 1; i <= 3; i++) scanf("%d", &p[i]);
	for (int i = 1; i <= 3; i++) scanf("%d", &x[i]);
	for (int i = 1; i <= 3; i++) scanf("%d", &y[i]);
	for (int i = 1; i <= 3; i++) scanf("%d", &z[i]);
}

int main() {
	freopen("ratios.in", "r", stdin);
	freopen("ratios.out", "w", stdout);
	scan();
	int t = x[1] * y[2] * z[3] + x[3] * y[1] * z[2] + x[2] * y[3] * z[1]
		- x[3] * y[2] * z[1] - x[2] * y[1] * z[3] - x[1] * y[3] * z[2];
	if (!t) {
		printf("NONE\n");	
		return 0;
	}
	int a = p[1] * y[2] * z[3] + y[1] * z[2] * p[3] + z[1] * p[2] * y[3]
		- p[1] * z[2] * y[3] - p[2] * y[1] * z[3] - p[3] * y[2] * z[1];
	int b = x[1]*p[2]*z[3] + p[1]*z[2]*x[3] + z[1]*x[2]*p[3]
		-x[1]*z[2]*p[3] - p[1]*x[2]*z[3] - z[1]*p[2]*x[3];
	int c = x[1] * y[2] * p[3] + y[1] * p[2] * x[3] + y[3] * x[2] * p[1]
		-x[2] * y[1] * p[3] - y[2] * x[3] * p[1] - x[1] * y[3] * p[2];
	if (t < 0) {
		t = -t;
		a = -a;
		b = -b;
		c = -c;
	}
	if (a < 0 || b < 0 || c < 0) {
		printf("NONE\n");
		return 0;
	}
	int d = gcd(t, gcd(a, gcd(b, c)));
	printf("%d %d %d %d\n", a / d, b / d, c / d, t / d);
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值