前言
依旧是白嫖账号,只打了一些题/kk
正题
1002 Buying Snacks
题目大意
nnn个物品,每个可以买一次也可以不买,如果买需要选择1/21/21/2块钱的,然后也可以相邻两个一起买并且减少一块的花销,求恰好用掉mmm块钱的方案。
1≤n≤109,1≤m≤20000,1≤T≤51\leq n\leq 10^9,1\leq m\leq 20000,1\leq T\leq 51≤n≤109,1≤m≤20000,1≤T≤5
解题思路
设fi,jf_{i,j}fi,j表示iii个物品花jjj块钱的方案那么有
fi,j=fi−1,j+fi−1,j−1+fi−1,j−2+fi−2,j−1+2×fi−2,j−2+fi−2,j−3f_{i,j}=f_{i-1,j}+f_{i-1,j-1}+f_{i-1,j-2}+f_{i-2,j-1}+2\times f_{i-2,j-2}+f_{i-2,j-3}fi,j=fi−1,j+fi−1,j−1+fi−1,j−2+fi−2,j−1+2×fi−2,j−2+fi−2,j−3
然后化成生成函数就是
Fi(x)=(1+x+x2)Fi−1(x)+(x+2x2+x3)Fi−2(x)F_i(x)=(1+x+x^2)F_{i-1}(x)+(x+2x^2+x^3)F_{i-2}(x)Fi(x)=(1+x+x2)Fi−1(x)+(x+2x2+x3)Fi−2(x)
和之前一道倍增FFTFFTFFT很像,然后考虑分割位置填两个就是F1(x)=1+x+x2F_{1}(x)=1+x+x^2F1(x)=1+x+x2,然后
Fa+b(x)=Fa(x)Fb(x)+(x+2x2+x3)Fa−1(x)Fb−1(x)F_{a+b}(x)=F_{a}(x)F_b(x)+(x+2x^2+x^3)F_{a-1}(x)F_{b-1}(x)Fa+b(x)=Fa(x)Fb(x)+(x+2x2+x3)Fa−1(x)Fb−1(x)
上倍增FFTFFTFFT就好了。
code
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=1<<16,P=998244353;
ll T,n,k,m,r[N],f[3][N],t[3][N],g[2][N];
void fm(ll &x){x+=x>>31&P;}
ll power(ll x,ll b){ll ans=1;while(b){if(b&1)ans=ans*x%P;x=x*x%P;b>>=1;}return ans;
}
void NTT(ll *f,ll op){for(ll i=0;i<n;i++)if(i<r[i])swap(f[i],f[r[i]]);for(ll p=2;p<=n;p<<=1){ll tmp=power(3,(P-1)/p),len=p>>1;if(op==-1)tmp=power(tmp,P-2);for(ll k=0;k<n;k+=p){ll buf=1;for(ll i=k,tt;i<(k|len);i++){tt=1ll*buf*f[i|len]%P;fm(f[i|len]=f[i]-tt);fm(f[i]=f[i]+tt-P);buf=1ll*buf*tmp%P;}}}if(op==-1){ll invn=power(n,P-2);for(ll i=0;i<n;i++)f[i]=1ll*f[i]*invn%P;}return;
}
void print(ll x)
{if(x>9)print(x/10);putchar(x%10+'0');return;}
signed main()
{scanf("%lld",&T);while(T--){memset(f,0,sizeof(f));memset(g,0,sizeof(g));memset(t,0,sizeof(t)); scanf("%lld%lld",&m,&k);k++;n=1;while(n<(k*2))n<<=1;for(ll i=0;i<n;i++)r[i]=(r[i>>1]>>1)|((i&1)?(n>>1):0);f[0][0]=f[0][1]=f[0][2]=f[1][0]=g[0][0]=1;for(ll d=1;d<=m;d<<=1){if(m&d){for(ll j=0;j<3;j++){for(ll i=0;i<n;i++)t[j][i]=(i<k)?f[j][i]:0;NTT(t[j],1);}NTT(g[0],1);NTT(g[1],1);for(ll i=0;i<n;i++){ll b0=g[0][i],b1=g[1][i];g[0][i]=1ll*b0*t[0][i]%P;g[1][i]=1ll*b0*t[1][i]%P;t[0][i]=1ll*t[1][i]*b1%P;t[1][i]=1ll*t[2][i]*b1%P;}NTT(g[0],-1);NTT(g[1],-1);NTT(t[0],-1);NTT(t[1],-1);for(ll i=0;i<k;i++){(g[0][i+1]+=t[0][i])%=P;(g[0][i+2]+=2ll*t[0][i])%=P;(g[0][i+3]+=t[0][i])%=P;(g[1][i+1]+=t[1][i])%=P;(g[1][i+2]+=2ll*t[1][i])%=P;(g[1][i+3]+=t[1][i])%=P;}for(ll i=k;i<n;i++)g[0][i]=g[1][i]=0;}if(d*2>m)break;for(ll j=0;j<3;j++){for(ll i=0;i<n;i++)t[j][i]=(i<k)?f[j][i]:0;NTT(t[j],1);}for(ll i=0;i<n;i++){f[0][i]=1ll*t[0][i]*t[0][i]%P;f[1][i]=1ll*t[1][i]*t[0][i]%P;f[2][i]=1ll*t[1][i]*t[1][i]%P;t[0][i]=1ll*t[1][i]*t[1][i]%P;t[1][i]=1ll*t[1][i]*t[2][i]%P;t[2][i]=1ll*t[2][i]*t[2][i]%P;}for(ll j=0;j<3;j++)NTT(f[j],-1),NTT(t[j],-1);for(ll i=0;i<k-1;i++){(f[0][i+1]+=t[0][i])%=P;(f[0][i+2]+=2ll*t[0][i])%=P;(f[0][i+3]+=t[0][i])%=P;(f[1][i+1]+=t[1][i])%=P;(f[1][i+2]+=2ll*t[1][i])%=P,(f[1][i+3]+=t[1][i])%=P;(f[2][i+1]+=t[2][i])%=P;(f[2][i+2]+=2ll*t[2][i])%=P;(f[2][i+3]+=t[2][i])%=P;}for(ll i=k;i<n;i++)f[0][i]=f[1][i]=f[2][i]=0;}for(ll i=1;i<k;i++)print(g[0][i]),putchar(' ');putchar('\n');}return 0;
}
1004 Counting Stars
题目大意
nnn个数字要求支持操作
- 将一个区间所有数字的最低位减去
- 将一个区间非零数字的最高位提高一位
每次操作玩输出操作区间的和
1≤T≤5,1≤n,q≤105,1≤ai≤1091\leq T\leq 5,1\leq n,q\leq 10^5,1\leq a_i\leq 10^91≤T≤5,1≤n,q≤105,1≤ai≤109
解题思路
每个数字操作一最多执行logai\log a_ilogai次所以这个操作暴力做就好了
然后第二个操作我们把每个数字的最高位分离出来就变成了区间乘二。
代码是lemondinosaur\text{lemondinosaur}lemondinosaur爷写的
1006 GCD Game
题目大意
nnn个数字每次操作的人可以选择一个数字xxx然后找一个1≤k<x1\leq k<x1≤k<x将xxx变为gcd(x,k)gcd(x,k)gcd(x,k)。
不能操作者输,求是否先手必胜
1≤T≤100,∑n≤106,1≤ai≤1071\leq T\leq 100,\sum n\leq 10^6,1\leq a_i\leq 10^71≤T≤100,∑n≤106,1≤ai≤107
解题思路
其实就是去掉任意个质因子,所以用线性筛筛出每个数字的最小质因子,然后递推出每个数字的质因子数就好了。
时间复杂度O(ai+∑n)O(a_i+\sum n)O(ai+∑n)
code
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1e7+10;
int T,n,cnt,pri[N/10],v[N],c[N];
void Prime(){v[1]=1;for(int i=2;i<N;i++){if(!v[i])pri[++cnt]=i,v[i]=i;for(int j=1;j<=cnt&&i*pri[j]<N;j++){v[i*pri[j]]=pri[j];if(i%pri[j]==0)break;}}return;
}
int main()
{Prime();for(int i=2;i<N;i++)c[i]=c[i/v[i]]+1;scanf("%d",&T);while(T--){scanf("%d",&n);int ans=0;for(int i=1,x;i<=n;i++)scanf("%d",&x),ans=ans^c[x];if(ans)puts("Alice");else puts("Bob");}return 0;
}
1009 Singing Superstar
题目大意
给一个字符串SSS,nnn次询问字符串aia_iai在SSS中出现的最大不重次数。
1≤T≤5,1≤n≤105,1≤∣ai∣≤30,∑∣S∣,∑∣ai∣≤4×1051\leq T\leq 5,1\leq n\leq 10^5,1\leq |a_i|\leq 30,\sum |S|,\sum|a_i|\leq 4\times 10^51≤T≤5,1≤n≤105,1≤∣ai∣≤30,∑∣S∣,∑∣ai∣≤4×105
解题思路
根据贪心的思想是能够匹配尽量匹配
把所有SSS长度不超过303030的子串建一棵TrieTrieTrie,然后记下每个子串的最后出现位置然后统计ansansans就好了。
code
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ull unsigned long long
using namespace std;
const int N=1e5+10;
int n,T,l,cnt,ans[N*30],t[N*30][26],last[N*30];
char s[N],st[N];
int main()
{scanf("%d",&T);while(T--){scanf("%s",s+1);l=strlen(s+1);cnt=1;last[1]=ans[1]=0;memset(t[1],0,sizeof(t[1]));for(int i=1;i<=l;i++){int x=1;for(int j=1;j<=30;j++){if(i-j+1<1)break;int c=s[i-j+1]-'a';if(!t[x][c]){cnt++;memset(t[cnt],0,sizeof(t[cnt]));t[x][c]=cnt;last[cnt]=ans[cnt]=0;}x=t[x][c];if(i-last[x]>=j)last[x]=i,ans[x]++;}}scanf("%d",&n);while(n--){scanf("%s",st+1);int tl=strlen(st+1),x=1;for(int i=tl;i>=1;i--)x=t[x][st[i]-'a'];printf("%d\n",ans[x]);}}return 0;
}