Submission #1202318

#TimeUsernameProblemLanguageResultExecution timeMemory
1202318t9unkubjVirus Experiment (JOI19_virus)C++20
20 / 100
2096 ms124648 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 struct scc{ vector<vector<int>>g; vector<vector<int>>inv; int n; scc(int n):n(n),g(n),inv(n){} void add(int x,int y){ g[x].emplace_back(y); inv[y].emplace_back(x); } vector<vector<int>> do_scc(){ vector<int>seen(n),seen2(n); vector<vector<int>>ret; vector<int>num; for(int i=0;i<n;i++){ if(!seen[i]){ auto dfs=[&](auto&dfs,int now,int pre)->void{ seen[now]=1; for(auto x:g[now]){ if(x==pre)continue; if(!seen[x])dfs(dfs,x,now); } num.emplace_back(now); }; dfs(dfs,i,-1); } } for(int i=n-1;i>=0;i--){ int x=num[i]; if(!seen2[x]){ vector<int>res; seen2[x]=1; auto dfs=[&](auto&dfs,int now,int pre)->void{ seen2[now]=1; for(auto x:inv[now]){ if(x==pre)continue; if(!seen2[x])dfs(dfs,x,now); } res.emplace_back(now); }; dfs(dfs,x,-1); ret.push_back(res); } } return ret; } }; 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>; 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){ scc scc(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))scc.add(x,dsu.leader(i*c+j)); } } int update=0; for(auto&x:scc.do_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(scc.g[j].size()==0){ chmin(mi,dsu.size(j)); } } } int ans=0; rep(j,r*c){ if(dsu.leader(j)==j&&u[j/c][j%c]){ if(scc.g[j].size()==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...