期末考炸飞了啊QAQ求安慰。
首先特判不合法的情况;然后考虑没有权值为0的边的时候,一条边(i,j)如果可有可无(也就是存在k使得d(i,k)+d(k,j)=d(i,j)),说明i->j可以不经过(i,j),也就是这条边只要取>=当前权值且<=k的所有值即可;否则这条边不能变化。
有0边的时候,把用0边连的全部缩成一个团然后像上述类似地处理不同团的连边;0边连的就是一个dp就好了。
AC代码如下:
#include<bits/stdc++.h>
#define ll long long
#define mod 998244353
#define N 405
using namespace std;
int n,m,a[N][N],cbn[N][N],fa[N],sz[N],f[N],g[N];
bool ok(){
int i,j,k;
for (i=1; i<=n; i++) if (a[i][i]) return 0;
for (i=1; i<=n; i++)
for (j=1; j<=n; j++){
if (a[i][j]>m) return 0; if (a[i][j]!=a[j][i]) return 0;
for (k=1; k<=n; k++) if (a[i][j]+a[j][k]<a[i][k]) return 0;
}
return 1;
}
int getfa(int x){ return (x==fa[x])?x:fa[x]=getfa(fa[x]); }
int pw(int x,int y){
int z=1; for (; y; y>>=1,x=(ll)x*x%mod) if (y&1) z=(ll)z*x%mod; return z;
}
int main(){
scanf("%d%d",&n,&m);
int i,j,k;
for (i=1; i<=n; i++)
for (j=1; j<=n; j++) scanf("%d",&a[i][j]);
if (!ok()){ puts("0"); return 0; }
for (i=1; i<=n; i++) fa[i]=i;
for (i=1; i<n; i++)
for (j=i+1; j<=n; j++) if (!a[i][j]) fa[getfa(i)]=getfa(j);
for (i=1; i<=n; i++) sz[getfa(i)]++;
int ans=1;
for (i=1; i<n; i++) if (fa[i]==i)
for (j=i+1; j<=n; j++) if (fa[j]==j && i!=j){
for (k=1; k<=n; k++) if (fa[k]==k && i!=k && k!=j && a[i][k]+a[k][j]==a[i][j]) break;
if (k>n) ans=(ll)ans*(pw(m-a[i][j]+1,sz[i]*sz[j])-pw(m-a[i][j],sz[i]*sz[j])+mod)%mod;
else ans=(ll)ans*pw(m-a[i][j]+1,sz[i]*sz[j])%mod;
}
for (i=0; i<=n; i++)
for (j=cbn[i][0]=1; j<=i; j++) cbn[i][j]=(cbn[i-1][j-1]+cbn[i-1][j])%mod;
for (i=1; i<=n; i++){
f[i]=g[i]=pw(m+1,i*(i-1)>>1);
for (j=1; j<i; j++) f[i]=(f[i]-(ll)f[j]*g[i-j]%mod*cbn[i-1][j-1]%mod*pw(m,j*(i-j))%mod+mod)%mod;
}
for (i=1; i<=n; i++) if (fa[i]==i) ans=(ll)ans*f[sz[i]]%mod;
printf("%d\n",ans);
return 0;
}
by lych
2017.1.15