标签：Platinum Subsets int rep USACO ans y1 x1 dp

「USACO 2021 US Open Platinum」Balanced Subsets

考虑题目给出的定义对应怎样的图形，显然是一个凸的封闭图形

不妨通过左右边线描述，

1.左边线先左移再右移

2.右边线先右移再左移

不妨直接令\(dp_{i,l,r,f1,f2}\)表示当前第\(i\)行，当前左右边线为\(l,r(l\leq r,\forall j\in[l,r],a_{i,j}=G)\)

\(f1,f2\)表示当前左右边线处于左移还是右移状态

以左边线为例，定义右移开始的时刻为第一个\(l>l'\)的时刻

容易得到转移，是一个前/后缀和的形式

那么对于\([l,r]\)两维分别做前缀和，然后\(O(n^3)\)转移即可

注意转移过程中要确保\([l',r'],[l,r]\)有交

以下是暴力二维前缀和+手艹9种转移的代码

const int N=170,P=1e9+7;

int n;
int c[N];
char s[N];
int dp[N][N][2][2];
int F[N][N][2][2];
int S(int a,int b,int x1,int x2,int y1,int y2){
	x1--,y1--;
	return (0ll+F[x2][y2][a][b]-F[x1][y2][a][b]-F[x2][y1][a][b]+F[x1][y1][a][b])%P;
}

int main(){
	n=rd();
	int ans=0;
	rep(i,1,n) {
		scanf("%s",s+1);
		rep(j,1,n) c[j]=c[j-1]+(s[j]=='G');
		memset(F,0,sizeof F);
		rep(i,1,n) rep(j,1,n) rep(a,0,1) rep(b,0,1) {
			F[i][j][a][b]=(0ll+F[i-1][j][a][b]+F[i][j-1][a][b]-F[i-1][j-1][a][b]+dp[i][j][a][b])%P;
		}
		memset(dp,0,sizeof dp);
		rep(l,1,n) rep(r,l,n) if(c[r]-c[l-1]==r-l+1) {
			(dp[l][r][0][0]+=S(0,0,l,r,l,r))%=P;
			(dp[l][r][1][0]+=S(0,0,1,l-1,l,r))%=P;
			(dp[l][r][0][1]+=S(0,0,l,r,r+1,n))%=P;
			(dp[l][r][1][1]+=S(0,0,1,l-1,r+1,n))%=P;
	
			(dp[l][r][0][1]+=S(0,1,l,r,r,n))%=P;
			(dp[l][r][1][1]+=S(0,1,1,l-1,r,n))%=P;
			
			(dp[l][r][1][0]+=S(1,0,1,l,l,r))%=P;
			(dp[l][r][1][1]+=S(1,0,1,l,r+1,n))%=P;
			
			(dp[l][r][1][1]+=S(1,1,1,l,r,n))%=P;
		}
		rep(l,1,n) rep(r,l,n) if(c[r]-c[l-1]==r-l+1) dp[l][r][0][0]++;
		rep(l,1,n) rep(r,l,n) rep(a,0,1) rep(b,0,1) if(dp[l][r][a][b]) 
			ans=(ans+dp[l][r][a][b])%P;
	}
	Mod2(ans),printf("%d\n",ans);
}

标签：Platinum,Subsets,int,rep,USACO,ans,y1,x1,dp
来源： https://www.cnblogs.com/chasedeath/p/14764103.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。