Codeforces Round #289 (Div. 2, ACM ICPC Rules) D. Restoring Numbers 构造 数学

18 篇文章 0 订阅
D. Restoring Numbers
time limit per test
1 second
memory limit per test
256 megabytes
input
standard input
output
standard output

Vasya had two arrays consisting of non-negative integers: a of size n and b of size m. Vasya chose a positive integer k and created ann × m matrix v using the following formula:

Vasya wrote down matrix v on a piece of paper and put it in the table.

A year later Vasya was cleaning his table when he found a piece of paper containing an n × m matrix w. He remembered making a matrix one day by the rules given above but he was not sure if he had found the paper with the matrix v from those days. Your task is to find out if the matrix w that you've found could have been obtained by following these rules and if it could, then for what numbersk, a1, a2, ..., an, b1, b2, ..., bm it is possible.

Input

The first line contains integers n and m (1 ≤ n, m ≤ 100), separated by a space — the number of rows and columns in the found matrix, respectively.

The i-th of the following lines contains numbers wi, 1, wi, 2, ..., wi, m (0 ≤ wi, j ≤ 109), separated by spaces — the elements of the i-th row of matrix w.

Output

If the matrix w could not have been obtained in the manner described above, print "NO" (without quotes) in the single line of output.

Otherwise, print four lines.

In the first line print "YES" (without quotes).

In the second line print an integer k (1 ≤ k ≤ 1018). Note that each element of table w should be in range between 0 and k - 1inclusively.

In the third line print n integers a1, a2, ..., an (0 ≤ ai ≤ 1018), separated by spaces.

In the fourth line print m integers b1, b2, ..., bm (0 ≤ bi ≤ 1018), separated by spaces.

Sample test(s)
input
2 3
1 2 3
2 3 4
output
YES
1000000007
0 1 
1 2 3 
input
2 2
1 2
2 0
output
YES
3
0 1 
1 2 
input
2 2
1 2
2 1
output
NO
Note

By  we denote the remainder of integer division of b by c.

It is guaranteed that if there exists some set of numbers k, a1, ..., an, b1, ..., bm, that you could use to make matrix w, then there also exists a set of numbers that meets the limits 1 ≤ k ≤ 10181 ≤ ai ≤ 10181 ≤ bi ≤ 1018 in the output format. In other words, these upper bounds are introduced only for checking convenience purposes.

题目给出公式aij = (aai + bbj)% g; 给出vij,要求是否存在这样的数列,若存在则求出aai bbj 和g

首先每行每列,通过 (a[i+1][j+1] - a[i+1][j]) - (a[i][j+1] - a[i][j]) =k1g - k2g = k3g ; (a[i+1][j+1] - a[i][j+1]) - (a[i+1][j] - a[i][j])) = k4g - k5g = k6g,这样求

一下gcd就可以得出g了,若g为0,说明没有冲突,直接置为最大值。若a有数比g大的数,则说明求出的g不能满足题意,直接输出no,

否则,说明不冲突,则要构造出a b即可!如何构造a b呢,直接令aa[0] = 0(因为最后都要取模的,每一个数是几并不影响后面的值),则bb[0] = a[0][0],通过a[i][0] - a[i-1][0]  = a[i-1]%g,即可构造出来,且构造出

来的一定成立,因为,若有一个数不成立,则这个数与上面的数相减后,就不会是g的倍数,这与前面是矛盾的。这样就ok了!


#define N 105
#define M 100005
#define maxn 205
#define MOD 1000000007
int n,m;
ll a[N][N],aa[N],bb[N];
int main()
{
    while(S2(n,m)!=EOF)
    {
       FI(n){
            FJ(m){
                cin>>a[i][j];
            }
       }
       ll g = 0;
       FI(n - 1){
            FJ(m - 1){
                ll g1 = abs((a[i+1][j+1] - a[i+1][j]) - (a[i][j+1] - a[i][j]));
                ll g2 = abs((a[i+1][j+1] - a[i][j+1]) - (a[i+1][j] - a[i][j]));
                g = gcd(g,gcd(g1,g2));
            }
       }
       if(g == 0) g = MOD;
       bool flag = true;
       for(int i=0;i<n && flag;i++){
            for(int j = 0;j<m && flag;j++){
                if(a[i][j] >= g) flag = false;
            }
       }
       if(!flag){printf("NO\n");continue;}
       aa[0] = 0;bb[0] = a[0][0];
       for(int i=1;i<n;i++){aa[i] = ( aa[i-1] + a[i][0] - a[i-1][0] + g) % g;}
       for(int i=1;i<m;i++){bb[i] = ( bb[i-1] + a[0][i] - a[0][i-1] + g) % g;}
       cout<<"YES\n"<<g<<endl;
       for(int i=0;i<n;i++){cout<<aa[i]<<" ";}cout<<endl;
       for(int i=0;i<m;i++){cout<<bb[i]<<" ";}cout<<endl;
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值