結果

問題 No.1479 Matrix Eraser
ユーザー ayoiyouyoeyoayoiyouyoeyo
提出日時 2021-04-16 21:35:48
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 559 ms / 3,000 ms
コード長 2,751 bytes
コンパイル時間 4,179 ms
コンパイル使用メモリ 192,396 KB
実行使用メモリ 30,848 KB
最終ジャッジ日時 2024-07-03 01:00:14
合計ジャッジ時間 14,930 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,812 KB
testcase_02 AC 3 ms
6,944 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 3 ms
6,940 KB
testcase_06 AC 2 ms
6,944 KB
testcase_07 AC 54 ms
6,944 KB
testcase_08 AC 94 ms
8,832 KB
testcase_09 AC 222 ms
14,336 KB
testcase_10 AC 448 ms
22,656 KB
testcase_11 AC 267 ms
16,128 KB
testcase_12 AC 68 ms
7,424 KB
testcase_13 AC 96 ms
8,832 KB
testcase_14 AC 74 ms
7,680 KB
testcase_15 AC 15 ms
6,940 KB
testcase_16 AC 82 ms
8,064 KB
testcase_17 AC 549 ms
25,984 KB
testcase_18 AC 550 ms
25,984 KB
testcase_19 AC 558 ms
26,112 KB
testcase_20 AC 546 ms
25,984 KB
testcase_21 AC 545 ms
25,984 KB
testcase_22 AC 549 ms
25,984 KB
testcase_23 AC 542 ms
25,984 KB
testcase_24 AC 553 ms
25,984 KB
testcase_25 AC 552 ms
26,112 KB
testcase_26 AC 559 ms
25,984 KB
testcase_27 AC 168 ms
9,472 KB
testcase_28 AC 168 ms
9,344 KB
testcase_29 AC 169 ms
9,216 KB
testcase_30 AC 169 ms
9,088 KB
testcase_31 AC 171 ms
9,216 KB
testcase_32 AC 102 ms
20,052 KB
testcase_33 AC 103 ms
20,112 KB
testcase_34 AC 102 ms
20,060 KB
testcase_35 AC 102 ms
20,060 KB
testcase_36 AC 102 ms
20,056 KB
testcase_37 AC 45 ms
11,328 KB
testcase_38 AC 195 ms
7,552 KB
testcase_39 AC 416 ms
30,848 KB
testcase_40 AC 2 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <string>
#include <cmath>
#include<algorithm>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<iomanip>
#include<bitset>
#define _USE_MATH_DEFINES
#include <math.h>
#include <functional>
#include<complex>
#include<cassert>
#include<random>
using namespace std;
#include<atcoder/all>
using namespace atcoder;


#define rep(i,x) for(ll i=0;i<x;i++)
#define repn(i,x) for(ll i=1;i<=x;i++)

typedef long long ll;
typedef pair<ll, ll> P;
const ll INF = 1e17;
const ll MAX = 4000001;
const long double eps = 1E-14;

ll max(ll a, ll b) {
	if (a > b) { return a; }
	return b;
}

ll min(ll a, ll b) {
	if (a > b) { return b; }
	return a;
}

ll gcd(ll a, ll b) {
	if (b == 0) { return a; }
	if (a < b) { return gcd(b, a); }
	return gcd(b, a % b);
}

ll lcm(ll a, ll b) {
	return a / gcd(a, b) * b;
}

struct edge {
	ll id;
	ll fr;
	ll to;
	ll d;
};

using mint = modint;

typedef vector<ll> vll;
typedef vector <vector<ll>> vvll;
typedef vector<vector<vector<ll>>> vvvll;

typedef vector<mint> vmint;
typedef vector<vector<mint>> vvmint;
typedef vector<vector<vector<mint>>> vvvmint;

vmint f, finv, inv;

void cominit(ll N) {
	ll MOD = modint::mod();//デフォルトが998244353に注意

	f.assign(N + 1, 1);
	finv.assign(N + 1, 1);
	inv.assign(N + 1, 1);
	inv[1] = 1;

	repn(i, N) {
		f[i] = f[i - 1] * i;
		if (i > 1)inv[i] = -inv[MOD % i] * (MOD / i);
		finv[i] = finv[i - 1] * inv[i];
	}
}

mint com(ll a, ll b) {
	if (a < 0 || b < 0 || a < b) { return 0; }
	return f[a] * finv[b] * finv[a - b];
}


/////////////////////////////////////


int main() {
	ll H, W;
	cin >> H >> W;

	map<ll, vector<P>> mp;

	rep(i, H)rep(j, W) {
		ll a;
		cin >> a;
		mp[a].push_back({ i,j });
	}

	ll ans = 0;

	for (auto p : mp) {
		if (p.first == 0) { continue; }

		vector<P> v = p.second;
		vll lisx, lisy;
		rep(i, v.size()) {
			lisx.push_back(v[i].first);
			lisy.push_back(v[i].second);
		}
		sort(lisx.begin(), lisx.end());
		sort(lisy.begin(), lisy.end());

		vll zx, zy;
		rep(i, lisx.size()) {
			if (i == 0 || lisx[i] != lisx[i - 1]) { zx.push_back(lisx[i]); }
			if (i == 0 || lisy[i] != lisy[i - 1]) { zy.push_back(lisy[i]); }
		}

		ll X = zx.size();
		ll Y = zy.size();

		mf_graph<ll> g(X + Y + 2);
		rep(i, X) { g.add_edge(X + Y, i,1); }
		rep(i, Y) { g.add_edge(X + i, X + Y + 1, 1); }

		rep(i, v.size()) {
			//cout << v[i].first << " " << v[i].second << endl;

			auto itr = lower_bound(zx.begin(), zx.end(), v[i].first);
			ll x = distance(zx.begin(), itr);

			auto itr2 = lower_bound(zy.begin(), zy.end(), v[i].second);
			ll y = distance(zy.begin(), itr2);

			g.add_edge(x, X + y, 1);
		}

		ll now = g.flow(X + Y, X + Y + 1);
		//cout << now << endl;

		ans += now;
	}

	cout << ans << endl;
}
0