POJ - 3734 Blocks 指数生成函数

传送门

文章目录

  • 题意:
  • 思路:

题意:

一段长度为nnn的序列,你有红黄蓝绿四种颜色的砖块,问你铺砖的方案数,每块砖长度为111,其中红黄颜色个数必须为偶数。

思路:

考虑多重集合排列数:
在这里插入图片描述
所以我们构造四种颜色的指数生成函数,并转换成封闭式:
f红(x)=f黄(x)=1+x22!+x44!+...=ex+e−x2f_{红}(x)=f_{黄}(x)=1+\frac{x^2}{2!}+\frac{x^4}{4!}+...=\frac{e^x+e^{-x}}{2}f(x)=f(x)=1+2!x2+4!x4+...=2ex+ex
f蓝(x)=f绿(x)=1+x1!+x22!+...=exf_{蓝}(x)=f_{绿}(x)=1+\frac{x}{1!}+\frac{x^2}{2!}+...=e^xf(x)=f绿(x)=1+1!x+2!x2+...=ex

将其乘起来得封闭式
G=e2x∗e2x+e−2x+24=e4x+2e2x+14G=e^{2x}*\frac{e^{2x}+e^{-2x}+2}{4}=\frac{e^{4x}+2e^{2x}+1}{4}G=e2x4e2x+e2x+2=4e4x+2e2x+1
再将封闭式展开
ex=∑i≥0xii!e^x=\sum_{i\ge0}\frac{x^i}{i!}ex=i0i!xi,得ekx=∑i≥0kixii!e^{kx}=\sum_{i\ge0}\frac{k^ix^i}{i!}ekx=i0i!kixi,即G=14+∑i≥0(4i+2i+1)∗xii!4G=\frac{1}{4}+\frac{\sum_{i\ge0}(4^i+2^{i+1})*\frac{x^i}{i!}}{4}G=41+4i0(4i+2i+1)i!xi,取指数为nnn的项,得系数4n+2n+14∗n!\frac{4^n+2^{n+1}}{4*n!}4n!4n+2n+1,可以将常数项忽略,再乘上多重集合排列数的分子n!n!n!4n+2n+14=4n−1+2n−1\frac{4^n+2^{n+1}}{4}=4^{n-1}+2^{n-1}44n+2n+1=4n1+2n1,答案即为4n−1+2n−14^{n-1}+2^{n-1}4n1+2n1,快速幂直接算一下即可。

// Problem: Blocks
// Contest: Virtual Judge - POJ
// URL: https://vjudge.net/problem/POJ-3734
// Memory Limit: 65 MB
// Time Limit: 1000 ms
// 
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=10007,INF=0x3f3f3f3f;
const double eps=1e-6;LL qmi(LL a,LL b) {LL ans=1;while(b) {if(b&1) ans=ans*a%mod;b>>=1;a=a*a%mod;}return ans%mod;
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);	int _; cin>>_;while(_--) {int n; cin>>n;cout<<((1ll*qmi(4,n-1)+qmi(2,n-1))%mod)<<endl;}return 0;
}
/**/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314739.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

理想化的DevOps团队里只需要有Dev就够了?

&#xff08;图片来源于网络&#xff09;几天前&#xff0c;本公众号发布的一篇译文列举了9种DevOps团队结构适用类型与7种反型&#xff08;点击查看原文&#xff09;。文章转发到朋友圈之后&#xff0c;很多DevOps同行留言&#xff08;吐槽&#xff09;了自己团队的现状&#…

hdu 1521 排列组合 多重集排列 + 指数生成函数

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 显然是多重集排列数&#xff0c;我们考虑构造指数生成函数&#xff0c;让后模拟一下多项式乘法即可啦。 由于存在分数&#xff0c;所以直接用doubledoubledouble即可。 //#pragma GCC opt…

架构杂谈《七》

容器VS虚拟机一、什么是虚拟机虚拟机&#xff08;Virtual Machine&#xff09;指通过软件模拟的具有完整硬件系统功能的、运行在一个完全隔离环境中的完整计算机系统。虚拟系统通过生成现有操作系统的全新虚拟镜像&#xff0c;它具有真实windows系统完全一样的功能&#xff0c;…

#3027. [Ceoi2004]Sweet 生成函数 + 组合数学

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 这个题显然可以容斥来写&#xff0c;刚学生成函数就来水一下。 对于每一堆iii我们写出其生成函数Fi(x)∑k0mi(1xx2...xmi)1−x1mi1−xF_i(x)\sum_{k0}^{m_i}(1xx^2...x^{m_i})\frac{1-x^{1…

坚持一个好习惯该有多难?

引子这段时间以来我都在坚持写博客&#xff0c;但是更新得比较少&#xff0c;每周才能一更&#xff0c;不过好歹也勉勉强强也能坚持下去了。我们社区还组建了一个博客互勉群&#xff0c;不过除了个别人能经常打卡外&#xff0c;大部分人都是混个脸熟&#xff0c;凑个热闹。在这…

#3771. Triple 生成函数 + FFT + 容斥

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 注意到这个题是求若干个数的组合数&#xff0c;(a,b),(b,a)(a,b),(b,a)(a,b),(b,a)视为一种方案&#xff0c;所以我们考虑生成一个普通型生成函数。 考虑到每个数只能选一次&#xff0c;但…

gRPC in ASP.NET Core 3.0 -- 前言

