703 删数 差分+倍增+动态规划 [代码源][namomo spring camp]每日一题div1

这肯定到银的程度了

//acmer mxc
#include<bits/stdc++.h>
#define mst(s,x) memset(s,x,sizeof(s));
#define sr(x) scanf("%lld",&x); 
#define sr2(a,b) scanf("%d%d",&a,&b); 
#define sr3(a,b,c) scanf("%d%d%d",&a,&b,&c); 
#define f(i,a,n) for(int i=a;i<=n;i++)
#define sc(x) cout<<#x<<" : "<<x<<endl;
#define hh cout<<endl;
#define pii pair<int,int>
#define pb(x) push_back(x)
#define mk(a,b) make_pair(a,b)
#define ls(x) x<<1
#define rs(x) x<<1|1
#define sz(s) (int)(s.size())
#define fi first
#define se second
#define inf 0x3f3f3f3f
#define ll __int128

using namespace std;
inline ll read()
{
	ll x=0,w=1; char c=getchar();
	while(c<'0'||c>'9') {if(c=='-') w=-1; c=getchar();}
	while(c<='9'&&c>='0') x=(x<<1)+(x<<3)+c-'0',c=getchar();
	return w==1?x:-x;
}
inline void write(ll x)
{
	if(x>=10) write(x/10);
	putchar(x%10+'0');
}
//------------------------------------------------
// #define int __int128
// #define int long long
#define LL long long
#define dd double
const int N=3e5+5;
int n;
int a[N],d[N],cnt[N];
int g[N][55];
int dp[N];
inline void init()
{
	for(int i=1;i<=n;i++)a[i]=d[i]=cnt[i]=0;
	for(int i=1;i<=n;i++)for(int j=0;j<53;j++)g[i][j]=-1;
	for(int i=1;i<=n;i++)dp[i]=inf;
}
inline void solve()
{
	scanf("%d",&n);
	init();
	for(int i=1;i<=n;i++)scanf("%d",&a[i]);
	for(int i=1;i<n;i++)d[i]=a[i+1]-a[i];//差分
	for(int i=1;i<n;i++)//只有二进制的时候才可能被后面的消去，0的时候特殊考虑
	{
		if(d[i]==0)d[i]=inf;
		else 
		{
			while(d[i]%2==0) cnt[i]++,d[i]>>=1;
		}
	}
	//从后往前考虑g[i][cnt[i]]表示从区间[i,g[i][cnt[i]]]可以得到一个2^{cnt[i]}的值
	//考虑每一位初始化都是区间[i.i+1]得到2^{cnt[i]}
	for(int i=n-1;i>=1;i--)
	{
		g[i][cnt[i]]=i+1;
		for(int j=cnt[i]+1;j<=53;j++)//当前是2^cnt[i]很显然我们只能合成比他大的 
		{
			if(g[i][j-1]==-1) break;//如果能合成2^j那么一定可以有2^{j-1}
			if(g[i][j-1]>=n) break;//合成2^{j-1}肯定不能超范围
			if(d[g[i][j-1]]!=d[i])break;//我们在之前把d[i]的数组预处理了，如果他们都不符合二进制的倍数关系肯定没法合呀！
			int k=g[i][j-1];
			g[i][j]=g[k][j-1];//合一个2^j需要两个2^j-1需要从在[i,g[i][j-1]]合一个再从[g[i][j-1],j[g[i][j-1]][j-1]]合一个
		}
	}
	//到现在我们就知道了如果g[i][j]不等于-1，那么就可以从[i,g[i][j]]合成一个数，这样就可以DP了
	//dp[i]表示前1~i-1的所有数字最小长度
	//那么就有dp[g[i][j]]=min(dp[i][j],dp[i]+1);
	dp[1]=1;
	for(int i=1;i<n;i++)
	{
		if(d[i]==inf&&d[i+1]==inf) //差分为0的单独考虑
			dp[i+1]=min(dp[i+1],dp[i]);
		else 
		{
			for(int j=0;j<=53;j++)
			{
				if(g[i][j]!=-1) dp[g[i][j]]=min(dp[i]+1,dp[g[i][j]]);
			}
		}
	}
	//我们要考虑所有的数全用那么ans就是dp[n]
	printf("%d\n",dp[n]);
	return;
}

//------------------------------------------------
signed main()
{
	clock_t c1=clock();
#ifdef LOCAL
    freopen("in.in","r",stdin);
    freopen("out.out","w",stdout);
#endif
//==================================
    int T=1;
	cin>>T;
	for(int i=1;i<=T;i++)
	{
		solve();
	}
//==================================
end:
    cerr<< endl << "Time used:" << clock() - c1 << endl;
	return 0;
}