結果

問題 No.1479 Matrix Eraser
ユーザー ayoiyouyoeyoayoiyouyoeyo
提出日時 2021-04-16 21:35:48
言語 C++14
(gcc 13.3.0 + boost 1.87.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

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