Submission #99710

#	Time	Username	Problem	Language	Result	Execution time	Memory
99710		kjp4155	Coin Collecting (JOI19_ho_t4)	C++17	8 / 100	660 ms	275456 KiB

joi2019_ho_t4

#include <bits/stdc++.h>

using namespace std;

#define VA_NUM_ARGS(...) VA_NUM_ARGS_IMPL_((0,__VA_ARGS__, 6,5,4,3,2,1))
#define VA_NUM_ARGS_IMPL_(tuple) VA_NUM_ARGS_IMPL tuple
#define VA_NUM_ARGS_IMPL(_0,_1,_2,_3,_4,_5,_6,N,...) N
#define macro_dispatcher(macro, ...) macro_dispatcher_(macro, VA_NUM_ARGS(__VA_ARGS__))
#define macro_dispatcher_(macro, nargs) macro_dispatcher__(macro, nargs)
#define macro_dispatcher__(macro, nargs) macro_dispatcher___(macro, nargs)
#define macro_dispatcher___(macro, nargs) macro ## nargs
#define Debug1(a)           cout<<#a<<"="<<(a)<<"\n"
#define Debug2(a,b)         cout<<#a<<"="<<(a)<<", "<<#b<<"="<<(b)<<"\n"
#define Debug3(a,b,c)       cout<<#a<<"="<<(a)<<", "<<#b<<"="<<(b)<<", "<<#c<<"="<<(c)<<"\n"
#define Debug4(a,b,c,d)     cout<<#a<<"="<<(a)<<", "<<#b<<"="<<(b)<<", "<<#c<<"="<<(c)<<", "<<#d<<"="<<(d)<<"\n"
#define Debug5(a,b,c,d,e)   cout<<#a<<"="<<(a)<<", "<<#b<<"="<<(b)<<", "<<#c<<"="<<(c)<<", "<<#d<<"="<<(d)<<", "<<#e<<"="<<(e)<<"\n"
#define Debug6(a,b,c,d,e,f) cout<<#a<<"="<<(a)<<", "<<#b<<"="<<(b)<<", "<<#c<<"="<<(c)<<", "<<#d<<"="<<(d)<<", "<<#e<<"="<<(e)<<", "<<#f<<"="<<(f)<<"\n"
#define Debug(...) macro_dispatcher(Debug, __VA_ARGS__)(__VA_ARGS__)
#define DA(a,s,n) cout<<#a<<"=["; printarray(a,s,n); cout<<"]\n"

#define TT1 template<class T>
#define TT1T2 template<class T1, class T2>
#define TT1T2T3 template<class T1, class T2, class T3>
template<class T, size_t N> ostream& operator << (ostream& os, const array<T, N>& v);
TT1T2 ostream& operator << (ostream& os, const pair<T1, T2>& p){ return os <<"("<<p.first<<", "<< p.second<<")"; }
TT1 ostream& operator << (ostream& os, const vector<T>& v){       bool f=1;os<<"[";for(auto& i : v) { if (!f)os << ", ";os<<i;f=0;}return os << "]"; }
template<class T, size_t N> ostream& operator << (ostream& os, const array<T, N>& v) {     bool f=1;os<<"[";for(auto& i : v) { if (!f)os << ", ";os<<i;f=0;}return os << "]"; }
TT1T2 ostream& operator << (ostream& os, const set<T1, T2>&v){    bool f=1;os<<"[";for(auto& i : v) { if (!f)os << ", ";os<<i;f=0;}return os << "]"; }
TT1T2 ostream& operator << (ostream& os, const multiset<T1,T2>&v){bool f=1;os<<"[";for(auto& i : v) { if (!f)os << ", ";os<<i;f=0;}return os << "]"; }
TT1T2T3 ostream& operator << (ostream& os, const map<T1,T2,T3>& v){ bool f = 1; os << "["; for (auto& ii : v) { if (!f)os << ", "; os << "(" << ii.first << " -> " << ii.second << ") "; f = 0; }return os << "]"; }
TT1T2 ostream& operator << (ostream& os, const multimap<T1, T2>& v){ bool f = 1; os << "["; for (auto& ii : v) { if (!f)os << ", "; os << "(" << ii.first << " -> " << ii.second << ") "; f = 0; }return os << "]"; }
TT1T2 ostream& operator << (ostream& os, priority_queue<T1, T2> v) { bool f = 1; os << "["; while (!v.empty()) { auto x = v.top(); v.pop(); if (!f) os << ", "; f = 0; os << x; } return os << "]"; }
TT1T2 void printarray(const T1& a, T2 l, T2 r){ for (T2 i = l; i<=r; i++) cout << a[i] << (i<r?" ":""); }

void cio(){
    ios::sync_with_stdio(false);
    cin.tie(NULL);
    //cout << setprecision(15);
}

