zoukankan      html  css  js  c++  java
  • UR#6 懒癌

    超级结论题。。。

    首先建出原图的补图,再考虑缩点。

    缩完点后是一个有向无环图。
    这个图的每个点代表一个强联通分量。
    考虑一个大小大于(1)的强联通分量中的一个点(i),那么(i)永远不可能开枪。

    因为(i)开枪需要假设自己的狗没病,然后枚举自己看不到的狗的状态,然而由于有他的转移依赖的状态 依赖他,那么就会循环转移。
    导致他一直到inf都不会开枪。

    那么任何可以到达大于1的强联通分量的强联通分量都废了(包括本身),因为他们中所有点都只能不是病狗。

    现在每个强联通分量都是单点强联通分量。

    考虑一个状态,每个病狗的主人都非常聪明但也非常谨慎,它从他的狗没病,并且他不知道的狗状态随便的状态 中 开枪时间最大 的一种 转移过来。
    然后所有主人的答案取(min)
    这样一次最多在 可到达点集 中去掉一个点,且肯定可以去掉一个点,因为拓扑序最小的点肯定去掉后再也不能被染黑了。
    那么可以容易得出,状态的贡献就是它可以到的点集大小。再计算点的贡献。

    考虑只有凸集上的点才能贡献答案(因为你取min的时候只能选凸集上的点),第二问也可以用点的贡献做。

    #include <bits/stdc++.h>
    
    
    using namespace std;
    const int N=3010;
    typedef long long ll;
    int dfn[N],low[N],clk;
    int isc[N],sz[N],nscc;
    stack<int> st;
    vector<int> e[N],g[N];
    void Tajan(int x){
    	//cerr<<"Tajan"<<x<<endl;
    	dfn[x]=low[x]=++clk;
    	st.push(x);
    	for (auto i:e[x]){
    		if (!dfn[i]) Tajan(i);
    		if (!isc[i]) low[x]=min(low[x],low[i]);
    	}
    	if (dfn[x]==low[x]){
    		++nscc;
    		while (1){
    			int u=st.top(); st.pop();
    			isc[u]=nscc;
    			++sz[nscc];
    			//cerr<<"pop"<<nscc<<" "<<u<<endl;
    			if (u==x) break;
    		}
    	}
    }
    bitset<N> b[N];
    
    const int M=998244353;
    
    int add(int x,int y){
    	return (x+=y)>=M?x-M:x;
    }
    int sub(int x,int y){
    	return (x-=y)<0?x+M:x;
    }
    int fp(int x,int y){
    	int ret=1;
    	for (; y; y>>=1,x=(ll)x*x%M) if (y&1) ret=(ll)ret*x%M;
    	return ret;
    }
    char s[N];
    int n;
    int fc[N];
    vector<int> ups[N];
    int f[N];
    int main(){
    	scanf("%d",&n);
    	for (int i=1; i<=n; ++i){
    		scanf("%s",s+1);
    		for (int j=1; j<=n; ++j)
    			if (j!=i&&s[j]=='0'){
    				//cerr<<i<<" "<<j<<endl;
    				e[i].push_back(j);
    			}
    	}
    
    	for (int i=1; i<=n; ++i){
    		//cerr<<"????"<<i<<endl;
    		if (!dfn[i]) Tajan(i);
    	}
    	//cerr<<nscc<<endl;
    	for (int i=1; i<=n; ++i){
    		for (auto j:e[i]){
    			int x=isc[i];
    			int y=isc[j];
    			ups[y].push_back(x);
    			if (sz[x]>1||sz[y]>1) continue;
    			g[x].push_back(y);
    		}
    	}
    	queue<int> q;
    	for (int i=1; i<=nscc; ++i)
    		if (sz[i]>1) q.push(i),f[i]=0; else f[i]=1000000;
    	while (!q.empty()){
    		int x=q.front(); q.pop();
    		for (auto j:ups[x])
    			if (f[j]==1000000){
    				f[j]=f[x]+1;
    				q.push(j);
    			}
    	}
    	int tot=0,ans1=0,ans2=0;
    	for (int i=1; i<=nscc; ++i) tot+=(sz[i]==1&&f[i]==1000000);
    	for (int i=1; i<=n; ++i) fc[isc[i]]=i;
    	for (int i=nscc; i>=1; --i)
    		if (sz[i]==1&&f[i]==1000000){
    			//cerr<<fc[i]<<endl;
    			//cerr<<i<<" "<<fc[i]<<endl;
    			b[i].set(fc[i]);
    			for (auto j:g[i]) b[j]|=b[i];
    			ans1=add(ans1,(ll)sub(fp(2,b[i].count()),1)*fp(2,tot-b[i].count())%M);
    			ans2=add(ans2,fp(2,tot-b[i].count()));
    		}
    	cout<<ans1<<" "<<ans2<<endl;
    }
    
  • 相关阅读:
    第一课:人人站安装与使用教程
    linux 多个python版本的切换
    参考文章
    windows下scrapy 的安装
    vi 使用入门
    linux 源码安装
    C语言—第二次作业
    C语言第0次作业
    c语言博客园作业03函数
    tomcat对于web.xml的securityconstraint使用的处理机制
  • 原文地址:https://www.cnblogs.com/Yuhuger/p/10560512.html
Copyright © 2011-2022 走看看