結果

問題 No.1078 I love Matrix Construction
ユーザー leaf_1415
提出日時 2020-06-12 22:55:54
言語 C++11(廃止可能性あり)
(gcc 13.3.0)
結果
AC  
実行時間 250 ms / 2,000 ms
コード長 2,475 bytes
コンパイル時間 798 ms
コンパイル使用メモリ 86,064 KB
実行使用メモリ 58,436 KB
最終ジャッジ日時 2024-06-24 05:41:21
合計ジャッジ時間 5,675 ms
ジャッジサーバーID
(参考情報)
judge3 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <iostream>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <cstdlib>
#include <cassert>
#include <vector>
#include <list>
#include <stack>
#include <queue>
#include <deque>
#include <map>
#include <set>
#include <bitset>
#include <string>
#include <algorithm>
#include <utility>
#define llint long long
#define inf 1e18
#define rep(x, s, t) for(llint (x) = (s); (x) < (t); (x)++)
#define Rep(x, s, t) for(llint (x) = (s); (x) <= (t); (x)++)
#define chmin(x, y) (x) = min((x), (y))
#define chmax(x, y) (x) = max((x), (y))
using namespace std;
typedef pair<llint, llint> P;
llint n;
llint s[505], t[505], u[505];
vector<llint> G[600005], revG[600005];
vector<int> topo;
bool used[600005];
int scc[600005], pos[600005];
int ans[505][505];
void tpsort(int v)
{
used[v] = true;
for(int i = 0; i < G[v].size(); i++){
if(!used[G[v][i]]) tpsort(G[v][i]);
}
topo.push_back(v);
}
void sccdfs(int v, int id)
{
used[v] = true;
scc[v] = id;
for(int i = 0; i < revG[v].size(); i++){
if(!used[revG[v][i]]) sccdfs(revG[v][i], id);
}
}
int main(void)
{
ios::sync_with_stdio(0);
cin.tie(0);
cin >> n;
for(int i = 1; i <= n; i++) cin >> s[i];
for(int i = 1; i <= n; i++) cin >> t[i];
for(int i = 1; i <= n; i++) cin >> u[i];
llint N = (n+1)*n;
for(int i = 1; i <= n; i++){
for(int j = 1; j <= n; j++){
llint S = s[i]*n+j, T = j*n+t[i];
llint a = u[i]%2, b = u[i]/2;
G[a*N+S].push_back((1-b)*N+T);
G[b*N+T].push_back((1-a)*N+S);
}
}
for(int i = 1; i <= 2*N; i++){
for(int j = 0; j < G[i].size(); j++){
llint u = G[i][j];
revG[u].push_back(i);
}
}
for(int i = 1; i <= 2*N; i++) if(!used[i]) tpsort(i);
reverse(topo.begin(), topo.end());
int id = 1;
for(int i = 1; i <= 2*N; i++) used[i] = false;
for(int i = 0; i < topo.size(); i++) if(!used[topo[i]]) sccdfs(topo[i], id++);
for(int i = 0; i < topo.size(); i++) pos[topo[i]] = i;
for(int i = 1; i <= n; i++){
for(int j = 1; j <= n; j++){
int v = i*n+j;
//cout << scc[v] << " " << scc[N+v] << endl;
if(scc[v] == scc[N+v]){
cout << -1 << endl;
return 0;
}
if(scc[v] < scc[N+v]) ans[i][j] = 1;
else ans[i][j] = 0;
}
}
for(int i = 1; i <= n; i++){
for(int j = 1; j <= n; j++){
if(ans[s[i]][j] + ans[j][t[i]] * 2 == u[i]){
cout << -1 << endl;
return -1;
}
}
}
for(int i = 1; i <= n; i++){
for(int j = 1; j <= n; j++){
cout << ans[i][j];
if(j < n) cout << " ";
}
cout << endl;
}
return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0