2018年10月31日提高组

前言

爆零的恐惧

---

JZOJ 5931 冒泡排序

题目

对于一个等概率随机的长度为$n$的排列，期望最少交换次数是多少

---

分析

对于n的排列，插到后面不影响的有$(n-1)!$种排列，那么增加1的就有$n!-(n-1)!$种，也就是$f[n](最少交换次数)=f[n-1]\ast (n-1)!+n!-(n-1)!$，除以$n!$求逆元即可

---

代码

#include <cstdio>
#define rr register
using namespace std;
const int N=10000000,mod=998244353;
int inv[N+5],f[N+5]={0,0,1},jc=2;
inline signed iut(){
	rr int ans=0; rr char c=getchar();
	while (c<48||c>57) c=getchar();
	while (c>47&&c<58) ans=(ans<<3)+(ans<<1)+c-48,c=getchar();
	return ans;
}
void print(int ans){
	if (ans>9) print(ans/10);
	putchar(ans%10+48);
}
inline signed ksm(long long x,long long y){
	rr long long ans=1ll;
	while (y){
		if (y&1) ans=ans*x%mod;
		x=x*x%mod; y>>=1;
	}
	return ans;
}
signed main(){
	for (rr long long i=3;i<=N;++i)
        f[i]=(i*f[i-1]+(i-1)*jc+mod)%mod,jc=i*jc%mod;
    inv[N]=ksm(jc,mod-2);
    for (rr long long i=N;i>=1;--i) inv[i-1]=i*inv[i]%mod;//倒着求逆元
	for (rr int t=iut();t;--t){
		rr int x=iut();
		print(1ll*f[x]*inv[x]%mod); putchar(10);
	}
	return 0;
}

---

JZOJ 5925 naive的瓶子

题目

有$n$个瓶子，排成一排。第$i$个瓶子的颜色为$ci$，每次操作可以选择两个相邻的瓶子，消耗他们颜色的数值乘积的代价将其中一个瓶子的颜色变成另一个瓶子的颜色。让所有瓶子的颜色都一样，操作次数不限，但要使得操作的总代价最小。

---

分析

那么在一个区间内如何变成同一种颜色，有两种方法，直接变和先变成最小的颜色，所以说可以枚举每一个瓶子颜色,然后$O(n)$判断即可

---

代码

#include <cstdio>
#define rr register
#define min(a,b) ((a)<(b))?(a):(b)
using namespace std;
struct rec{long long x[301];}a,b;
int n; long long ans;
inline signed iut(){
	rr int ans=0; rr char c=getchar();
	while (c<48||c>57) c=getchar();
	while (c>47&&c<58) ans=(ans<<3)+(ans<<1)+c-48,c=getchar();
	return ans;
}
void print(long long ans){
	if (ans>9) print(ans/10);
	putchar(ans%10+48);
}
void sol(int t){
	rr long long rest=0,sum=0; rec b=a;
	for (rr int i=1;i<n;++i)
	if (b.x[i]!=t&&b.x[i+1]!=t&&b.x[i]>b.x[i+1]&&b.x[i]*b.x[i+1]+b.x[i+1]*t<b.x[i]*t)//也就是说左边比右边大且直接合并比变颜色合并更好
	    rest+=b.x[i]*b.x[i+1],b.x[i]=b.x[i+1];
	for (rr int i=2;i<=n;++i)
	if (b.x[i]!=t&&b.x[i-1]!=t&&b.x[i]>b.x[i-1]&&b.x[i]*b.x[i-1]+b.x[i-1]*t<b.x[i]*t)//右边比左边大且直接合并比变颜色合并更好
	    rest+=b.x[i]*b.x[i-1],b.x[i]=b.x[i-1];
	for (rr int i=1;i<=n;++i) if (t!=b.x[i]) sum+=b.x[i];
	ans=min(ans,rest+sum*t);
}
signed main(){
	for (rr int t=iut();t;--t){
		n=iut(),ans=1e18;
		for (rr int i=1;i<=n;++i) a.x[i]=iut();
		for (rr int i=1;i<=n;++i) sol(a.x[i]);
		if (ans) print(ans); else putchar(48); putchar(10);
	}
	return 0;
}

---

JZOJ 5922 sequence

题目

给你一个长度为$n$的数列，开始时每个数字都是0。
有$m$个操作，每个操作形如$(l,r,k)$，对于区间$[l,r]$中的一个位置$i$，增加$C(i+k-l,k)$，求最终数列的每一项，对$10^9+7$取模。

---

分析

首先，暴力40我就不多说了
直接讲正解，根据$c(i,j)=c(i-1,j-1)+c(i-1,j)$，可以发现可以把k阶的前缀和分开，最后再通过杨辉三角全部加起来，then

---

代码

#include <cstdio>
#define rr register
using namespace std;
const int N=500001,mod=1000000007;
int c[N][22],n,f[N+1][22],m;
inline signed iut(){
	rr int ans=0; rr char c=getchar();
	while (c<48||c>57) c=getchar();
	while (c>47&&c<58) ans=(ans<<3)+(ans<<1)+c-48,c=getchar();
	return ans;
}
void print(int ans){
	if (ans>9) print(ans/10);
	putchar(ans%10+48);
}
signed main(){
	c[0][0]=1; n=iut();
	for (rr int i=1;i<N;++i)
	for (rr int j=0;j<21;++j) c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;//组合数不多说
    for (rr int m=iut();m;--m){
    	rr int l=iut(),r=iut(),k=iut();
    	for (rr int j=0;j<=k;++j) f[l][j]=(f[l][j]+c[k][k-j])%mod,
		f[r+1][j]=(f[r+1][j]-c[k+r-l+1][k-j]+mod)%mod;//差分不多说
	}
	for (rr int i=1;i<n;++i)
	for (rr int j=21;j>=1;--j) f[i+1][j-1]=(1ll*f[i+1][j-1]+f[i][j-1]+f[i][j]+mod)%mod;//最后再用杨辉三角合起来
	for (rr int i=1;i<=n;++i,putchar(10)) if (f[i][0]) print(f[i][0]); else putchar(48);
	return 0;
}