MG loves string

 

 Accepts: 30    Submissions: 67

 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Others)

问题描述

MGMG是一个很忙碌的男孩子。今天他沉迷于这样一个问题:

对于一个长度为NN的由小写英文字母构成的随机字符串,当它进行一次变换,所有字符ii都会变成a[i]a[i]。

MGMG规定所有a[i]a[i]构成了2626个字母组成的排列。

MGMG现在需要知道这个随机串变换到自身的期望变换次数。请你输出期望答案乘上26^n26

​n

​​ 以后模 10000000071000000007 的结果。

MGMG认为这件事非常容易,不屑于用计算机解决,于是运用他高超的人类智慧开始进行计算。作为一名旁观者,你也想挑战MGMG智慧,请你写个程序,计算答案。

输入描述

第一行一个整数TT,代表数据组数(1 <=T<=101<=T<=10)。

接下来,对于每组数据——

第一行一个整数NN,表示给定的随机串长度(1<=N<=10000000001<=N<=1000000000)。

第二行2626个字母,表示a_ia

​i

​​ 序列

输出描述

对于每一组数据,输出一行。

显然,这个期望是一个实数。请你输出它乘上26^N26

​N

​​ 以后模 10000000071000000007 的结果

输入样例

2

2

abcdefghijklmnpqrstuvwxyzo

1

abcdefghijklmnopqrstuvwxyz

输出样例

5956

26

题解请参考这位大神:

https://blog.csdn.net/junior19/article/details/71599520

他的代码:

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<string>
#include<vector>
#include<deque>
#include<queue>
#include<algorithm>
#include<set>
#include<map>
#include<stack>
#include<ctime>
#include<string.h>
#include<math.h>
#include<list> using namespace std; #define ll long long
#define pii pair<int,int>
const int inf = 1e9 + 7; char to[27];
bool vis[27];
vector<pii>loop;//fist 循环节长度为first 有second个循环节长度为first void calLoop(){
for(int i=0;i<26;++i){
to[i]-='a';
}
map<int,int>mp;
fill(vis,vis+27,0);
loop.clear();
for(int i=0;i<26;++i){
if(vis[i]==0){
vis[i]=1;
int ans=1;
int x=to[i];
while(x!=i){
vis[x]=1;
++ans;
x=to[x];
}
mp[ans]+=1;
}
}
for(auto it=mp.begin();it!=mp.end();++it){
loop.push_back(*it);
}
} ll lcm(ll a,ll b){
ll tmp=__gcd(a,b);
return a/tmp*b;
} ll quickMulti(ll a,ll n){
ll ans=1;
ll t=a%inf;
while(n){
if(n&1){
ans=(ans*t)%inf;
}
t=(t*t)%inf;
n>>=1;
}
return ans;
} inline ll mod(ll x){
return (x%inf+inf)%inf;
} ll f(vector<int>&vec,int n){//容斥
ll ans=0;
int nv=vec.size();
for(int i=1,end=1<<nv;i<end;++i){
ll flag=-1;
int sum=0,num=0;
for(int j=0;j<nv;++j){
if((1<<j)&i){
sum+=vec[j];
++num;
}
}
if(num%2==nv%2){
flag=1;
}
ll t=quickMulti(sum,n);
ans=(ans+mod(flag*t))%inf;
}
return ans;
} ll slove(int n){
calLoop();
ll ans=0;
vector<int>vec;//当前选择了的循环节包含的字母数
for(int i=1,end=1<<loop.size();i<end;++i){//枚举循环节的状态
ll ans1=1;//选择了这几个循环节的贡献
vec.clear();
for(int j=0;j<loop.size();++j){
if((1<<j)&i){
ans1=lcm(ans1,loop[j].first);
vec.push_back(loop[j].first*loop[j].second);
}
}
if(vec.size()>n){
continue;
}
ll ans2=f(vec,n);
ans=(ans+ans1*ans2)%inf;
}
return ans;
} int main()
{
//freopen("in.txt","r",stdin);
//freopen("/home/lu/Documents/w.txt","w",stdout);
int T;
scanf("%d",&T);
while(T--){
int n;
scanf("%d%s",&n,to);
printf("%lld\n",slove(n));
}
return 0;
}