现如今微服务很流行&#xff0c;而微服务很有可能是使用不同语言进行构建的。而微服务之间通常需要相互通信&#xff0c;所以微服务之间必须在以下几个方面达成共识&#xff1a;需要使用某种API数据格式错误的模式负载均衡。。。现在最流行的一种API风格可能是REST&#xff0c;…

P3246 [HNOI2016]序列 莫队 + ST表 + 单调栈

传送门 文章目录题意&#xff1a;思路&#xff1a;Update题意&#xff1a; 思路&#xff1a; 比较神奇的一个题&#xff0c;这里先介绍莫队的离线解法。 不难发现&#xff0c;用莫队来做最大的难点就是在进行区间移动的时候如何快速计算贡献。 比如[l,r]−>[l,r1][l,r]->…

WPF 渲染原理

对于开发者来说&#xff0c;WPF 中最主要的知识点就是渲染&#xff0c;因为 WPF 是一个界面框架。想用一篇博客就能告诉大家完整的 WPF 渲染原理是不可能的。本文向大家介绍从开发者执行绘图指令到在屏幕显示的过程。本文是从一个全局的角度来看渲染的过程&#xff0c;在本文之…

HDU - 6989 Didn‘t I Say to Make My Abilities Average in the Next Life?! 莫队/单调栈 + 线段树/ST表在线

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 考虑将贡献分开来算&#xff0c;先计算最大值&#xff0c;再算个最小值&#xff0c;之后答案就是((maxmin)/2)/(len∗(len1)/2)((maxmin)/2)/(len*(len1)/2)((maxmin)/2)/(len∗(len1)/2)。…

ASP.NET Core 中的管道机制

首先&#xff0c;很感谢在上篇文章 C# 管道式编程 中给我有小额捐助和点赞的朋友们&#xff0c;感谢你们的支持与肯定。希望我的每一次分享都能让彼此获得一些收获&#xff0c;当然如果我有些地方叙述的不正确或不当&#xff0c;还请不客气的指出。好了&#xff0c;下面进入正文…

微软宣布SQL Server 2019免费支持Java

在 2018 年 9 月的时候&#xff0c;微软就宣布与领先的 Java 开源贡献者和发行商 Azul Systems 建立新的合作伙伴关系。这一关系允许所有 Azure 客户在微软和 Azul 联合提供的支持下&#xff0c;免费使用 Azul 的 Zulu for Azure-Enterprise Java 发行版 。最近&#xff0c;微软…

P4231 三步必杀 二次差分

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 考虑给[2,6][2,6][2,6]加上s2,e10s2,e10s2,e10的等差数列&#xff0c;变成2,4,6,8,102,4,6,8,102,4,6,8,10&#xff0c;考虑差分数组2,2,2,2,2,−102,2,2,2,2,-102,2,2,2,2,−10&#xff0…

创新的迷思——2019 Microsoft Hackathon 苏州站回顾

每年&#xff0c;微软在7月份都会在内部举办一次黑客松&#xff08;Hachathon&#xff09;活动&#xff0c;所有的员工&#xff0c;包括实习生都可以自由组队参赛&#xff0c;在为期一周的时间内完成构想、设计和开发的一系列过程&#xff0c;并且在最后一天进行项目路演和评选…

P4062 [Code+#1]Yazid 的新生舞会 树状数组维护三阶差分

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个序列aaa&#xff0c;让你求有多少个子区间满足存在一个数是这个区间的绝对众数&#xff0c;绝对众数指该数在区间内出现的次数严格大于r−l12\frac{r-l1}{2}2r−l1​。 n≤5e5,0≤ai≤n−1n\le5e5,0…

一句话概括4本管理著作

之前618活动的时候买了不少书&#xff0c;最近一段时间集中看了四本管理相关的书籍&#xff0c;收获颇丰。在这里分享给大家。如果你不是管理者也没关系&#xff0c;可以换个视角来看看管理者眼中的世界&#xff0c;毕竟&#xff0c;“赤兔马”也需要让“伯乐”看到才行嘛。虽然…

[译].Net中的内存-什么分配在了哪里

原文链接&#xff1a;https://jonskeet.uk/csharp/memory.html人们在理解值类型和引用类型之间的差异时因为“值类型在栈上分配&#xff0c;引用类型在堆上分配”这句话造成了很多混乱。这完全是不对的&#xff0c;本文试图澄清这个问题。变量中有什么&#xff1f;理解.NET中内…

上车时机已到--.NETCore是适应时代发展的雄鹰利剑

要起飞了.NET Core 3.0-prevew7&#xff1a;https://dotnet.microsoft.com/download/dotnet-core/3.0随着 .NET Core 3.0-prevew7 的发布&#xff0c;开源社区的一支重要力量重要即将起飞&#xff1b;官方指出&#xff0c;该预览版为可能为最终版本&#xff0c;在很长一段时间内…

HDU - 7028 Decomposition 无向完全图构造欧拉回路

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一张无向完全图&#xff0c;让你构造kkk个长度分别为lil_ili​的路径&#xff0c;这些路径不相交&#xff0c;且∑lin∗(n−1)2\sum l_i\frac{n*(n-1)}{2}∑li​2n∗(n−1)​。 n≤1000,nmod21n\le1000,…

「Azure」数据分析师有理由爱Azure之二-立即申请帐号开始学习之旅

目前关于Azure的学习资料不多&#xff0c;除了官方的文档和Microsoft Learn频道外&#xff0c;几乎没有什么中文性资料可学习&#xff0c;就算有&#xff0c;也是以IT的思维方式来展开介绍&#xff0c;对没有IT背景的数据分析师来说&#xff0c;非常难于适应。本篇给大家带来一…