首页 > 其他分享 >CF506E Mr. Kitayuta's Gift 思考--zhengjun

CF506E Mr. Kitayuta's Gift 思考--zhengjun

时间:2023-07-24 19:45:57浏览次数:37  
标签:matrix Gift -- zhengjun int 端点 ans 转移

妙妙题。

首先可以有一个 \(O(kn^2)\) 的 dp,但是显然不行。

但是,发现其中的大多数转移都浪费在自环上了,所以考虑不要这个东西。

这个 dp 一共有三种转移:

  1. 左右端点一起向内移动一格;

  2. 左端点或右端点单独移动;

  3. 左右端点都不动。

所以考虑加一维 \(k\) 表示走了 \(k\) 次转移 1,这样可以计算出转移 2 的次数。

最后对于转移 1 走了 \(i\) 次,转移 2 走了 \(j\) 次的方案数一起计算贡献:

\[g_{i,j}=\sum\limits_{x=0}^{m-i-j}26^x\sum\limits_{y=0}^{m-i-j-x}24^y\times 25^{m-i-j-x}\times \binom{k+j-1}{j-1}\times \binom{m-j-x-1}{i-1} \]

如果使用生成函数的话是可以解决的,但是想想细节太多,不想写。

不过值得注意的是,\(g_{i,j}=25g_{i,j-1}-24g_{i-1,j}\)。

所以换一种思路,这种东西也是可以矩阵乘法得到的,复杂度为 \(O(n^4\log k)\),过不去。

但是发现其实可以把这些所有的 \(i,j\) 放在同一个矩阵里面转移,所以复杂度变成 \(O(n^3\log k)\),可以过去。

注意有个细节,就是如果最后长度为奇数,那么最后一步转移不能是转移 3,把贡献减去即可。

代码

#include<bits/stdc++.h>
using namespace std;
using ll=long long;
const int N=2e2+10,M=N*1.5,mod=1e4+7;
int n,m,f[N][N][N];
char a[N];
int cnt,p[N],q[N];
struct matrix{
	int a[M][M];
	matrix(){
		memset(a,0,sizeof a);
	}
	matrix operator * (const matrix &x)const{
		matrix b;
		for(int k=0;k<cnt;k++)
			for(int i=0;i<=k;i++)
				for(int j=k;j<cnt;j++)
					b.a[i][j]=(b.a[i][j]+a[i][k]*x.a[k][j])%mod;
		return b;
	}
}base;
struct vec{
	int a[M];
	vec(){
		memset(a,0,sizeof a);
	}
	vec operator * (const matrix &x)const{
		vec b;
		for(int i=0;i<cnt;i++)
			for(int j=i;j<cnt;j++)
				b.a[j]=(b.a[j]+a[i]*x.a[i][j])%mod;
		return b;
	}
}ans;
int main(){
	freopen(".in","r",stdin);
	//freopen(".out","w",stdout);
	scanf("%s%d",a+1,&m),n=strlen(a+1);
	bool flag=(m+n)&1;
	m=(m+n+1)/2;
	f[1][n][0]=1;
	for(int l=1;l<=n;l++){
		for(int r=n;r>=l;r--){
			for(int i=0;i<=n;i++){
				if(!f[l][r][i])continue;
				if(a[l]==a[r]){
					(f[l+1][r-1][i+1]+=f[l][r][i])%=mod;
				}
				else{
					(f[l+1][r][i]+=f[l][r][i])%=mod;
					(f[l][r-1][i]+=f[l][r][i])%=mod;
				}
			}
		}
	}
	for(int i=0;i<=n;i++)p[i]=cnt++;
	for(int i=(n+1)/2;i>=0;i--)q[i]=cnt++;
	for(int l=1;l<=n+1;l++){
		for(int r=0;r<l;r++){
			for(int i=0;i<=n;i++){
				if(!f[l][r][i])continue;
				int j=(l-1)+(n-r)-i*2;
//				fprintf(stderr,"%d %d %d\n",i,j,f[l][r][i]);
				base.a[p[j]][q[i]]=(base.a[p[j]][q[i]]+f[l][r][i])%mod;
			}
		}
	}
	for(int i=1;i<=n;i++)base.a[p[i]][p[i]]=24;
	for(int i=0;i<n;i++)base.a[p[i]][p[i+1]]=1;
	for(int i=1;i<=(n+1)/2;i++)base.a[q[i]][q[i]]=25,base.a[q[i]][q[i-1]]=1;
	base.a[q[0]][q[0]]=26;
	ans.a[p[0]]=1;
	for(int y=m+1;y;base=base*base,y>>=1)if(y&1)ans=ans*base;
	int res=ans.a[q[0]];
	if(flag){
		base=matrix(),ans=vec();
		for(int l=0;l<=n+1;l++){
			for(int r=0;r<l;r++){
				for(int i=1;i<=n;i++){
					if(!f[l][r][i]||l!=r+1)continue;
					int j=(l-1)+(n-r)-i*2;
					base.a[p[j]][q[i]]=(base.a[p[j]][q[i]]+f[l][r][i])%mod;
				}
			}
		}
		for(int i=1;i<=n;i++)base.a[p[i]][p[i]]=24;
		for(int i=0;i<n;i++)base.a[p[i]][p[i+1]]=1;
		for(int i=1;i<=(n+1)/2;i++)base.a[q[i]][q[i]]=25,base.a[q[i]][q[i-1]]=1;
		base.a[q[0]][q[0]]=26;
		ans.a[p[0]]=1;
		for(int y=m;y;base=base*base,y>>=1)if(y&1)ans=ans*base;
		(res+=mod-ans.a[q[1]])%=mod;
	}
	cout<<res;
	return 0;
}