hdu 6021 MG loves string (一道容斥原理神题)(转)的更多相关文章

  1. hdu 6021 MG loves string

    MG loves string Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Others ...

  2. ●HDU 6021 MG loves string

    题链: http://acm.hdu.edu.cn/showproblem.php?pid=6021 题解: 题意:对于一个长度为 N的由小写英文字母构成的随机字符串,当它进行一次变换,所有字符 i ...

  3. 【HDU 6021】 MG loves string (枚举+容斥原理)

    MG loves string  Accepts: 30  Submissions: 67  Time Limit: 2000/1000 MS (Java/Others)  Memory Limit: ...

  4. MG loves string

    MG loves string Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Others ...

  5. hdu 6020 MG loves apple 恶心模拟

    题目链接:点击传送 MG loves apple Time Limit: 3000/1500 MS (Java/Others)    Memory Limit: 262144/262144 K (Ja ...

  6. hdu6021[BestCoder #93] MG loves string

    这场BC实在是有趣啊,T2是个没有什么算法但是细节坑的贪心+分类讨论乱搞,T3反而码起来很顺. 然后出现了T2过的人没有T3多的现象(T2:20人,T3:30人),而且T2的AC率是惨烈的不到3% ( ...

  7. 【CZY选讲·一道图论神题】

    题目描述 LYK有一张无向图G={V,E},这张无向图有n个点m条边组成.并且这是一张带权图,只有点权. LYK想把这个图删干净,它的方法是这样的.每次选择一个点,将它删掉,但删这个点是需要代价的 ...

  8. 一道背包神题-Petrozavodsk Winter-2018. Carnegie Mellon U Contest Problem I

    题目描述 有\(n\)个物品,每个物品有一个体积\(v_i\),背包容量\(s\).要求选一些物品恰好装满背包且物品个数最少,并在这样的方案中: (1)求出中位数最小的方案的中位数(\(k\)个元素的 ...

  9. 清北学堂模拟赛d2t1 一道图论神题(god)

    题目描述 LYK有一张无向图G={V,E},这张无向图有n个点m条边组成.并且这是一张带权图,只有点权. LYK想把这个图删干净,它的方法是这样的.每次选择一个点,将它删掉,但删这个点是需要代价的.假 ...

随机推荐

  1. typora的基本使用技巧汇总

    typora的基本使用技巧汇总 链接: https://www.jianshu.com/p/380005c8f104

  2. VS2017生成一个简单的DLL文件 和 LIB文件——C语言

    下面我们将用两种不同的姿势来用VS2017生成dll文件(动态库文件)和lib文件(静态库文件),这里以C语言为例,用最简单的例子,来让读者了解如何生成dll文件(动态库文件) 生成动态库文件 姿势一 ...

  3. 阿里云 负载均衡 HTTP转HTTPS

    一.相关文档 1.证书服务 2.简单路由-HTTP 协议变为 HTTPS 协议 二.阿里云操作界面 1.云盾证书服务管理控制台(查询CA证书服务) 2.负载均衡管理控制台 三.相关文档 1.Syman ...

  4. Visual Studio container tools require Docker to be running

    处理项目在生成时报错"Visual Studio container tools require Docker to be running" 最初win10上安装docker,项目 ...

  5. GoBelieve service部署问题整理

  6. TypeScript入门六:TypeScript的泛型

    泛型函数 泛型类 一.泛型函数 在泛型函数之前,先简单的描述一下泛型,将变量定义成泛型可以在使用变量时来决定它的类型.什么意思呢?假如现在有一个函数,可能出现参数和返回值出现多种情况的现象,只有在调用 ...

  7. C和指针课后问答题答案翻译

    第11章 动态内存分配 1.在你的系统中,你能够声明的静态数组最大长度能达到多少?使用动态内存分配,你最大能够获取的内存块有多大? 英文答案原文: This will vary from system ...

  8. Linux学习笔记(三)Linux常用命令:链接命令和文件查找命令

    一.链接命令 ln -s [原文件] [目标文件] (link) -s意为创建软连接 硬链接和软连接 硬链接的特点: (1)拥有相同的 i 结点和block块,可以看作是同一个文件 (2)可以通过 i ...

  9. Hadoop_14_MapReduce框架结构及其运行流程

    1.MapReduce原理篇 Mapreduce是一个分布式运算程序的编程框架,是用户开发“基于hadoop的数据分析应用”的核心框架: Mapreduce核心功能是将用户编写的业务逻辑代码和自带默认 ...

  10. Ubuntu系统---安装搜狗输入法

    Ubuntu_搜狗输入法 第一步: 1,下载搜狗输入法的安装包. 下载地址为:http://pinyin.sogou.com/linux/ ,如下图,要选择与自己系统位数一致的安装包,我的系统是64位 ...