#define tt() printf("%.4f sec\n", (double) clock() / CLOCKS_PER_SEC )

typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;

#define Fi first
#define Se second
#define pb(x) push_back(x)
#define sz(x) (int)x.size()
#define rep(i, n) for(int i=0;i<n;i++)
#define repp(i, n) for(int i=1;i<=n;i++)
#define all(x) x.begin(), x.end()
#define geti1(X) cin >> X
#define geti2(X,Y) cin >> X >> Y
#define geti3(X,Y,Z) cin >> X >> Y >> Z
#define geti4(X,Y,Z,W) cin >> X >> Y >> Z >> W
#define GET_MACRO(_1,_2,_3,_4,NAME,...) NAME
#define geti(...) GET_MACRO(__VA_ARGS__, geti4, geti3, geti2, geti1) (__VA_ARGS__)

#define endl '\n'

typedef tuple<int,int,int> t3;
typedef tuple<int,int,int,int> t4;
void die(){printf("-1\n"); exit(0);}
int xx[8] = {1,-1,0,0,1,1,-1,-1}, yy[8] = {0,0,1,-1,1,-1,1,-1};
const ll mod = 1e9+7;

const ll INF = 1e18;
struct MCMF{
	struct edge{
		ll to, inv, cap, flow, cost;
		ll res(){
			return cap - flow;
		}
	};

	vector<vector<edge>> graph;
	vector<int> pv, pe;
	vector<ll> dist, inq;
	
	void init(int x){
		graph.resize(x+10);
		pv.resize(x+10); pe.resize(x+10);
		dist.resize(x+10);
		inq.resize(x+10);
	}

	void make_edge(int from, int to, ll cap, ll cost){
		//printf("%d -> %d | cost = %lld\n",from,to,cost);
		edge forward = {to, (int)graph[to].size(), cap, 0, cost};
		edge backward = {from, (int)graph[from].size(), 0, 0, -cost};
		graph[from].push_back(forward);
		graph[to].push_back(backward);
	}

	ll solve(int source, int sink){
		ll ans = 0;
		ll totalflow = 0;
		while(true){
			for(auto& e : dist) e = INF;
			for(auto& e : inq) e = 0;
			queue<int> q;
			q.push(source); inq[source] = 1;
			dist[source] = 0;

			while(!q.empty()){
				int cur = q.front(); q.pop();
				inq[cur] = 0;
				for(int i=0;i<(int)graph[cur].size();i++){
					auto& e = graph[cur][i];
					if( e.res() > 0 && dist[e.to] > dist[cur] + e.cost ){
						dist[e.to] = dist[cur] + e.cost;
						pv[e.to] = cur; pe[e.to] = i;
						if( inq[e.to] == 0 ){
							q.push(e.to); inq[e.to] = 1;
						}
					}
				}
			}

			if( dist[sink] == INF ) break;
			
			ll mnflow = INF;
			for( int v = sink; v != source; v = pv[v] ){
				mnflow = min( mnflow, graph[pv[v]][pe[v]].res() );
			}

			for( int v = sink; v != source; v = pv[v] ){
				int tmp = graph[pv[v]][pe[v]].inv;
				graph[pv[v]][pe[v]].flow += mnflow;
				graph[v][tmp].flow -= mnflow;
			}
			totalflow += mnflow;
			ans += dist[sink] * mnflow;
		}
		return ans;
	}

};

int N;
ll x[3000], y[3000];
int main(){
    scanf("%d",&N);
    repp(i,2*N) scanf("%lld%lld",&x[i],&y[i]);
    int source = 0, sink = 4 * N + 1;
    MCMF mf; mf.init(sink+1);
    repp(i,2*N) mf.make_edge(source, i, 1, 0);
    repp(i,2*N) mf.make_edge(2*N+i, sink, 1, 0);
    for(int i=1;i<=2*N;i++){
        for(int j=1;j<=N;j++){
            mf.make_edge(i, 2*N+j, 1, abs(x[i]-j) + abs(y[i]-1));
        }
        for(int j=1;j<=N;j++){
            mf.make_edge(i, 3*N+j, 1, abs(x[i]-j) + abs(y[i]-2));
        }
    }
    printf("%lld\n",mf.solve(source,sink));
}

Compilation message (stderr)

joi2019_ho_t4.cpp: In function 'int main()':
joi2019_ho_t4.cpp:146:10: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
     scanf("%d",&N);
     ~~~~~^~~~~~~~~
joi2019_ho_t4.cpp:147:22: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
     repp(i,2*N) scanf("%lld%lld",&x[i],&y[i]);
                 ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~

• Subtask 1: 8 / 8
• Subtask 2: 0 / 29
• Subtask 3: 0 / 63