Problem M – Marblecoin(后缀数组,贪心)

首先:不得不说,后缀数组,永远滴神!
题意:有n个栈,每次可以从栈顶取一个元素,获得价值为v*365 t ^t t,其中v是这个元素的价值,是这个元素将拥有的天数,求最小总价值。
思路:这是个可以用贪心可以解决的问题,越早拿,系数越高,应该拿的元素的v也应该小,然而当两个数相同,就无法判断该从哪个栈中取。仔细想不难发现,这个问题的本质就是每次从字典序小的栈中拿出一个元素。然后就可以用后缀数组做了,两个栈之间用一个大的数隔开,每次取出后缀排名最高(rank最小的)的栈的栈顶。

收获:学到了维护多个串字典序关系的算法。

const int N = 6e5 + 5,M=22;
int s[N];
int sa[N], t[N << 1], t2[N << 1], c[N], height[N], rk[N];
int n, m;
void get_sa()
{
	//s从1开始,排名从1开始
	int i, *x = t, *y = t2, m = 310;
	for (i = 1;i <= m;i++)c[i] = 0;
	for (i = 1;i <= n;i++)c[x[i] = s[i]]++;
	for (int i = 2;i <= m;i++)c[i] += c[i - 1];
	for (i = n;i;i--)sa[c[x[i]]--] = i;
	for (int k = 1;k <= n;k <<= 1)
	{
		int num = 0;
		for (i = n - k + 1;i <= n;i++)y[++num] = i;
		for (i = 1;i <= n;i++)if (sa[i] > k)y[++num] = sa[i] - k;
		for (i = 1;i <= m;i++)c[i] = 0;
		for (i = 1;i <= n;i++)c[x[i]]++;
		for (i = 2;i <= m;i++)c[i] += c[i - 1];
		for (i = n;i;i--)sa[c[x[y[i]]]--] = y[i];
		swap(x, y);
		x[sa[1]] = 1;num = 1;
		for (int i = 2;i <= n;i++)
			x[sa[i]] = (y[sa[i]] == y[sa[i - 1]] && y[sa[i] + k] == y[sa[i - 1] + k]) ? num : ++num;
		if (num == n)break;
		m = num;
	}
}
void get_height()
{
	int k = 0;
	for (int i = 1;i <= n;i++)rk[sa[i]] = i;
	for (int i = 1;i <= n;i++)
	{
		if (rk[i] == 1)continue;
		if (k)k--;
		int j = sa[rk[i] - 1];
		while (j + k <= n && i + k <= n && s[i + k] == s[j + k])k++;
		height[rk[i]] = k;
	}
}
ll base[N + 1];
void init_base()
{
	base[0] = 1;
	f(i, 1, N)base[i] = base[i - 1] * 365 % mod;
}
struct node
{
	int idx, rk;
	bool operator<(const node b)const
	{
		if (this->rk == b.rk)return this->idx < b.idx;
		return this->rk < b.rk;
	}
};
int pos[N], ed[N];
int main()
{
	//freopen("in.txt", "r", stdin);
	init_base();n = 0;
	int t= in();
	int T = 0;
	f(i, 1,t )
	{
		int k = in();
		T += k;
		pos[i] = n + 1;//每个栈顶下标
		f(j, 1, k)
		{
			int x = in();
			s[++n] = x;
		}
		s[++n] = 305;//隔开每个栈顶,之后就可以识别每个栈的字典序
		ed[n] = 1;//标记每个栈结束位置
	}
	get_sa();
	get_height();
	//贪心,每次选每个栈字典序小的
	multiset<node> st;
	f(i, 1, t)st.insert({ pos[i],rk[pos[i]] });
	ll ans = 0,num=T;
	while (!st.empty())
	{
		auto it = *st.begin();
		st.erase(st.begin());
		ans = (ans+ base[T--]*s[it.idx]%mod)%mod;
		if (ed[it.idx + 1] == 0)
			st.insert({ it.idx + 1,rk[it.idx + 1] });
	}
	cout << ans << endl;
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值