牛客暑期多校训练营2020第3场

E Two Matchings

题意：
  给一个数组a，求排列p、q，满足p_i != q_i。p_pi=i，q_qi=i。(n为偶数) 求$\frac{1}{2}$${\sum }_{i=1}^n$a_I-a_pi与$\frac{1}{2}$${\sum }_{i=1}^n$a_I-a_qi的和的最小值。
思路：
  p_i != q_i即是p中均为二元环。
  n=4时，最小值为max-min；
  n=6时，最小值为max-min；
  n $\ge$ 8时，整块不会比切成4、6小块考虑更优。
  然后动态规划就可以了。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define LL long long
using namespace std;
const int N=2e5+10;
const LL inf=1e16;
LL a[N],f[N];
int cmp(LL a,LL b){return a>b;}
int main()
{
	int T,n;
	scanf("%d",&T);
	while(T)
	{
		T--;
		scanf("%d",&n);
		for(int i=1;i<=n;i++)	scanf("%lld",&a[i]);
		sort(a+1,a+n+1,cmp);
		f[4]=a[1]-a[4];
		f[6]=a[1]-a[6];
		for(int i=8;i<=n;i+=2)
		{
			f[i]=f[i-4]+a[i-4+1]-a[i];
			if(i-6>=4)	f[i]=min(f[i],f[i-6]+a[i-6+1]-a[i]);
		}
		printf("%lld\n",f[n]*2);
	}
	return 0;
 } 

H Sort the Strings Revision

题意：
  给一个串s，一个排列p，一个数组d。第i次操作为：将s_i-1的第p_i位换成d_i，求s₀、s₁、……s_n的排名。
思路：
  对于原串的第i位，假设它是在第k次操作中被改变的，如果：

1. d[k] > s[i] ，那么s_k，……s_n排在s₀……s_k-1后面
2. d[k] < s[i] ，那么s_k，……s_n排在s₀……s_k-1前面
3. d[k] > s[i] ，那就当什么都没发生过

  所以我们可以从第一位开始考虑，看它是在什么时候被改变的就可以了。但是，我们如何处理才能做到O(n)呢？按照笛卡尔树dfs序！
笛卡尔树

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define LL long long
#define int long long
using namespace std;
const int N = 2e6+10;
const LL base = 10000019;
const LL md = 1000000007;
LL p[N], d[N], c[N], ha[N];
int t[N], stk[N], l[N], r[N];
int num;
void init(int n)
{
	LL pseed, pa, pb, pmod;
	scanf("%lld%lld%lld%lld",&pseed,&pa,&pb,&pmod);
	for(int i = 0; i < n; i++)	p[i] = i;
	for(int i = 1; i < n; i++)
	{
		swap(p[pseed % (i + 1)], p[i]);
		pseed = (pseed * pa + pb) % pmod;
	}
	LL dseed, da, db, dmod;
	scanf("%lld%lld%lld%lld", &dseed, &da, &db, &dmod);
	for(int i = 0; i < n; i++)
	{
		d[i] = dseed % 10;
		dseed = (dseed * da + db) % dmod;
	}
}
void dfs(int x, int ll, int rr)
{
	if(ll>rr)	return;
	if(x==-1)
	{
		for(int i = ll; i <= rr; i++)	c[i] = num++;
		return;
	}
	if(d[x] > p[x] % 10)
	{
		dfs(l[x], ll, x);
		dfs(r[x], x+1, rr);
	} 
	else
	{
		dfs(r[x], x+1, rr);
		dfs(l[x], ll, x); 
	}
}
void solve(int n)
{
	int cnt=0, rt;
	for(int i = 0; i <= n; i++)	c[i] = 0;
	for(int i = 0; i < n; i++)
		if(d[i] != p[i] % 10)
			t[++cnt] = i;
	if(!cnt)	rt=-1;
	else
	{
		rt = t[1];
		stk[1] = t[1];
		int top = 1;
		l[t[1]] = r[t[1]] = -1;
		for(int i = 2; i <= cnt; i++)
		{
			l[t[i]] = r[t[i]] = -1;
			while(top && p[stk[top]] > p[t[i]])
			{
				l[t[i]] = stk[top];
				top--;
			} 
			if(top) r[stk[top]] = t[i];
			else rt = t[i];
			stk[++top] = t[i];
		}
	}
	num=0;
	dfs(rt, 0, n);
}
void print(int n)
{
	LL ans=0;
	for(int i = 0; i <= n; i++)	ans = (ans + c[i] * ha[i] % md) % md; 
	printf("%lld\n", ans);
}
signed main()
{
	int T, n;
	ha[0] = 1;
	for(int i = 1; i <= 2e6; i++)	ha[i] = ha[i-1] * base % md;
	scanf("%lld",&T);
	while(T)
	{
		T--;
		scanf("%lld",&n);
		init(n);
		solve(n);
		print(n);
	}
	return 0;
 }