2021CCPC广州 H. Three Integers-程序员宅基地

技术标签: c++  矩阵  线性代数  # 思维  

题目链接
You are given three non-negative integers a a a, b b b, and c c c. Find three positive integers x x x, y y y, and z z z that satisfy x   m o d   y = a x\bmod y=a xmody=a, y   m o d   z = b y\bmod z=b ymodz=b, and z   m o d   x = c z\bmod x=c zmodx=c.

Input

The first line contains an integer t t t ( 1 ≤ t ≤ 1 0 5 1\le t\le 10^5 1t105) — the number of test cases.

Each test case contains 3 3 3 integers a a a, b b b, c c c ( 0 ≤ a , b , c ≤ 1 0 9 0\le a,b,c\le 10^9 0a,b,c109) on a single line.

Output

For each test case, if there are such three integers satisfying the condition, output “YES”, then output the three integers x x x, y y y, z z z ( 1 ≤ x , y , z ≤ 1 0 18 1\le x,y,z\le 10^{18} 1x,y,z1018) on the following line, or “NO” otherwise.

Example

input

4
0 0 0
1 2 3
6 6 6
11 3 3

output

YES
1 1 1
YES
5 2 8
NO
YES
11 45 14

首先如果 a = b = c a=b=c a=b=c ,则当且仅当 a = b = c = 0 a=b=c=0 a=b=c=0 时有解,否则无解。
否则根据条件可知
{ x = k 1 y + a y = k 2 z + b z = k 3 x + c \left\{\begin{matrix} x=k_1y+a\\ y=k_2z+b\\ z=k_3x+c \end{matrix}\right. x=k1y+ay=k2z+bz=k3x+c
其中 k 1 , k 2 , k 3 ≥ 0 k_1,k_2,k_3\geq 0 k1,k2,k30
上式转化一下可以得到:
{ x = k 1 k 2 k 3 x + k 1 k 2 c + k 1 b + a y = k 2 k 3 x + k 2 c + b z = k 3 x + c \left\{\begin{matrix} x=k_1k_2k_3x+k_1k_2c+k_1b+a\\ y=k_2k_3x+k_2c+b\\ z=k_3x+c \end{matrix}\right. x=k1k2k3x+k1k2c+k1b+ay=k2k3x+k2c+bz=k3x+c
由于 x , y , z x,y,z x,y,z 均为正数且此时没有 a = b = c = 0 a=b=c=0 a=b=c=0 的情况,因此要使方程组中第一个式子成立,则 k 1 , k 2 , k 3 k_1,k_2,k_3 k1,k2,k3 中至少有一个为 0 0 0
这里讨论 k 1 , k 2 , k 3 k_1,k_2,k_3 k1,k2,k3 中有一个为 0 0 0 的情况:

  • k 1 = 0 k_1=0 k1=0 时,可将上式转化为:
    { x = a y = k 2 k 3 a + k 2 c + b z = k 3 a + c \left\{\begin{matrix} x=a\\ y=k_2k_3a+k_2c+b\\ z=k_3a+c \end{matrix}\right. x=ay=k2k3a+k2c+bz=k3a+c
    如果要满足题目条件,则
    { a > c k 3 a + c > b k 2 k 3 a + k 2 c + b > a \left\{\begin{matrix} a>c\\ k_3a+c>b\\ k_2k_3a+k_2c+b>a \end{matrix}\right. a>ck3a+c>bk2k3a+k2c+b>a
    为了确保上述不等式组恒成立,令
    { k 2 = 2 k 3 = ⌊ max ⁡ ( 0 , b − c ) a ⌋ + 1 \left\{\begin{matrix} k_2=2\\ k_3=\left \lfloor \frac{\max(0,b-c)}{a} \right \rfloor+1 \end{matrix}\right. { k2=2k3=amax(0,bc)+1
  • k 2 = 0 k_2=0 k2=0 时,可将上式转化为:
    { x = k 1 b + a y = b z = k 3 k 1 b + k 3 a + c \left\{\begin{matrix} x=k_1b+a\\ y=b\\ z=k_3k_1b+k_3a+c \end{matrix}\right. x=k1b+ay=bz=k3k1b+k3a+c
    如果要满足题目条件,则
    { b > a k 1 b + a > c k 3 k 1 b + k 3 a + c > b \left\{\begin{matrix} b>a\\ k_1b+a>c\\ k_3k_1b+k_3a+c>b \end{matrix}\right. b>ak1b+a>ck3k1b+k3a+c>b
    为了确保上述不等式组恒成立,令
    { k 3 = 2 k 1 = ⌊ max ⁡ ( 0 , c − a ) b ⌋ + 1 \left\{\begin{matrix} k_3=2\\ k_1=\left \lfloor \frac{\max(0,c-a)}{b} \right \rfloor+1 \end{matrix}\right. { k3=2k1=bmax(0,ca)+1
  • k 3 = 0 k_3=0 k3=0 时,可将上式转化为:
    { x = k 1 k 2 c + k 1 b + a y = k 2 c + b z = c \left\{\begin{matrix} x=k_1k_2c+k_1b+a\\ y=k_2c+b\\ z=c \end{matrix}\right. x=k1k2c+k1b+ay=k2c+bz=c
    如果要满足题目条件,则
    { c > b k 2 c + b > a k 1 k 2 c + k 1 b + a > c \left\{\begin{matrix} c>b\\ k_2c+b>a\\ k_1k_2c+k_1b+a>c \end{matrix}\right. c>bk2c+b>ak1k2c+k1b+a>c
    为了确保上述不等式组恒成立,令
    { k 1 = 2 k 2 = ⌊ max ⁡ ( 0 , a − b ) c ⌋ + 1 \left\{\begin{matrix} k_1=2\\ k_2=\left \lfloor \frac{\max(0,a-b)}{c} \right \rfloor+1 \end{matrix}\right. { k1=2k2=cmax(0,ab)+1

显然,上述三种情况的前提条件分别是 a > c a>c a>c b > a b>a b>a c > b c>b c>b ,这三种情况可以覆盖除 a , b , c a,b,c a,b,c 相等外的所有情况。因此针对每种情况将对应 k k k 值代入方程即可计算出 x , y , z x,y,z x,y,z 的其中一种方案。

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;

int main() {
    
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int T;
    cin >> T;
    while (T--) {
    
        int a, b, c;
        cin >> a >> b >> c;
        if (a == b && b == c) {
    
            if (!a) {
    
                cout << "YES\n";
                cout << 1 << " " << 1 << " " << 1 << "\n";
            } else cout << "NO\n";
            continue;
        }
        ll x, y, z, k;
        if (b > a) {
    
            k = max(0, (c - a) / b) + 1;
            x = k * b + a;
            y = b;
            z = (k * b + a) * 2 + c;
        } else if (a > c) {
    
            k = max(0, (b - c) / a) + 1;
            x = a;
            y = (k * a + c) * 2 + b;
            z = k * a + c;
        } else if (c > b) {
    
            k = max(0, (a - b) / c) + 1;
            x = (k * c + b) * 2 + a;
            y = k * c + b;
            z = c;
        } else
            assert(false);
        cout << "YES\n";
        cout << x << " " << y << " " << z << "\n";
        assert(x % y == a && y % z == b && z % x == c);
        assert(x <= 1e18);
        assert(y <= 1e18);
        assert(z <= 1e18);
    }
    return 0;
}
版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/qq_45323960/article/details/121368980

智能推荐

使用nginx解决浏览器跨域问题_nginx不停的xhr-程序员宅基地

文章浏览阅读1k次。通过使用ajax方法跨域请求是浏览器所不允许的,浏览器出于安全考虑是禁止的。警告信息如下:不过jQuery对跨域问题也有解决方案,使用jsonp的方式解决,方法如下:$.ajax({ async:false, url: 'http://www.mysite.com/demo.do', // 跨域URL ty..._nginx不停的xhr

在 Oracle 中配置 extproc 以访问 ST_Geometry-程序员宅基地

文章浏览阅读2k次。关于在 Oracle 中配置 extproc 以访问 ST_Geometry,也就是我们所说的 使用空间SQL 的方法,官方文档链接如下。http://desktop.arcgis.com/zh-cn/arcmap/latest/manage-data/gdbs-in-oracle/configure-oracle-extproc.htm其实简单总结一下,主要就分为以下几个步骤。..._extproc

Linux C++ gbk转为utf-8_linux c++ gbk->utf8-程序员宅基地

文章浏览阅读1.5w次。linux下没有上面的两个函数,需要使用函数 mbstowcs和wcstombsmbstowcs将多字节编码转换为宽字节编码wcstombs将宽字节编码转换为多字节编码这两个函数,转换过程中受到系统编码类型的影响,需要通过设置来设定转换前和转换后的编码类型。通过函数setlocale进行系统编码的设置。linux下输入命名locale -a查看系统支持的编码_linux c++ gbk->utf8

IMP-00009: 导出文件异常结束-程序员宅基地

文章浏览阅读750次。今天准备从生产库向测试库进行数据导入,结果在imp导入的时候遇到“ IMP-00009:导出文件异常结束” 错误,google一下,发现可能有如下原因导致imp的数据太大,没有写buffer和commit两个数据库字符集不同从低版本exp的dmp文件,向高版本imp导出的dmp文件出错传输dmp文件时,文件损坏解决办法:imp时指定..._imp-00009导出文件异常结束

python程序员需要深入掌握的技能_Python用数据说明程序员需要掌握的技能-程序员宅基地

文章浏览阅读143次。当下是一个大数据的时代,各个行业都离不开数据的支持。因此,网络爬虫就应运而生。网络爬虫当下最为火热的是Python,Python开发爬虫相对简单,而且功能库相当完善,力压众多开发语言。本次教程我们爬取前程无忧的招聘信息来分析Python程序员需要掌握那些编程技术。首先在谷歌浏览器打开前程无忧的首页,按F12打开浏览器的开发者工具。浏览器开发者工具是用于捕捉网站的请求信息,通过分析请求信息可以了解请..._初级python程序员能力要求

Spring @Service生成bean名称的规则(当类的名字是以两个或以上的大写字母开头的话,bean的名字会与类名保持一致)_@service beanname-程序员宅基地

文章浏览阅读7.6k次,点赞2次,收藏6次。@Service标注的bean,类名:ABDemoService查看源码后发现,原来是经过一个特殊处理:当类的名字是以两个或以上的大写字母开头的话,bean的名字会与类名保持一致public class AnnotationBeanNameGenerator implements BeanNameGenerator { private static final String C..._@service beanname

随便推点

二叉树的各种创建方法_二叉树的建立-程序员宅基地

文章浏览阅读6.9w次,点赞73次,收藏463次。1.前序创建#include&lt;stdio.h&gt;#include&lt;string.h&gt;#include&lt;stdlib.h&gt;#include&lt;malloc.h&gt;#include&lt;iostream&gt;#include&lt;stack&gt;#include&lt;queue&gt;using namespace std;typed_二叉树的建立

解决asp.net导出excel时中文文件名乱码_asp.net utf8 导出中文字符乱码-程序员宅基地

文章浏览阅读7.1k次。在Asp.net上使用Excel导出功能,如果文件名出现中文,便会以乱码视之。 解决方法: fileName = HttpUtility.UrlEncode(fileName, System.Text.Encoding.UTF8);_asp.net utf8 导出中文字符乱码

笔记-编译原理-实验一-词法分析器设计_对pl/0作以下修改扩充。增加单词-程序员宅基地

文章浏览阅读2.1k次,点赞4次,收藏23次。第一次实验 词法分析实验报告设计思想词法分析的主要任务是根据文法的词汇表以及对应约定的编码进行一定的识别,找出文件中所有的合法的单词,并给出一定的信息作为最后的结果,用于后续语法分析程序的使用;本实验针对 PL/0 语言 的文法、词汇表编写一个词法分析程序,对于每个单词根据词汇表输出: (单词种类, 单词的值) 二元对。词汇表:种别编码单词符号助记符0beginb..._对pl/0作以下修改扩充。增加单词

android adb shell 权限,android adb shell权限被拒绝-程序员宅基地

文章浏览阅读773次。我在使用adb.exe时遇到了麻烦.我想使用与bash相同的adb.exe shell提示符,所以我决定更改默认的bash二进制文件(当然二进制文件是交叉编译的,一切都很完美)更改bash二进制文件遵循以下顺序> adb remount> adb push bash / system / bin /> adb shell> cd / system / bin> chm..._adb shell mv 权限

投影仪-相机标定_相机-投影仪标定-程序员宅基地

文章浏览阅读6.8k次,点赞12次,收藏125次。1. 单目相机标定引言相机标定已经研究多年,标定的算法可以分为基于摄影测量的标定和自标定。其中,应用最为广泛的还是张正友标定法。这是一种简单灵活、高鲁棒性、低成本的相机标定算法。仅需要一台相机和一块平面标定板构建相机标定系统,在标定过程中,相机拍摄多个角度下(至少两个角度,推荐10~20个角度)的标定板图像(相机和标定板都可以移动),即可对相机的内外参数进行标定。下面介绍张氏标定法(以下也这么称呼)的原理。原理相机模型和单应矩阵相机标定,就是对相机的内外参数进行计算的过程,从而得到物体到图像的投影_相机-投影仪标定

Wayland架构、渲染、硬件支持-程序员宅基地

文章浏览阅读2.2k次。文章目录Wayland 架构Wayland 渲染Wayland的 硬件支持简 述: 翻译一篇关于和 wayland 有关的技术文章, 其英文标题为Wayland Architecture .Wayland 架构若是想要更好的理解 Wayland 架构及其与 X (X11 or X Window System) 结构;一种很好的方法是将事件从输入设备就开始跟踪, 查看期间所有的屏幕上出现的变化。这就是我们现在对 X 的理解。 内核是从一个输入设备中获取一个事件,并通过 evdev 输入_wayland

推荐文章

热门文章

相关标签