题目链接

题面：

题意：
v [ i ] [ j] = ( a [ i ] + b [ j ] ) %k

现在给定一个n*m的矩阵w，问w能不能由以上方式构造出来，若不能输出NO。
若能输出YES，并且输出k，a1–an，b1–bm

贴一下官方题解：

题解翻译：
我们发现，如果 ai ，bi 是其的一个解，那么对于任意一个整数p来说，ai - p 和 ai + p也是他的一个解。
所以我们可以将 a1 赋值为0，那么 bi 就可以通过 w 的第一行求出来。
然后通过 b1 和 w 的第一列将 ai 求出来，此时允许 ai 和 bi 有负值（最后取模加模即可）。

我们令 e [ i ] [ j ] = abs ( a [ i ] + b [ j ] - w [ i ] [ j ] )
若 e [ i ] [ j ] 全部为0，那么 k > max（ w [ i ] [ j ] ）即可。
若 e [ i ] [ j ] 并非全部为0，那么令 g = gcd （e [ i ] [ j ]），显然 g 是最大的 k（因为要使e矩阵全部变为0），若 g > max（ w [ i ] [ j ] ），那么 我们令 k 为 g，否则无解。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<string>
#include<queue>
#include<bitset>
#include<map>
#include<set>
#define ll long long
#define llu unsigned ll
#define ld long double
#define ui unsigned int
#define pr make_pair
#define pb push_back
#define ui unsigned int
#define lc (cnt<<1)
#define rc (cnt<<1|1)
#define len(x) (t[(x)].r-t[(x)].l+1)
#define tmid ((l+r)>>1)
#define forhead(x) for(int i=head[(x)];i;i=nt[i])
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)>(y)?(y):(x))
using namespace std;
const int inf=0x3f3f3f3f;
const ll lnf=0x3f3f3f3f3f3f3f3f;
const double dnf=1e18;
const int mod=1e9+7;
const double eps=1e-8;
const double pi=acos(-1.0);
const int maxn=110;
const int maxm=100100;
const int up=100000;
const int hashp=13331;

ll w[maxn][maxn],e[maxn][maxn],a[maxn],b[maxn];
ll gcd(ll a,ll b)
{
    if(a==0||b==0) return a|b;
    return gcd(b,a%b);
}

int main(void)
{
    ll n,m,k;
    scanf("%lld%lld",&n,&m);
    ll maxx=0;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
            scanf("%lld",&w[i][j]),maxx=max(maxx,w[i][j]);
    }
    for(int i=1;i<=m;i++)
        b[i]=w[1][i];
    for(int i=2;i<=n;i++)
        a[i]=w[i][1]-b[1];
    ll g=0;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
            e[i][j]=abs(a[i]+b[j]-w[i][j]),g=gcd(g,e[i][j]);
    }
    if(g==0)
        k=maxx+1;
    else
    {
        if(g<=maxx)
        {
            printf("NO\n");
            return 0;
        }
        k=g;
    }

    printf("YES\n");
    printf("%lld\n",k);
    for(int i=1;i<=n;i++)
        printf("%lld ",(a[i]%k+k)%k);
    putchar('\n');
    for(int i=1;i<=m;i++)
        printf("%lld ",(b[i]%k+k)%k);
    putchar('\n');
    return 0;
}