#990699 (C++14) No.957 植林

提出ソース

結果

問題	No.957 植林
ユーザー	Katu2ou
提出日時	2024-06-24 22:47:01
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	335 ms / 2,000 ms
コード長	5,633 bytes
コンパイル時間	4,336 ms
コンパイル使用メモリ	240,860 KB
実行使用メモリ	9,760 KB
最終ジャッジ日時	2024-06-24 22:47:15
合計ジャッジ時間	13,052 ms
ジャッジサーバーID （参考情報）	judge1 / judge3

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 6,816 KB
testcase_01	AC	2 ms 6,944 KB
testcase_02	AC	2 ms 6,944 KB
testcase_03	AC	21 ms 8,956 KB
testcase_04	AC	20 ms 8,540 KB
testcase_05	AC	21 ms 8,844 KB
testcase_06	AC	23 ms 9,236 KB
testcase_07	AC	21 ms 8,812 KB
testcase_08	AC	17 ms 8,924 KB
testcase_09	AC	17 ms 8,924 KB
testcase_10	AC	17 ms 9,120 KB
testcase_11	AC	17 ms 9,044 KB
testcase_12	AC	17 ms 8,904 KB
testcase_13	AC	14 ms 7,524 KB
testcase_14	AC	17 ms 9,360 KB
testcase_15	AC	15 ms 8,944 KB
testcase_16	AC	14 ms 7,872 KB
testcase_17	AC	15 ms 8,548 KB
testcase_18	AC	248 ms 8,976 KB
testcase_19	AC	258 ms 9,072 KB
testcase_20	AC	261 ms 9,356 KB
testcase_21	AC	285 ms 9,088 KB
testcase_22	AC	300 ms 9,352 KB
testcase_23	AC	297 ms 9,480 KB
testcase_24	AC	317 ms 9,640 KB
testcase_25	AC	335 ms 9,760 KB
testcase_26	AC	329 ms 9,632 KB
testcase_27	AC	333 ms 9,640 KB
testcase_28	AC	329 ms 9,720 KB
testcase_29	AC	329 ms 9,636 KB
testcase_30	AC	329 ms 9,664 KB
testcase_31	AC	246 ms 8,952 KB
testcase_32	AC	258 ms 9,072 KB
testcase_33	AC	259 ms 9,228 KB
testcase_34	AC	281 ms 9,212 KB
testcase_35	AC	298 ms 9,356 KB
testcase_36	AC	297 ms 9,484 KB
testcase_37	AC	312 ms 9,512 KB
testcase_38	AC	328 ms 9,632 KB
testcase_39	AC	329 ms 9,760 KB
testcase_40	AC	328 ms 9,640 KB
testcase_41	AC	14 ms 9,632 KB
testcase_42	AC	14 ms 9,636 KB
testcase_43	AC	19 ms 9,636 KB
testcase_44	AC	21 ms 9,704 KB
testcase_45	AC	2 ms 6,944 KB
testcase_46	AC	2 ms 6,940 KB
testcase_47	AC	2 ms 6,944 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;

#define rep2(i, m, n) for (int i = (m); i < (n); ++i)
#define rep(i, n) rep2(i, 0, n)
#define drep2(i, m, n) for (int i = (m)-1; i >= (n); --i)
#define drep(i, n) drep2(i, n, 0)
#define all(...) std::begin(__VA_ARGS__), std::end(__VA_ARGS__)
#define rall(...) std::rbegin(__VA_ARGS__), std::rend(__VA_ARGS__)
#define FOR(i, a, b) for (int i = (a), i##_len = (b); i <= i##_len; ++i)
#define REV(i, a, b) for (int i = (a); i >= (b); --i)
#define CLR(a, b) memset((a), (b), sizeof(a))
#define DUMP(x) cout << #x << " = " << (x) << endl;
#define INF 1001001001001001001ll
#define inf (int)1001001000
#define MOD 998244353
#define MOD1 1000000007
#define PI 3.14159265358979
#define Dval 1e-12
#define fcout cout << fixed << setprecision(12)
#define Mp make_pair
#define pb push_back
#define fi first
#define se second
#define SORT(x) sort(x.begin(),x.end())
#define ERASE(x) x.erase(unique(x.begin(),x.end()),x.end())
#define POSL(x,v) (lower_bound(x.begin(),x.end(),v)-x.begin())
#define POSU(x,v) (upper_bound(x.begin(),x.end(),v)-x.begin())

using ll = long long;
using ld = long double;
using vi = vector<int>;
using vl = vector<long long>;
using vs = vector<string>;
using vd = vector<double>;
using vld = vector<long double>;
using vc = vector<char>;
using vb = vector<bool>;
using vpii = vector<pair<int, int>>;
using vpil = vector<pair<int, long long>>;
using vpll = vector<pair<long long, long long>>;
using vvi = vector<vector<int>>;
using vvl = vector<vector<long long>>;
using vvd = vector<vector<double>>;
using vvld = vector<vector<long double>>;
using vvc = vector<vector<char>>;
using vvb = vector<vector<bool>>;
using vvpii = vector<vector<pair<int,int>>>;
using vvpll = vector<vector<pair<long long,long long>>>;
using vvvi = vector<vector<vector<int>>>;
using vvvl = vector<vector<vector<long long>>>;
using pii = pair<int, int>;
using pll = pair<long long, long long>;
using LL = __int128_t;


