CF93D Flags

我的blog

题目链接：CF93D Flags

洛谷翻译：CF93D Flags

\[solution \]

我们先不考虑条件 \(4\) 。

设满足条件 \(1,2,3\) 时长度在 \([1,x]\) 之间的串的个数 \(g(x)\)

将 \(x\) 按奇偶性划分，通过简单的计算可以求出

\[g(x)=\left\{\begin{matrix}19\times3^{\frac{x}{2}-1}-7 & (x\equiv0\pmod2) \\ 11\times3^{\frac{x-1}2}-7 & (x\equiv1\pmod2) \end{matrix}\right. \]

我们考虑一下加上条件 \(4\) 第一反应为将答案除 \(2\) ，但是这种处理时不正确的，一个回文串本来只会被计算一次却被除 \(2\) 了，设满足题目所给的所有条件时长度在 \([1,x]\) 之间的串的个数 \(f(x)\) 。

因为 \(1,2,3\) 时长度在 \([1,x]\) 之间的回文串的个数 \(g(\left \lceil \frac{x}{2} \right \rceil)\) 。

所以 \(f(x)=\frac{g(x)+g(\left \lceil \frac{x}{2} \right \rceil)}{2}\)

答案为 \(f(r)-f(l-1)\)

\[code \]

#include <bits/stdc++.h>
#define lowbit(x) ((x)&(-(x)))
#define re register
#define ll long long
#define ull unsigned long long
#define rep(i,a,b) for(re int i=a;i<=b;++i)
#define per(i,a,b) for(re int i=a;i>=b;--i)
template<typename T>
inline void read(T&x)
{
	x=0;
	char s=(char)getchar();
	bool f=false;
	while(!(s>='0'&&s<='9'))
	{
		if(s=='-')
			f=true;
		s=(char)getchar();
	}
	while(s>='0'&&s<='9')
	{
		x=(x<<1)+(x<<3)+s-'0';
		s=(char)getchar();
	}
	if(f)
		x=(~x)+1;
}
template<typename T,typename ...T1>
inline void read(T&x,T1&...x1)
{
	read(x);
	read(x1...);
}
template<typename T>
inline void clear(T*array,int l,int r,int val)
{
	memset(&array[l],val,sizeof(T)*(r-l+1));
}
using namespace std;
const int mod=1e9+7;
const int inv2=(mod+1)/2;
inline ll quickpow(ll a,int b)
{
	ll res=1;
	while(b)
	{
		if(b&1)
			res=res*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return res;
}
inline ll f(int x)
{
	if(x==0)
		return 0;
	if(x&1)
		return (11*quickpow(3,x/2)-7)%mod;
	else
		return (19*quickpow(3,x/2-1)-7)%mod;
}
inline ll g(int x)
{
	return (f(x)+f((x+1)>>1))*inv2%mod;
}
int n;
signed main()
{
	int l,r;
	read(l,r);
	printf("%lld\n",(g(r)-g(l-1)+mod)%mod);
	return 0;
}