把博客园图标替换成自己的图标
把博客园图标替换成自己的图标end

CF1628D Game on Sum

题面传送门
听说和THUPC的切切糕是一个题目?
首先考虑暴力dp:设\(dp_{i,j}\)表示做了\(i\)次,正了\(j\)次的双方最优策略情况下的答案。
显然有状态转移方程:\(dp_{i,j}=\max\limits_{0\leq d\leq k}{\min(dp_{i-1,j}-d,dp_{i-1,j-1}+d)}\)
观察发现其实转移方程就是\(dp_{i,j}=\frac{dp_{i-1,j}+dp_{i-1,j-1}}{2}\)
但是特别的有\(dp_{i,i}=ik\),就可以写出\(O(n^2)\)的暴力。
然后看上去就特别像一个坐标系走路的问题。所以可以考虑组合意义。
发现一旦平走了一次就永远走不到\(y=x\)这条直线上了,所以可以枚举平走出去的点,然后计算贡献即可。
时间复杂度\(O(n\log p)\)
code:

#include<bits/stdc++.h>
#define I inline
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define abs(x) ((x)>0?(x):-(x))
#define re register
#define RI re int
#define ll long long
#define db double
#define lb long db
#define N (1000000+5)
#define M (N*200)+5
#define mod 1000000007
#define Mod (mod-1)
#define eps (1e-9)
#define U unsigned int
#define it iterator
#define Gc() getchar() 
#define Me(x,y) memset(x,y,sizeof(x))
#define Mc(x,y) memcpy(x,y,sizeof(x))
#define d(x,y) (n*(x-1)+(y))
#define R(n) (rand()*rand()%(n)+1)
#define Pc(x) putchar(x)
#define LB lower_bound
#define UB upper_bound
#define PB push_back
using namespace std;
int n,m,k,T;ll Ans,frc[N+5],Inv[N+5];const int Inv2=(mod+1)/2;
I ll mpow(ll x,int y=mod-2){ll Ans=1;while(y) y&1&&(Ans=Ans*x%mod),y>>=1,x=x*x%mod;return Ans;}
I ll C(int x,int y){return frc[x]*Inv[y]%mod*Inv[x-y]%mod;}
I void Solve(){
	RI i,j;scanf("%d%d%d",&n,&m,&k);for(frc[0]=Inv[0]=i=1;i<=n;i++)frc[i]=frc[i-1]*i%mod,Inv[i]=mpow(frc[i]);if(n==m){printf("%lld\n",1ll*n*k%mod);return;}
	Ans=0;for(i=1;i<=min(m,n-1);i++) Ans+=C(n-i-1,m-i)*k%mod*i%mod*mpow(Inv2,n-i)%mod;printf("%lld\n",Ans%mod);
}
int main(){
//freopen("1.in","r",stdin);
	RI i;scanf("%d",&T);while(T--) Solve();
}
posted @ 2022-03-04 19:21  275307894a  阅读(34)  评论(0编辑  收藏  举报
浏览器标题切换
浏览器标题切换end