loj6077. 「2017 山东一轮集训 Day7」逆序对

题目描述：

loj

题解：

容斥+生成函数。

考虑加入的第$i$个元素对结果的贡献是$[0,i-1]$，我们可以列出生成函数。

长这样：$(1)*(1+x)*(1+x+x^2)*……*(1+x+x^2+……+x^{n-1})=\frac{\prod_{i=1}^{n}1-x^i}{(1-x)^n}$

把分母提出来：$\frac{1}{(1-x)^n} = (1+x+x^2+……)^n = \sum_{i=0}^{k} C_{i+n-1}^{n-1}$，日常小球放盒。

现在还剩$\prod_{i=1}^n 1-x^i$，可以考虑将该式理解为从$1$到$n$选$i$个数，总和为$j$，对该项系数贡献为$(-1)^i$。

由于从$1$到$n$不重复，可以发现$\sum_{i=1}^{447}>100000$，那么$i \le 447$。

那么就可以$dp$了。状态为$f[i][j]$，表示当前选了$i$个数总和为$j$，且最后一项不大于$n$的方案数。

要求选数不重复怎么办？

考虑将其构造成一个上升序列。我们用枚举差值的思想，保证前后差值大于$0$。

转移有三种：

1.将$i$个数集体+1，此时$f[i][j]+=f[i][j-i]$；

2.将$i$个数集体+1再在最前面加入一个1，此时$f[i][j]+=f[i-1][j-i]$；

3.考虑我们每次都让每个数+1，所以当$j>n$时会出现最后一项为$n+1$的情况，此时$f[i][j]-=f[i-1][j-n-1]$；

这样转移就可以了。

最后卷积卷出第$k$位的值就可以了。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int MOD = 1000000007;
const int N = 100050;
const int M = 450;
int fastpow(int x,int y)
{int ret = 1;while(y){if(y&1)ret=1ll*ret*x%MOD;x=1ll*x*x%MOD;y>>=1;}return ret;
}
template<typename T>
inline void Mod(T&x){if(x>=MOD)x-=MOD;}
int n,m;
int jc[N<<1],jny[N<<1],f[M][N];
void init()
{jc[0] = 1;for(int i=1;i<=n+m;i++)jc[i]=1ll*jc[i-1]*i%MOD;jny[n+m] = fastpow(jc[n+m],MOD-2);for(int i=n+m;i;i--)jny[i-1]=1ll*jny[i]*i%MOD;
}
int C(int x,int y){return 1ll*jc[x]*jny[y]%MOD*jny[x-y]%MOD;}int main()
{scanf("%d%d",&n,&m);init();f[0][0] = 1;for(int i=1;i<M;i++)for(int j=i;j<=m;j++){Mod(f[i][j] = f[i-1][j-i]+f[i][j-i]);if(j>n)Mod(f[i][j]+=MOD-f[i-1][j-n-1]);}int ans = 0;for(int i=0;i<=m;i++){int tmp = 0;for(int j=0;j<M;j++)if(j&1)Mod(tmp+=MOD-f[j][i]);else Mod(tmp+=f[j][i]);Mod(ans+=1ll*tmp*C(n+m-i-1,n-1)%MOD);}printf("%d\n",ans);return 0;
}

View Code

转载于:https://www.cnblogs.com/LiGuanlin1124/p/11136070.html

loj6077. 「2017 山东一轮集训 Day7」逆序对相关推荐

Loj #6077. 「2017 山东一轮集训 Day7」逆序对
Loj #6077. 「2017 山东一轮集训 Day7」逆序对 Solution 令fi,jf_{i,j}fi,j表示前iii个数产生jjj个逆序对的方案数,每次考虑把i+1i+1i+1加入,有i ...
容斥问卷调查反馈——Co-prime，Character Encoding，Tree and Constraints，「2017 山东一轮集训 Day7」逆序对
文章目录 Co-prime source solution code Character Encoding source solution code Tree and Constraints sour ...
LOJ6079「2017 山东一轮集训 Day7」养猫
养ImmortalCO k可重区间问题的增强版:有上下界! 直接都选择s[i],然后再把一些调整到e[i] 考虑通过最大流的"最大",使得至少每k个有me个e, 通过最大流的&q ...
#6073. 「2017 山东一轮集训 Day5」距离（树链剖分 + 永久标记主席树）
#6073. 「2017 山东一轮集训 Day5」距离给定一颗有nnn个节点带边权的树,以及一个排列ppp,path(u,v)path(u, v)path(u,v)为u,vu, vu,v路径上的点集 ...
LOJ#6072. 「2017 山东一轮集训 Day5」苹果树解题报告
LOJ#6072. 「2017 山东一轮集训 Day5」苹果树解题报告好苹果会组成连通块,整棵树的权值为 ∑ i = 1 n c i [ c i ≥ 0 ] [ s i z n u m ( c i ...
LOJ#6074. 「2017 山东一轮集训 Day6」子序列
LOJ#6074. 「2017 山东一轮集训 Day6」子序列先考虑全局询问怎么做,设 f ( i , c ) f(i,c) f(i,c) 表示在 S 1 ⋯ i S_{1\cdots i} S1⋯ ...
LOJ#6103. 「2017 山东二轮集训 Day2」第一题解题报告
LOJ#6103. 「2017 山东二轮集训 Day2」第一题解题报告前置知识:闭区间上的连续函数的零点存在性定理: 我们定义这样的函数: 定义域为 [ l , r ] ∩ Z [l,r]\cap ...
「2017 山东一轮集训 Day5」距离
/* 写完开店再写这个题目顿时神清气爽, 腰也不疼了, 眼也不花了首先考虑将询问拆开, 就是查询一些到根的链和点k的关系根据我们开店的结论, 一个点集到一个定点的距离和可以分三部分算那么就很简单了吧 ...
[LOJ#6068]. 「2017 山东一轮集训 Day4」棋盘[费用流]
题意题目链接分析考虑每个棋子对对应的横向纵向的极大区间的影响:记之前这个区间中的点数为 $x$ ,那么此次多配对的数量即 $x$ . 考虑费用流,\(S\rightarrow 横向区间 ...
「2017 山东一轮集训 Day2」Pair (霍尔定理+线段树)
题目描述给出一个长度为的数列和一个长度为的数列 ,求有多少个长度为的连续子数列能与匹配. 两个数列可以匹配,当且仅当存在一种方案,使两个数列中的数可以两两配对,两个数可以配对 ...

loj6077. 「2017 山东一轮集训 Day7」逆序对

loj6077. 「2017 山东一轮集训 Day7」逆序对相关推荐

最新文章

热门文章