文章目录
Sereja and Equality
简单组合数学+无脑前缀和优化dp计数。
f i , j f_{i,j} fi,j表示 i i i个数中逆序对数 ≤ j \le j ≤j的排列数
略卡常。。。
早年写过的代码:
#include<bits/stdc++.h>
#define N 505
#define mod 1000000007
using namespace std;
inline int read(){
int ans=0;
char ch=getchar();
while(!isdigit(ch))ch=getchar();
while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
return ans;
}
int C[N][N],f[N][N*N],T,n,k,ans,fac[N];
inline void init(){
fac[0]=1;
for(int i=1;i<=500;++i)fac[i]=1ll*fac[i-1]*i%mod;
for(int i=1;i<=500;++i)C[i][0]=C[i][i]=1,C[i][1]=i;
for(int i=2;i<=500;++i)for(int j=1;j<=i;++j){
C[i][j]=C[i-1][j]+C[i-1][j-1];
if(C[i][j]>=mod)C[i][j]-=mod;
}
int sum=0;
f[1][0]=1;
for(int i=2;i<=500;++i){
sum=0;
for(int j=0;j<=(i-1)*i/2;++j){
sum+=f[i-1][j];
if(sum>=mod)sum-=mod;
f[i][j]=sum;
if(j+1-i>=0)sum=(sum-f[i-1][j+1-i]+mod)%mod;
}
}
for(int i=2;i<=500;++i)for(int j=1;j<=(i-1)*i/2;++j){
f[i][j]+=f[i][j-1];
if(f[i][j]>=mod)f[i][j]-=mod;
}
}
int main(){
init(),T=read();
while(T--){
n=read(),k=read(),ans=0;
for(int i=1;i<=n;++i){
int tmp=1ll*C[n][i]*C[n][i]%mod*f[i][min(i*(i-1)/2,k)]%mod*fac[n-i]%mod*fac[n-i]%mod*(n-i+1)%mod;
ans+=tmp;
if(ans>=mod)ans-=mod;
}
printf("%d\n",ans);
}
return 0;
}
Division Game
传送门
不知道算不算博弈论,手玩一下推 % 4 = 0 , 1 , 2 , 3 \%4=0,1,2,3 %4=0,1,2,3的堆数的奇偶性对应的答案即可。
代码:
#include<bits/stdc++.h>
#define ri register int
using namespace std;
const int rlen=1<<18|1;
inline char gc(){
static char buf[rlen],*ib,*ob;
(ib==ob)&&(ob=(ib=buf)+fread(buf,1,rlen,stdin));
return ib==ob?-1:*ib++;
}
inline int read(){
int ans=0;
char ch=gc();
while(!isdigit(ch))ch=gc();
while(isdigit(ch))ans=((ans<<2)+ans<<1)+(ch^48),ch=gc();
return ans;
}
const int N=1e6+5;
int cnt[4],n;
int main(){
for(ri tt=read();tt;--tt){
n=read();
cnt[0]=cnt[1]=cnt[2]=cnt[3]=0;
for(ri i=1;i<=n;++i)cnt[read()%4]^=1;
int ans=cnt[0]*8+cnt[1]*4+cnt[2]*2+cnt[3];
if(ans==8)puts("Avnish");
else if(!ans)puts("Avnish");
else puts("Bishal");
}
return 0;
}
Reach Equilibrium
传送门
首先要首推一个结论:模长最长的向量长度不超过总长的一半。
然后从低维向高维推找规律即可。
然后官方题解给出了证明
代码:
#include<bits/stdc++.h>
#define ll long long
#define mod 1000000007
using namespace std;
inline ll ksm(ll x,ll p){
ll ret=1;
while(p){
if(p&1)ret=ret*x%mod;
x=x*x%mod;
p>>=1;
}
return ret;
}
int main(){
ll n,k;
scanf("%lld%lld",&n,&k);
printf("%lld",(ksm(2,n-1)-n+mod)%mod*ksm(2,(n-1)*(mod-2))%mod);
return 0;
}
Chef and Soccer
传送门
这题我用了一个挺复杂的 d p dp dp,貌似可以更简单???然而因为水过去了就懒得看其它做法了
设 f i , 0 / 1 / 2 , 0 / 1 , 0 / 1 f_{i,0/1/2,0/1,0/1} fi,0/1/2,0/1,0/1表示走到 i i i点,这一步向右走 0 / 1 / 2 0/1/2 0/1/2格( 0 0 0格表示在这里结束),第 i + 1 i+1 i+1个格子是否走过,第 i + 2 i+2 i+2个格子是否走过的档案数。
设 g i , 0 / 1 / 2 g_{i,0/1/2} gi,0/1/2表示走到 i i i点,从这一步开始一直向左走且这一步走了 0 / 1 / 2 0/1/2 0/1/2步的方案数。
然后简单讨论一下即可。
代码:
#include<bits/stdc++.h>
#define ri register int
using namespace std;
inline int read(){
int ans=0;
char ch=