template<typename T>
T dgt(T b, T n){ T cnt=0; while(n){ cnt++; n/=b;} return cnt;}
ll gcd(ll x, ll y) {	if (x == 0) return y;	return gcd(y%x, x);} ll lcm(ll x, ll y) { __int128_t xx,yy; xx=x; yy=y; __int128_t ans=xx * yy / gcd(x, y); ll ans2=ans; return ans; }
template<typename T>
T POW(T x, ll n){T ret=1;	while(n>0){		if(n&1) ret=ret*x;		x=x*x;		n>>=1;	}	return ret;}
template<typename T>
T modpow(T a, ll n, T p) {	if(n==0) return (T)1;  if (n == 1) return a % p;  if (n % 2 == 1) return (a * modpow(a, n - 1, p)) % p;  T t = modpow(a, n / 2, p);  return (t * t) % p;}
template<typename T>
T modinv(T a, T m) {	if(m==0)return (T)1;	T b = m, u = 1, v = 0;	while (b) {		T t = a / b;		a -= t * b; swap(a, b);		u -= t * v; swap(u, v);	}	u %= m;	if (u < 0) u += m;	return u;}
template<typename T>
T REM(T a, T b){ return (a % b + b) % b;}
template<typename T>
T QUO(T a, T b){ return (a - REM(a, b)) / b;}
const int MAXCOMB=510000;
ll MODCOMB = 998244353;
ll fac[MAXCOMB], finv[MAXCOMB], inv[MAXCOMB]; 
void COMinit() {	fac[0] = fac[1] = 1;	finv[0] = finv[1] = 1;	inv[1] = 1;	for (int i = 2; i < MAXCOMB; i++) {		fac[i] = fac[i - 1] * i % MODCOMB;		inv[i] = MODCOMB - inv[MODCOMB%i] * (MODCOMB / i) % MODCOMB;		finv[i] = finv[i - 1] * inv[i] % MODCOMB;	}}
ll COM(ll n, ll k) {	if (n < k) return 0;	if (n < 0 || k < 0) return 0;	return fac[n] * (finv[k] * finv[n - k] % MODCOMB) % MODCOMB;}
ll com(ll n,ll m){   if(n<m || n<=0 ||m<0){		return 0;	}	if( m==0 || n==m){		return 1;	}	ll k=1;	for(ll i=1;i<=m;i++){       k*=(n-i+1); 	   k%=MODCOMB;	   k*=modinv(i,MODCOMB);	   k%=MODCOMB;	}	return k;}
ll rad(ll u, ll p){ ll cnt=0;	while(u%p==0){		u/=p;		cnt++;	}	return cnt;}

template <typename T> inline bool chmax(T &a, T b) { return ((a < b) ? (a = b, true) : (false));}
template <typename T> inline bool chmin(T &a, T b) { return ((a > b) ? (a = b, true) : (false));}
template <class T> T BS(vector<T> &vec, T key) { auto itr = lower_bound(vec.begin(), vec.end(), key); return distance(vec.begin(), itr); }
template<class T> pair<T,T> RangeBS(vector<T> &vec, T lowv, T highv){auto itr_l = lower_bound(vec.begin(), vec.end(), lowv); auto itr_r = upper_bound(vec.begin(), vec.end(), highv); return make_pair(distance(vec.begin(), itr_l), distance(vec.begin(), itr_r)-1);}
void fail() { cout << "-1\n"; exit(0); } void no() { cout << "No\n"; exit(0); } void yes() { cout << "Yes\n"; exit(0); }
template<class T> void er(T a) { cout << a << '\n'; exit(0); }
int dx[] = { 1,0,-1,0,1,1,-1,-1,0 }; int dy[] = { 0,1,0,-1,1,-1,1,-1,0 };
bool range_in(int i, int j, int h, int w){ if(i<0 || j<0 || i>=h || j>=w) return false; return true;} 
int bitcount(int n){n=(n&0x55555555)+(n>>1&0x55555555); n=(n&0x33333333)+(n>>2&0x33333333); n=(n&0x0f0f0f0f)+(n>>4&0x0f0f0f0f); n=(n&0x00ff00ff)+(n>>8&0x00ff00ff); n=(n&0x0000ffff)+(n>>16&0x0000ffff); return n;}

void solve(int &FIN){
	int h,w;
	cin>>h>>w;
	vvl cost(h,vl(w));
	rep(i,h){
		rep(j,w){
			cin>>cost[i][j];
		}
	}
	vl gh(h);
	vl gw(w);
	ll ans=0;
	rep(i,h){
		cin>>gh[i];
		ans+=gh[i];
	}
	rep(i,w){
		cin>>gw[i];
		ans+=gw[i];
	}
	mf_graph<ll> G(h+w+2);
	
	rep(i,h){
        ll val=0;
		rep(j,w){
            val+=cost[i][j];
        }
        G.add_edge(h+w,i,gh[i]);
        G.add_edge(i,h+w+1,val);
	}
    rep(j,w){
        G.add_edge(h+w,h+j,gw[j]);
    }
    rep(i,h){
        rep(j,w){
            G.add_edge(h+j,i,cost[i][j]);
        }
    }
    ll f=G.flow(h+w,h+w+1);
	
	cout<<ans-f<<endl;

}

signed main(){
	cin.tie(0);
	ios::sync_with_stdio(0);
	cout<<fixed<<setprecision(20);
    int FIN=0;
	while(1){
        solve(FIN);
        FIN=1;
        if(FIN)return 0;
    }
}

yukicoder

結果

テストケース

ソースコード