CF954F Runner‘s Problem

最新推荐文章于 2024-09-03 22:46:42 发布

C202044zxy

最新推荐文章于 2024-09-03 22:46:42 发布

阅读量148

点赞数

分类专栏：矩阵加速

本文链接：https://blog.csdn.net/C202044zxy/article/details/107912776

版权

矩阵加速专栏收录该内容

24 篇文章 1 订阅

订阅专栏

一、题目

点此看题

二、解法

本题显著的特点就是线段很长，但是段数很少，而且用 $d p$ 来做转移也很简单，这很容易联想到分段矩阵加速。

重点是线段的离散化，要把所有的端点插入一个数组中排序，再用类似差分的方法就可以标记一个段是否某一个限制线段所覆盖。值得注意的是我们的分段矩阵加速实际上是左开右闭的，所以一开始要插入 $l - 1$ （左端点）

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int M = 10005;
const int MOD = 1e9+7;
#define int long long
int read()
{
	int x=0,f=1;char c;
	while((c=getchar())<'0' || c>'9') {if(c=='-') f=-1;}
	while(c>='0' && c<='9') {x=(x<<3)+(x<<1)+(c^48);c=getchar();}
	return x*f;
}
int n,m,k,a[M],l[M],r[M],b[2*M],c[4][2*M],nw[4];
struct Matrix
{
	int n,m,a[5][5];
	Matrix() {n=m=0;memset(a,0,sizeof a);}
	void clear() {memset(a,0,sizeof a);}
	Matrix operator * (const Matrix &b) const
	{
		Matrix r;
		r.n=n;r.m=b.m;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++)
				for(int k=1;k<=b.m;k++)
					r.a[i][k]=(r.a[i][k]+a[i][j]*b.a[j][k])%MOD;
		return r;
	}
}A,F;
Matrix qkpow(Matrix a,int b)
{
	Matrix r;
	r.n=r.m=a.n;
	for(int i=1;i<=a.n;i++)
		r.a[i][i]=1;
	while(b>0)
	{
		if(b&1) r=r*a;
		a=a*a;
		b>>=1;
	}
	return r;
}
signed main()
{
	n=read();m=read();
	for(int i=1;i<=n;i++)
	{
		a[i]=read();l[i]=read();r[i]=read();
		b[++k]=l[i]-1;//important
		b[++k]=r[i];
	}
	b[++k]=1;b[++k]=m;
	sort(b+1,b+1+k);k=unique(b+1,b+1+k)-b-1;
	for(int i=1;i<=n;i++)
	{
		l[i]=lower_bound(b+1,b+1+k,l[i])-b;
		r[i]=lower_bound(b+1,b+1+k,r[i])-b;
		c[a[i]][l[i]]++;c[a[i]][r[i]+1]--;
	}
	A.n=A.m=F.n=3;F.m=1;
	F.a[2][1]=1;
	for(int i=2;i<=k;i++)
	{
		A.clear();
		A.a[1][1]=A.a[1][2]=A.a[2][1]=
		A.a[2][2]=A.a[2][3]=A.a[3][2]=A.a[3][3]=1;
		for(int j=1;j<=3;j++)
		{
			nw[j]+=c[j][i];
			if(nw[j]) A.a[j][1]=A.a[j][2]=A.a[j][3]=0;
		}
		A=qkpow(A,b[i]-b[i-1]);
		F=A*F;
	}
	printf("%lld\n",F.a[2][1]);
}