标签:matrix,Gift,--,zhengjun,int,端点,ans,转移
From: https://www.cnblogs.com/A-zjzj/p/17578139.html

相关文章

  • 分享一个关于网页上无法复制内容的处理技巧
      最近有朋友吐槽CSDN连个框选代码、复制功能都要登录才能使用,之前也有听到过类似的吐槽,这对于很多面向百度开发的程序员并不友好,于是想了个办法。  以CSDN为例,当需要复制CSDN上的内容是,可以先在开发者工具的控制台中输入下面的代码执行一遍,就可以使用框选复制了:  doc......
  • Linux系列---【如何确定自己服务器的磁盘是机械硬盘还是固态硬盘?】
    如何确定自己服务器的磁盘是机械磁盘还是固态硬盘?使用lsblk命令#这将列出所有磁盘设备的名称和rota(旋转)属性lsblk-d-oname,rota查看输出结果的第二列,rota属性的值:如果值为0,则表示该磁盘是固态硬盘(SSD)。如果值为1,则表示该磁盘是机械硬盘(HDD)。......
  • chinese 以前版本
    /*0.0.0梦开始的地方*/#include<bits/stdc++.h>#include<windows.h>#include<conio.h>usingnamespacestd;intpx=10,py=5,ma=0;stringjie[100005]={"这是墙壁,你不能通过","这是地面,上面似乎布满了灰尘"};stringdui[100005]={"墙","地&......
  • 我真的想知道,AI编译器中的IR是什么?
    随着深度学习的不断发展,AI模型结构在快速演化,底层计算硬件技术更是层出不穷,对于广大开发者来说不仅要考虑如何在复杂多变的场景下有效的将算力发挥出来,还要应对AI框架的持续迭代。AI编译器就成了应对以上问题广受关注的技术方向,让用户仅需专注于上层模型开发,降低手工优化性能......
  • JS中文件相关的知识(一):MIME类型
    不知道有没有同学和我一样,写代码时一遇到文件操作就犯怵,必须要先去把知识补一遍再说;对于Content-Type、responseType、ArrayBuffer、buffer、blob、file等这些词汇,心里问号一大堆,从来都没有真正区分清楚过;这样下去不是办法呀,真的猛士,应该敢于...一百次浮于表面,不如一次深入骨髓。......
  • 基础模型自监督预训练的数据之谜:大量数据究竟是福还是祸?
    前言 在自监督预训练中,是否数据越多越好?数据增广是否始终有效?本文转载自PaperWeekly作者|诺亚方舟实验室仅用于学术分享,若侵权请联系删除欢迎关注公众号CV技术指南,专注于计算机视觉的技术总结、最新技术跟踪、经典论文解读、CV招聘信息。CV各大方向专栏与各个部署框架最全......
  • HAOI2011 Problem b
    Problemblink做法:莫比乌斯反演。思路:对于给出的\(n\)个询问,每次求有多少个数对\((x,y)\),满足\(a\lex\leb\),\(c\ley\led\),且\(\gcd(x,y)=k\),\(\gcd(x,y)\)函数为\(x\)和\(y\)的最大公约数。我们设\[\operatorname{f}(n)=\sum\limits_{i=1}^x\sum\limit......
  • 获取上一个交易日接口
    实际实现交易日放在txt中上一交易日即当前交易日所在行的上一行数据python文件名lastjyr.pyimportosimportlinecacheimporttimedeflastjyr(): string=time.strftime("%Y%m%d") count=1 f=open('jyr.txt','r',encoding='utf8') forlineinf.re......
  • 磁盘分区
    使用fdisk命令对磁盘进行分区实验:磁盘分区1)新建磁盘sdb 一直下一步即可 使用命令lsblk查看分区情况,确认已经 2)分区设置一个主分区 查看是否完成分配使用命令w来保存配置 3)设置扩展分区 4)设置虚拟分区设置两个虚拟分区,会在扩展分区中存在,且序号将从......
  • Tool-Gitlab-502-端口占用
    安装Gitlab后,登陆报错502,端口占用卸载一、卸载GitLab及其依赖1、首先停止GitLab的运行命令:sudogitlab-ctlstop2、卸载GitLab:sudoapt-getremovegitlab-ce3、卸载GitLab依赖:sudoapt-getautoremove-y4、移除GitLab相关的文件:sudorm-rf/opt/gitlabtail日志sud......