[nowcoder5666D]Quadratic Form

首先猜测$\sum_{i=1}^{n}b_{i}x_{i}$取到最小值时存在$x_{i}$满足$\sum_{i=1}^{n}\sum_{j=1}^{n}A_{i,j}x_{i}x_{j}=1$，否则一定可以适当调整某一个$x_{i}$

因此可以使用拉格朗日乘数法，构造函数$F(x_{1},x_{2},...,x_{n},\lambda)=\sum_{i=1}^{n}b_{i}x_{i}+\lambda(\sum_{i=1}^{n}\sum_{j=1}^{n}A_{i,j}x_{i}x_{j}-1)$，要保证$x_{i}$和$\lambda$的偏导数为0，由此可以列出$n+1$个方程：$，\begin{cases}\forall 1\le i\le n，b_{i}+2\lambda\sum_{j=1}^{n}A_{i,j}x_{j}=0\\ \sum_{i=1}^{n}\sum_{j=1}^{n}A_{i,j}x_{i}x_{j}=1 \end{cases}$

把A、b和x看成矩阵来表示，即得到$\begin{cases}b+2\lambda Ax=[0\ 0...0]\\ x^{T}Ax=[1]\end{cases}$，分别化简即得$\begin{cases}x=-\frac{1}{2\lambda}A^{-1}b\\ x^{T}x=A^{-1}\end{cases}$，同时由于$bx^{T}=xb^{T}$，则答案$(xb^{T})^{2}=(bx^{T})(xb^{T})=b(x^{T}x)b^{T}=bA^{-1}b^{T}$，矩阵求逆+矩阵乘法即可，复杂度$o(n^{3})$

#include<bits/stdc++.h>
using namespace std;
#define N 205
#define mod 998244353
int n,ans,a[N][N],b[N][N],c[N],d[N];
int ksm(int n,int m){
    if (!m)return 1;
    int s=ksm(n,m>>1);
    s=1LL*s*s%mod;
    if (m&1)s=1LL*s*n%mod;
    return s;
}
int main(){
    while (scanf("%d",&n)!=EOF){
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)scanf("%d",&a[i][j]);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)b[i][j]=(i==j);
        for(int i=1;i<=n;i++){
            for(int j=i;j<=n;j++)
                if (a[j][i]){
                    for(int k=i;k<=n;k++)swap(a[i][k],a[j][k]);
                    break;
                }
            int t=ksm(a[i][i],mod-2);
            for(int j=1;j<=n;j++){
                a[i][j]=1LL*t*a[i][j]%mod;
                b[i][j]=1LL*t*b[i][j]%mod;
            }
            for(int j=i+1;j<=n;j++){
                int t=a[j][i];
                for(int k=1;k<=n;k++){
                    a[j][k]=(a[j][k]-1LL*t*a[i][k]%mod+mod)%mod;
                    b[j][k]=(b[j][k]-1LL*t*b[i][k]%mod+mod)%mod;
                }
            }
        }
        for(int i=n;i;i--)
            for(int j=1;j<i;j++)
                for(int k=1;k<=n;k++)b[j][k]=(b[j][k]-1LL*b[i][k]*a[j][i]%mod+mod)%mod;
        for(int i=1;i<=n;i++)scanf("%d",&c[i]);
        int ans=0;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)ans=(ans+1LL*c[i]*c[j]%mod*b[j][i]%mod+mod)%mod;
        printf("%d\n",ans);
    }
}