結果
問題 | No.1078 I love Matrix Construction |
ユーザー |
|
提出日時 | 2020-06-12 21:56:55 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 256 ms / 2,000 ms |
コード長 | 2,111 bytes |
コンパイル時間 | 970 ms |
コンパイル使用メモリ | 79,088 KB |
実行使用メモリ | 49,284 KB |
最終ジャッジ日時 | 2024-06-24 04:56:36 |
合計ジャッジ時間 | 6,763 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 22 |
コンパイルメッセージ
main.cpp:64:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type] 64 | main() | ^~~~
ソースコード
#include<iostream>using namespace std;//Strongly Connected Components#include<vector>struct SCC{int n;vector<int>comp,order;vector<bool>used;vector<vector<int> >G,RG;SCC(int _n=0):n(_n),comp(_n,-1),used(_n,false),G(_n),RG(_n){}void add_edge(int from,int to){G[from].push_back(to);RG[to].push_back(from);}void copy(const vector<vector<int> >&H){for(int i=0;i<H.size();i++){for(int j=0;j<H[i].size();j++){G[i].push_back(H[i][j]);RG[H[i][j]].push_back(i);}}}int operator[](int u)const{return comp[u];}void dfs(int u){used[u]=true;for(int i=0;i<G[u].size();i++)if(!used[G[u][i]])dfs(G[u][i]);order.push_back(u);}void rdfs(int u,int cnt){comp[u]=cnt;for(int i=0;i<RG[u].size();i++)if(comp[RG[u][i]]==-1)rdfs(RG[u][i],cnt);}int build(){for(int i=0;i<n;i++)if(!used[i])dfs(i);int cnt=0;for(int i=n-1;i>=0;i--)if(comp[order[i]]==-1)rdfs(order[i],cnt++);return cnt;}int build(vector<vector<int> >&H){int ret=build();H.assign(ret,vector<int>());for(int i=0;i<n;i++){for(int j=0;j<G[i].size();j++){if(comp[i]!=comp[G[i][j]])H[comp[i]].push_back(comp[G[i][j]]);}}return ret;}};int N;int S[500],T[500],U[500];int A[500][500];main(){cin>>N;for(int i=0;i<N;i++){cin>>S[i];S[i]--;}for(int i=0;i<N;i++){cin>>T[i];T[i]--;}for(int i=0;i<N;i++)cin>>U[i];SCC P(N*N*2);for(int i=0;i<N;i++)for(int j=0;j<N;j++){int x=S[i]*N+j,y=j*N+T[i];if(U[i]==0){P.add_edge(N*N+x,y);P.add_edge(N*N+y,x);}else if(U[i]==1){P.add_edge(x,y);P.add_edge(N*N+y,N*N+x);}else if(U[i]==2){P.add_edge(N*N+x,N*N+y);P.add_edge(y,x);}else{P.add_edge(x,N*N+y);P.add_edge(y,N*N+x);}}P.build();bool out=false;for(int i=0;i<N;i++)for(int j=0;j<N;j++){int id=i*N+j;if(P[id]==P[N*N+id]){out=true;}else{A[i][j]=P[id]>P[N*N+id]?1:0;}}if(out)cout<<-1<<endl;else{for(int i=0;i<N;i++){for(int j=0;j<N;j++)cout<<A[i][j]<<(j+1==N?"\n":" ");}}}