Submission #1202327

#TimeUsernameProblemLanguageResultExecution timeMemory
1202327t9unkubjVirus Experiment (JOI19_virus)C++20
100 / 100
1718 ms147788 KiB
#ifdef t9unkubj #include"my_template.h" //#include"my_template_no_debug.h" #else #define dbg(...) 199958 #pragma GCC optimize("O3") #include"bits/stdc++.h" using namespace std; using uint=unsigned; using ll=long long; using ull=unsigned long long; using ld=long double; using i128=__int128; template<class T>using vc=vector<T>; template<class T>using vvc=vc<vc<T>>; template<class T>using vvvc=vc<vvc<T>>; template<class T>using vvvvc=vc<vvvc<T>>; template<class T>using smpq=priority_queue<T,vc<T>,greater<T>>; template<class T>using bipq=priority_queue<T>; using vi=vc<int>; using vvi=vvc<int>; using vl=vc<ll>; using vvl=vvc<ll>; #define rep(i,n) for(ll i=0;i<(ll)(n);i++) #define REP(i,j,n) for(ll i=(j);i<(ll)(n);i++) #define DREP(i,n,m) for(ll i=(n);i>=(m);i--) #define drep(i,n) for(ll i=(ll(n)-1);i>=0;i--) #define all(x) x.begin(),x.end() #define rall(x) x.rbegin(),x.rend() #define pb push_back template<class T,class F> bool chmin(T &x, F y){ if(x>y){ x=y; return true; } return false; } template<class T, class F> bool chmax(T &x, F y){ if(x<y){ x=y; return true; } return false; } template<class T> T sum(const vector<T>&v){ return accumulate(all(v),T(0)); } template<class T> T min(const vector<T>&v){ return *min_element(all(v)); } template<class T> T max(const vector<T>&v){ return *max_element(all(v)); } void YesNo(bool y){ cout<<(y?"Yes":"No")<<endl; } void Yes(){ cout<<"Yes"<<endl; } void No(){ cout<<"No"<<endl; } template<class T> void unique(vc<T>&a){ a.erase(unique(all(a)),a.end()); } vvi readgraph(int n,int m,int off = -1){ vvi g(n); rep(i, m){ int u,v; cin>>u>>v; u+=off,v+=off; g[u].push_back(v); g[v].push_back(u); } return g; } vvi readtree(int n,int off=-1){ return readgraph(n,n-1,off); } template<class T> vc<T> presum(vc<T> &a){ vc<T> ret(a.size()+1); rep(i,a.size())ret[i+1]=ret[i]+a[i]; return ret; } template<class T, class F> vc<T> &operator+=(vc<T> &a,F b){ for (auto&v:a)v += b; return a; } template<class T, class F> vc<T> &operator-=(vc<T>&a,F b){ for (auto&v:a)v-=b; return a; } ostream&operator<<(ostream&os,i128 num) { if(num==0){ os<<"0"; return os; } if(num<0){ os<<"-"; num=-num; } string res; while(num){ res+='0'+static_cast<int>(num%10); num/=10; } reverse(all(res)); os<<res; return os; } istream&operator>>(istream& is,i128&num) { string s; is>>s; bool neg=0; num=0; if(s[0]=='-'){ neg=1; s=s.substr(1); } for(auto&c:s){ if (c>='0'&&c<='9'){ num=num*10+(c-'0'); }else{ is.setstate(ios::failbit); return is; } } if(neg){ num=-num; } return is; } void scan(int&a) { cin >> a; } void scan(ll&a) { cin >> a; } void scan(string&a) { cin >> a; } void scan(char&a) { cin >> a; } void scan(uint&a) { cin >> a; } void scan(ull&a) { cin >> a; } void scan(bool&a) { cin >> a; } void scan(ld&a){ cin>> a; } void scan(i128&a){ cin>> a; } template<class T> void scan(vector<T>&a) { for(auto&x:a) scan(x); } void read() {} template<class Head, class... Tail> void read(Head&head, Tail&... tail) { scan(head); read(tail...); } #define INT(...) int __VA_ARGS__; read(__VA_ARGS__); #define LL(...) ll __VA_ARGS__; read(__VA_ARGS__); #define ULL(...) ull __VA_ARGS__; read(__VA_ARGS__); #define STR(...) string __VA_ARGS__; read(__VA_ARGS__); #define CHR(...) char __VA_ARGS__; read(__VA_ARGS__); #define DBL(...) double __VA_ARGS__; read(__VA_ARGS__); #define LD(...) ld __VA_ARGS__; read(__VA_ARGS__); #define I128(...) i128 __VA_ARGS__; read(__VA_ARGS__); #define VC(type, name, ...) vector<type> name(__VA_ARGS__); read(name); #define VVC(type, name, size, ...) vector<vector<type>> name(size, vector<type>(__VA_ARGS__)); read(name); void print(int a) { cout << a; } void print(ll a) { cout << a; } void print(string a) { cout << a; } void print(char a) { cout << a; } void print(uint a) { cout << a; } void print(bool a) { cout << a; } void print(ull a) { cout << a; } void print(double a) { cout << a; } void print(ld a){ cout<< a; } void print(i128 a){ cout<< a; } template<class T> void print(vector<T>a) { for(int i=0;i<(int)a.size();i++){if(i)cout<<" ";print(a[i]);}cout<<endl;} void PRT() { cout <<endl; return ; } template<class T> void PRT(T a) { print(a); cout <<endl; return; } template<class Head, class... Tail> void PRT(Head head, Tail ... tail) { print(head); cout << " "; PRT(tail...); return; } struct ioset{ ioset(){ cin.tie(0)->sync_with_stdio(0); #ifdef t9unkubj cout<<fixed<<setprecision(6); #else cout<<fixed<<setprecision(20); #endif } }ioset_______; struct MY_TIMER{ double start_time; MY_TIMER(){ start_time=clock(); } double out(){ return (clock()-start_time)/double(CLOCKS_PER_SEC); } }TIMER; #endif template<class T,auto op,int extra> struct base_dsu{ vector<int>par; vector<T>data; base_dsu(int n):par(n,-1){ static_assert(!extra,"e is needed"); } base_dsu(int n,T e):par(n,-1),data(n,e){} T& operator[](int i){ static_assert(extra,"No data"); return data[leader(i)]; } int root(int x){ if(par[x]<0)return x; else return par[x]=root(par[x]); } bool same(int x,int y){ return root(x)==root(y); } bool merge(int x,int y){ x=root(x); y=root(y); if(x==y)return false; if(par[x]>par[y])swap(x,y); par[x]+=par[y]; par[y]=x; if(extra){ data[x]=op(data[y],data[x]); } return true; } int size(int x){ return -par[root(x)]; } int leader(int x){ return root(x); } }; int tf323(int a,int b){return a;} using dsu=base_dsu<int,tf323,0>; template<class T,auto op> using extra_dsu=base_dsu<T,op,1>; #ifndef ATCODER_SCC_HPP #define ATCODER_SCC_HPP 1 #include <algorithm> #ifndef ATCODER_INTERNAL_SCC_HPP #define ATCODER_INTERNAL_SCC_HPP 1 #include <algorithm> #include <utility> #include <vector> namespace atcoder { namespace internal { template <class E> struct csr { std::vector<int> start; std::vector<E> elist; csr(){} csr(int n, const std::vector<std::pair<int, E>>& edges) : start(n + 1), elist(edges.size()) { for (auto e : edges) { start[e.first + 1]++; } for (int i = 1; i <= n; i++) { start[i] += start[i - 1]; } auto counter = start; for (auto e : edges) { elist[counter[e.first]++] = e.second; } } }; // Reference: // R. Tarjan, // Depth-First Search and Linear Graph Algorithms struct scc_graph { public: scc_graph(int n) : _n(n) {} int num_vertices() { return _n; } struct edge { int to; }; void add_edge(int from, int to) { edges.push_back({from, {to}}); } csr<edge>g; // @return pair of (# of scc, scc id) std::pair<int, std::vector<int>> scc_ids() { g = csr<edge>(_n, edges); int now_ord = 0, group_num = 0; std::vector<int> visited, low(_n), ord(_n, -1), ids(_n); visited.reserve(_n); auto dfs = [&](auto self, int v) -> void { low[v] = ord[v] = now_ord++; visited.push_back(v); for (int i = g.start[v]; i < g.start[v + 1]; i++) { auto to = g.elist[i].to; if (ord[to] == -1) { self(self, to); low[v] = std::min(low[v], low[to]); } else { low[v] = std::min(low[v], ord[to]); } } if (low[v] == ord[v]) { while (true) { int u = visited.back(); visited.pop_back(); ord[u] = _n; ids[u] = group_num; if (u == v) break; } group_num++; } }; for (int i = 0; i < _n; i++) { if (ord[i] == -1) dfs(dfs, i); } for (auto& x : ids) { x = group_num - 1 - x; } return {group_num, ids}; } std::vector<std::vector<int>> scc() { auto ids = scc_ids(); int group_num = ids.first; std::vector<int> counts(group_num); for (auto x : ids.second) counts[x]++; std::vector<std::vector<int>> groups(ids.first); for (int i = 0; i < group_num; i++) { groups[i].reserve(counts[i]); } for (int i = 0; i < _n; i++) { groups[ids.second[i]].push_back(i); } return groups; } int _n; std::vector<std::pair<int, edge>> edges; }; } // namespace internal } // namespace atcoder #endif // ATCODER_INTERNAL_SCC_HPP #include <cassert> #include <vector> namespace atcoder { struct scc_graph { public: scc_graph() : internal(0) {} scc_graph(int n) : internal(n) {} void add_edge(int from, int to) { int n = internal.num_vertices(); assert(0 <= from && from < n); assert(0 <= to && to < n); internal.add_edge(from, to); } std::vector<std::vector<int>> scc() { return internal.scc(); } internal::scc_graph internal; }; } // namespace atcoder #endif // ATCODER_SCC_HPP void solve(){ INT(m,r,c); STR(s); VVC(int,u,r,c); vc<int>ren(16); map<char,int>inv; // inv['E']=0; inv['S']=1; inv['W']=2; inv['N']=3; rep(i,16){ vc<int>v; while(v.size()<2e5){ for(auto&x:s){ v.push_back(i>>inv[x]&1); } } vc<int>val(v.size()); drep(i,v.size()){ if(i+1<v.size()&&v[i]){ chmax(val[i],val[i+1]+v[i]); } chmax(val[i],v[i]); } ren[i]=max(val); } dbg(ren); vc<int>dx{0,1,0,-1},dy{1,0,-1,0}; dsu dsu(r*c); while(1){ atcoder::scc_graph scc(r*c); vc<int>cnt(r*c); rep(i,r)rep(j,c){ unordered_map<int,int>ma; rep(k,4){ uint ni=i+dx[k]; uint nj=j+dy[k]; if(ni<r&&nj<c){ ma[dsu.leader(ni*c+nj)]|=1<<k; } } for(auto&[x,y]:ma){ if(u[i][j]&&ren[y]>=u[i][j]&&x!=dsu.leader(i*c+j))cnt[x]++,scc.add_edge(x,dsu.leader(i*c+j)); } } int update=0; for(auto&x:scc.scc()){ update|=x.size()>1; rep(i,x.size()-1)dsu.merge(x[i],x[i+1]); } if(update==0){ int mi=1e9; rep(j,r*c){ if(dsu.leader(j)==j&&u[j/c][j%c]){ if(cnt[j]==0){ chmin(mi,dsu.size(j)); } } } int ans=0; rep(j,r*c){ if(dsu.leader(j)==j&&u[j/c][j%c]){ if(cnt[j]==0&&dsu.size(j)==mi){ ans+=mi; } } } cout<<mi<<" "<<ans<<"\n"; return; } } } signed main(){ int t=1; //cin>>t; while(t--)solve(); } /* ある状況がvalidかは炉リハで事前計算 一人だけで感染させられる人々を集めて連結成分を作る */
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...