BZOJ3168: [Heoi2013]钙铁锌硒维生素

传送门
题意
给定一个满秩的矩阵 AAA ，另一个矩阵 BBB
对于 AAA 的每个行向量 AiA_iAi 找到一个匹配 BBB 的行向量 BpiB_{p_i}Bpi
使得 AiA_iAi 替换成 BpiB_{p_i}Bpi 后的矩阵 A′A'A′ 仍然满秩
Sol
如果 BpiB_{p_i}Bpi 能替换 AiA_iAi，那么说明存在一个向量 λ\lambdaλ 使得 λ×A′=Bpi\lambda \times A'=B_{p_i}λ×A′=Bpi
而 AAA 线性无关，只要 AAA 的元素表示 BpiB_{p_i}Bpi 的线性组合中 AiA_iAi 的系数不为 000 就好了
所以只要求出系数矩阵 CCC，使得 C×A=BC\times A=BC×A=B
然后对于不为 000 的 Ci,jC_{i,j}Ci,j，将 jjj 向 iii 连边，求这个二分图的字典序最小的完备匹配就好了
求 CCC 可以矩阵求逆实现，C=B×A−1C=B\times A^{-1}C=B×A−1
求二分图的字典序最小的完备匹配，可以先跑一遍匈牙利算法
然后在这个基础上再跑一次增广，贪心的匹配最小的，只改变匹配比当前点大的匹配边

# include <bits/stdc++.h>
using namespace std;
typedef long long ll;const int mod(10007);int a[305][605], b[305][305], d[305][305], c[305][305];
int idx, n, ans, match[305], chos[305], cnt, vis[305];inline int Pow(int x, int y) {register int ret = 1;for (; y; y >>= 1, x = x * x % mod)if (y & 1) ret = ret * x % mod;return ret;
}inline void GetInv() {register int i, j, k, inv, m = n + n;for (i = 1; i <= n; ++i) a[i][i + n] = 1;for (i = 1; i <= n; ++i) {for (j = i; j <= n; ++j)if (a[j][i]) {if (i != j) swap(a[i], a[j]);break;}inv = Pow(a[i][i], mod - 2);for (j = i; j <= m; ++j) a[i][j] = a[i][j] * inv % mod;for (j = 1; j <= n; ++j)if (i != j)for (inv = a[j][i], k = i; k <= m; ++k) a[j][k] = (a[j][k] - a[i][k] * inv % mod + mod) % mod;}for (i = 1; i <= n; ++i)for (j = 1; j <= n; ++j) d[i][j] = a[i][j + n];
}inline void GetEdge() {register int i, j, k;for (i = 1; i <= n; ++i)for (j = 1; j <= n; ++j)for (k = 1; k <= n; ++k)(c[i][k] += b[i][j] * d[j][k] % mod) %= mod;
}int Dfs(int u) {register int v;for (v = 1; v <= n; ++v)if (c[v][u] && vis[v] != idx) {vis[v] = idx;if (!match[v] || Dfs(match[v])) return chos[u] = v, match[v] = u, 1;}return 0;
}int Change(int u, int rt) {register int v;for (v = 1; v <= n; ++v)if (c[v][u] && vis[v] != idx) {vis[v] = idx;if (!match[v] || (match[v] > rt && Change(match[v], rt))) return chos[u] = v, match[v] = u, 1;}return 0;
}int main() {scanf("%d", &n);register int i, j;for (i = 1; i <= n; ++i)for (j = 1; j <= n; ++j) scanf("%d", &a[i][j]);for (i = 1; i <= n; ++i)for (j = 1; j <= n; ++j) scanf("%d", &b[i][j]);GetInv(), GetEdge();for (i = 1; i <= n; ++i) ++idx, ans += Dfs(i);if (ans != n) return puts("NIE"), 0;puts("TAK");for (i = 1; i <= n; ++i) ++idx, match[chos[i]] = 0, chos[i] = 0, Change(i, i);for (i = 1; i <= n; ++i) printf("%d\n", chos[i]);return 0;
}

BZOJ3168: [Heoi2013]钙铁锌硒维生素相关推荐

Bzoj3168 [Heoi2013]钙铁锌硒维生素
Time Limit: 20 Sec Memory Limit: 256 MB Submit: 466 Solved: 151 Description 银河队选手名单出来了!小林,作为特聘的营养师 ...
BZOJ3168. [HEOI2013]钙铁锌硒维生素（线性代数＋二分图匹配）
题目链接 https://www.lydsy.com/JudgeOnline/problem.php?id=3168 题解首先,我们需要求出对于任意的 \(i, j(1 \leq i, j \leq ...
bzoj3168 [Heoi2013]钙铁锌硒维生素（矩阵求逆+匈牙利）
如果把ai替换成bj,这n个向量还是线性无关即可以. 我们求一个C∗A=B,=>C=B∗A−1C∗A=B,=>C=B∗A−1C*A=B,=>C=B*A^{-1} 矩阵求逆可以高斯消元 ...
【BZOJ3168】[Heoi2013]钙铁锌硒维生素高斯消元求矩阵的逆+匈牙利算法
[BZOJ3168][Heoi2013]钙铁锌硒维生素 Description 银河队选手名单出来了!小林,作为特聘的营养师,将负责银河队选手参加宇宙比赛的饮食.众所周知,前往宇宙的某个星球,通常要花 ...
洛谷 P4100 [HEOI2013]钙铁锌硒维生素解题报告
P4100 [HEOI2013]钙铁锌硒维生素题目描述银河队选手名单出来了!小林,作为特聘的营养师,将负责银河队选手参加宇宙比赛的饮食. 众所周知,前往宇宙的某个星球,通常要花费好长好长的时间, ...
BZOJ 3168: [Heoi2013]钙铁锌硒维生素 [线性基 Hungary 矩阵求逆]
3168: [Heoi2013]钙铁锌硒维生素题意:给一个线性无关组A,再给一个B,要为A中每个向量在B中选一个可以代替的向量,替换后仍然线性无关.判断可行和求字典序最小的解 PoPoQQQ orz ...
BZOJ 3168 Luogu P4100 [HEOI2013]钙铁锌硒维生素 (矩阵求逆、二分图匹配)
线性代数+图论好题. 题目链接: (bzoj) https://www.lydsy.com/JudgeOnline/problem.php?id=3168 (luogu) https://www.lu ...
[HEOI2013] 钙铁锌硒维生素
Description 给定2N2N个行向量A1,-,AN,B1,-,BNA_1,\dots,A_N,B_1,\dots,B_N,需要设计一种替换方案,AiA_i被BjB_j替换(i∈[1,N],j∈ ...
【BZOJ】3168: [Heoi2013]钙铁锌硒维生素
题解 Ca Fe Zn Se 显然我们既然初始矩阵就能通过线性变换变成单位矩阵,则该矩阵一定有逆没有逆输出NIE 而且因为这些向量两两正交,则表示一个向量的时候表示方法唯一那么我们求一个逆可以求出 ...

BZOJ3168: [Heoi2013]钙铁锌硒维生素

BZOJ3168: [Heoi2013]钙铁锌硒维生素相关推荐

最新文章

热门文章