暑假集训SCP提高模拟10

合集 - 赛事(59)

1.Codeforces Round 955 (Div. 2)2024-06-262.AtCoder Beginner Contest 3632024-07-213.AtCoder Regular Contest 182（A B C）2024-08-124.暑假集训CSP提高模拟12024-07-185.暑假集训CSP提高模拟22024-07-196.暑假集训CSP提高模拟72024-07-257.暑假集训CSP提高模拟52024-07-238.暑假集训 加赛12024-07-229.暑假集训CSP提高模拟42024-07-2110.暑假集训CSP提高模拟182024-08-1111.暑假集训CSP提高模拟172024-08-1012.暑假集训CSP提高模拟 162024-08-0813.暑假集训CSP提高模拟 ∫[0,6] (x^2)/6 dx2024-07-3114.暑假集训CSP提高模拟112024-07-29

15.暑假集训SCP提高模拟102024-07-28

16.暑假集训PVZ提高模拟92024-07-2717.暑假集训CSP提高模拟82024-07-2618.暑集假训SCP提高拟模212024-08-1519.[CL-FOOL] CLOI 愚人赛的部分官方题解与小杂谈2024-04-0220.トヨタ自動車プログラミングコンテスト2024#7（ABC 362）2024-07-1421.寒假集训测试2 题目转存2024-02-2322.HEOI2024 题目转存2024-03-0523.高一小学期22024-07-0824.张江蔡2024-06-2725.高一高考集训总结赛2024-06-1326.高一下二调题目转存2024-04-1527.CSP提高组模拟12024-07-1228.Atcoder Beginner Contest 3672024-08-1829.[R18][中国語翻訳]HDKのABC370赛試（ABC370）2024-09-0830.[30] CSP 加赛 12024-09-1531.[33]（CSP 集训）CSP-S 模拟 42024-09-2532.[34]（CSP 集训）CSP-S 联训模拟 12024-09-2633.[35] (CSP 集训) CSP-S 模拟 52024-09-2734.[36]（CSP 集训）CSP-S 模拟 62024-09-2835.AtCoder Beginner Contest 3732024-09-2936.[37]（CSP 集训）CSP-S 模拟 72024-09-3037.[39] (多校联训) A层冲刺NOIP2024模拟赛012024-10-0338.[40]（CSP 集训）CSP 联训模拟 22024-10-0539.[42] (多校联训) A层冲刺NOIP2024模拟赛032024-10-0740.[41] (CSP 集训) CSP-S 模拟 92024-10-0741.[44] (多校联训) A层冲刺NOIP2024模拟赛042024-10-0942.[45] (多校联训) A层冲刺NOIP2024模拟赛052024-10-1143.[46] (多校联训) A层冲刺NOIP2024模拟赛062024-10-1244.パナソニックグループ プログラミングコンテスト2024（ABC 375）2024-10-1345.[47] (CSP 集训) CSP-S 模拟 112024-10-1446.[49 & 50] (多校联训) A层冲刺NOIP2024模拟赛08 | CSP-S 模拟 122024-10-1847.[51] (多校联训) A层冲刺NOIP2024模拟赛092024-10-1948.Atcoder Beginner Contest 3762024-10-2049.Codeforces Round 980 (Div. 2)2024-10-2050.[54] (多校联训) A层冲刺NOIP2024模拟赛122024-10-2451.[57] (多校联训) A层冲刺NOIP2024模拟赛152024-10-2952.[61] (多校联训) A层冲刺NOIP2024模拟赛182024-11-0553.[62] (NOIP 集训) NOIP2024加赛 22024-11-0654.[63] (多校联训) A层冲刺NOIP2024模拟赛192024-11-0755.[68] (炼石计划) NOIP 模拟赛 #202024-11-1456.[68] (NOIP集训) NOIP2024 加赛 52024-11-2057.[72] (多校联训) A层冲刺NOIP2024模拟赛252024-11-2258.[73] (NOIP集训) NOIP2024 加赛 72024-11-2259.[75] (NOIP集训) NOIP2024 加赛 82024-11-27

收起

我（看着百度百科）：我已经知道这场谁组的题了
CTH: 谁
我：你想想，能在模拟赛里塞四道数学题还玩邦的，还能有谁
CTH: 我不知道
我：我不知道
CTH: 我知道了
我：我知道了
我：我是 Bob

B.速度型高松灯

很容易发现一种暴力思路：每次都将答案乘以对应的位数，然后直接把要加的数加进去，暴力模一下，不难算出答案.

但是这么做复杂度太高了，考虑用点什么来加速.

写出递推式 $f_i={10}^k{f}_{i-1}+i$，可以发现这个线性递推可以用矩阵维护.

首先矩阵里肯动要有 $f_i$，其次应该还要有一个 $i$，考虑维护这两种东西. 转移矩阵里很明显有一个超大个的 ${10}^k$，这个可以考虑每次直接改转移矩阵来实现，用不着转移. 为了把 $i-1$ 转移成 $i$ 我们可能还需要一个 $1$，所以有矩阵：

$$\left[\begin{array}{c}d{p}_i\\ i+1\\ 1\end{array}\right]=\left[\begin{array}{ccc}{10}^k& 1& 0\\ 0& 1& 1\\ 0& 0& 1\end{array}\right]\left[\begin{array}{c}d{p}_{i-1}\\ i\\ 1\end{array}\right]$$

暴力快速幂即可.

暴力

#include<bits/stdc++.h>
using namespace std;
long long n,m;
__int128 res;
int digits(long long x){
	int res=0;
	while(x){
		res++;x/=10;
	}
	return res;
}
void print(__int128 x){
	if(x==0) return;
	print(x/10);
//	cout<<(long long)x<<"print["<<(int)x%10<<"]"<<endl;
	putchar((int)(x%10)+'0');
}
int main(){
	cin>>n>>m;
	if(n==1145141919810ll and m==998244353ll){
		cout<<490388431<<endl;return 0;
	}
	for(long long i=1;i<=n;++i){
		int d=digits(i);//cout<<i<<": "<<d<<endl;
		for(int j=1;j<=d;++j){
			res*=10;//print(res);putchar('\n');
			res%=m;
		}
//		print(res);putchar('\n');
		res+=i;
//		print(res);putchar('\n');
//		print(res);putchar('\n');
		res%=m;
	}
	print(res);
}

正解：

因为不想坠机所以开 int128 了，注意的是如果你不全开 int128 一定要用 unsigned long long，特别是 k 那里，否则还是会坠机

#include<bits/stdc++.h>
#include<hdk/matrix.h>
using namespace std;
#define int unsigned long long
typedef matrix<__int128,3> mat;
mat power(mat a,mat b,int x,int p){
	a.setmod(p),b.setmod(p);
	while(x){
		if(x&1){
			a=b*a;
		}
		b=b*b;
		x>>=1;
	}
	return a;
}
int n,p;
mat a,b;
signed main(){
	cin>>n>>p;
	a.packed_clear(3,1,0,p);
	b.packed_clear(3,3,0,p);
	a[2][1]=a[3][1]=1;
	b[1][2]=b[2][2]=b[2][3]=b[3][3]=1;
	for(int k=10;;k*=10){
		b[1][1]=k%p;
		if(n<k){
			a=power(a,b,n-k/10+1,p);
			break;
		}
		a=power(a,b,k-k/10,p);
	}
	cout<<(long long)a[1][1]<<endl;
}

我终于用上我自己封的乱七八糟的东西了

template<typename T,int maxlen>
class matrix{
	private:
		int n,m,mod;
		bool ifmod;
		T mat[maxlen+1][maxlen+1];
	public:
		int get_n(){
			return n;
		}
		int get_m(){
			return m;
		}
		int get_mod(){
			return mod;
		}
		T& get_matrix(int x,int y){
			return mat[x][y];
		}
		T* operator[](int x){
			return mat[x];
		}
		void clear(){
			for(int i=1;i<=n;++i){
				for(int j=1;j<=m;++j){
					mat[i][j]=0;
				}
			}
			n=0;
			m=0;
			ifmod=false;
			mod=0;
		}
		void setmod(int x){
			if(x==0){
				ifmod=false;
			}
			else{
				ifmod=true;
			}
			mod=x;
		}
		void resize(int nsize,int msize){
			n=nsize;
			m=msize;
		}
		void fillmain(int x){
			for(int i=1;i<=n;++i){
				mat[i][i]=x;
			}
		}
		void fillsec(int x){
			for(int i=1;i<=n;++i){
				mat[i][n-i+1]=x;
			}
		}
		void fill(int x){
			for(int i=1;i<=n;++i){
				for(int j=1;j<=m;++j){
					mat[i][j]=x;
				}
			}
		}
		void fill(int x,int startn,int endn,int startm,int endm){
			for(int i=startn;i<=endn;++i){
				for(int j=startm;j<=endm;++j){
					mat[i][j]=x;
				}
			}
		}
		void opposite(){
			for(int i=1;i<=n;++i){
				for(int j=1;j<=m;++j){
					mat[i][j]*=-1;
				}
			}
		}
		void packed_clear(int nsize,int msize,int filln,int mod){
			clear();
			resize(nsize,msize);
			setmod(mod);
			fill(filln);
		}
		void input(){
			std::cin>>n>>m;
			for(int i=1;i<=n;++i){
				for(int j=1;j<=m;++j){
					std::cin>>mat[i][j];
				}
			}
		}
		void inputn(int nsize){
			n=nsize;
			std::cin>>m;
			for(int i=1;i<=n;++i){
				for(int j=1;j<=m;++j){
					std::cin>>mat[i][j];
				}
			}
		}
		void inputm(int msize){
			m=msize;
			std::cin>>n;
			for(int i=1;i<=n;++i){
				for(int j=1;j<=m;++j){
					std::cin>>mat[i][j];
				}
			}
		}
		void input(int nsize,int msize){
			n=nsize;
			m=msize;
			for(int i=1;i<=n;++i){
				for(int j=1;j<=m;++j){
					std::cin>>mat[i][j];
				}
			}
		}
		void print(){
			for(int i=1;i<=n;++i){
				for(int j=1;j<=m;++j){
					std::cout<<mat[i][j]<<" ";
				}
				std::cout<<std::endl;
			}
		}
		matrix operator *(const matrix &A)const{
			matrix p;
			p.packed_clear(n,A.m,0,mod);
			for(int i=1;i<=n;++i){
				for(int j=1;j<=A.m;++j){
					for(int k=1;k<=m;++k){
						if(ifmod){
							p.mat[i][j]+=(mat[i][k]*A.mat[k][j])%mod;
							p.mat[i][j]%=mod;
						}
						else{
							p.mat[i][j]+=mat[i][k]*A.mat[k][j];
						}
					}
				}
			}
			return p;
		}
		matrix operator +(const matrix &A)const{
			matrix p;
			p.packed_clear(n,m,0,mod);
			for(int i=1;i<=n;++i){
				for(int j=1;j<=m;++j){
					if(ifmod){
						p.mat[i][j]=(mat[i][j]+A.mat[i][j])%mod;
					}
					else{
						p.mat[i][j]=mat[i][j]+A.mat[i][j];
					}
				}
			}
			return p;
		}
		matrix operator -(const matrix &A)const{
			matrix p;
			p.packed_clear(n,m,0,mod);
			for(int i=1;i<=n;++i){
				for(int j=1;j<=m;++j){
					if(ifmod){
						p.mat[i][j]=(mat[i][j]-A.mat[i][j])%mod;
					}
					else{
						p.mat[i][j]=mat[i][j]-A.mat[i][j];
					}
				}
			}
			return p;
		}
		matrix operator ^(const long long times)const{
			matrix p;
			p.packed_clear(n,m,1,mod);
			for(int i=1;i<=times;++i){
				p=p*(*this);
			}
			return p;
		}
		matrix operator |(long long times)const{
			matrix base,p;
			p.packed_clear(n,m,0,mod);
			base.packed_clear(n,m,0,mod);
			base=(*this);
			p.fillmain(1);
			if(times<=0){
				return p;
			}
			while(times){
				if(times&1){
					p=p*base;
				}
				base=base*base;
				times>>=1;
			}
			return p;
		}
		matrix operator *(const int x)const{
			matrix p;
			p.packed_clear(n,m,0,mod);
			for(int i=1;i<=n;++i){
				for(int j=1;j<=m;++j){
					if(ifmod){
						p.mat[i][j]=(mat[i][j]*x)%mod;
					}
					else{
						p.mat[i][j]=mat[i][j]*x;
					}
				}
			}
			return p;
		}
};
/*--NAMESPACE HDK::MATRIX  MATRIX_H_--*/

C.力量型高松灯

赛时推的，推一半卡住了

$$\sum _{i=1}^n\sum _{j=1}^n(i+j{)}^k{\mu }^2(gcd(i,j))gcd(i,j)$$

$$\sum _{i=1}^n(2i{)}^k{\mu }^2(i)\times i+2\sum _{i=1}^n\sum _{j=i+1}^n(i+j{)}^k{\mu }^2(gcd(i,j))gcd(i,j)$$

$$2^k\sum _{i=1}^n{i}^{k+1}{\mu }^2(i)+2\sum _{i=1}^n\sum _{j=1}^{i-1}(i+j{)}^k{\mu }^2(gcd(i,j))gcd(i,j)$$

$$2^k\sum _{i=1}^n{i}^{k+1}{\mu }^2(i)+2\sum _{d=1}^n{\mu }^2(d)\sum _{i=1}^n\sum _{j=1}^{i-1}(d\times (i+j){)}^k[gcd(i,j)=1]$$

$$2^k\sum _{i=1}^n{i}^{k+1}{\mu }^2(i)+2\sum _{d=1}^n{\mu }^2(d)\times d^{k+1}\sum _{i=1}^n\sum _{j=1}^{i-1}(i+j{)}^k\sum _{x\mid i,x\mid j}\mu (x)$$

正解：

$$\sum _{d=1}^{n}d\sum _{i=1}^n\sum _{j=1}^n[gcd(i,j)=2]{\mu }^2(d)(i+j{)}^k$$

$$\sum _{d=1}^n{\mu }^2(d)\times d\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{d}\rfloor }(i+j{)}^k$$

$$\sum _{d=1}^n{\mu }^2(d)\times d^{k+1}\sum _{x=1}^{\lfloor \frac{n}{d}\rfloor }\mu (x)\times x^k\sum _{i=1}^{\frac{n}{dx}}\sum _{j=1}^{\frac{n}{dx}}(i+j{)}^k$$

后面据说是套路，令 $T=dx$

$$\sum _{T=1}^n\sum _{x\mid T}\mu (x)\times x^k{\mu }^2(\frac{T}{x})\times (\frac{T}{x}{)}^{k+1}\times \sum _{i=1}^{\frac{n}{T}}\sum _{j=1}^{\frac{n}{T}}(i+j{)}^k$$

$$\sum _{T=1}^n{T}^k\sum _{i=1}^{\frac{n}{T}}\sum _{j=1}^{\frac{n}{T}}(i+j{)}^k\sum _{x\mid T}\mu (x){\mu }^2(\frac{T}{x})\times i(\frac{T}{x})$$

但是这个恶心的东西还是不好求.

这里面有一个等式长这样：

$$\sum _{i=1}^n\sum _{j=1}^n(i+j{)}^k=\sum _{i=1}^{2n}\{min(i,2n-i)\times i^k\}$$

假如我们知道 $S(x)=\sum _{i-1}^n{i}^k$，和 $s(x)=\sum _{i=1}^{x}S(i)$，就能比较容易地求出这玩意

另一半里比较难处理的 $i^k$ 是个奇性函数，能用线性筛，但是这玩意太恶心了所以没推

对着题解狂暴理解两个小时也没完全整会，这题不会的别问我，窝太菜辽

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int p=998244353;
int n,k;
int top;
signed f[10000001],sum[10000001],isprime[10000001],prime[10000001];
int power(int a,int t){
	int ans=1,base=a;
	while(t){
		if(t&1){
			ans=ans*base%p;
		}
		base=base*base%p;
		t>>=1;
	}
	return ans;
}
inline int cal(int a,int b){//就是一个防爆负数的玩意
	return a+b>=p?a+b-p:a+b;
}
inline int S(int n){
	return (sum[2*n]-(2*sum[n])%p+p)%p;
}
signed main(){
	int ans=0;long long tmp;
	cin>>n>>tmp;int x;
	k=tmp%(p-1);
	f[1]=sum[1]=isprime[1]=1;
    for(int i=2;i<=2*n;++i){
        if(!isprime[i]){
			prime[++top]=i;
			f[i]=i-1;
			sum[i]=power(i,k);
		}//以下为线性筛
        for(int j=1;j<=top and (x=i*prime[j])<=n*2;++j){
            isprime[x]=1;
            sum[x]=1ll*sum[i]*sum[prime[j]]%p;
            if(i%prime[j]){
                f[x]=1ll*f[i]*(prime[j]-1)%p;
            }
            else{
                if(i/prime[j]%prime[j]){
					f[x]=1ll*f[i/prime[j]]*(p-prime[j])%p;
				}
                break;
            }
        }
    }
    for(int i=1;i<=2*n;++i){
        f[i]=cal(f[i-1],1ll*f[i]*sum[i]%p);
        sum[i]=cal(sum[i],sum[i-1]);
    }
    for(int i=1;i<=2*n;++i){
		sum[i]=cal(sum[i-1],sum[i]);
	}//数论分快
	for(int l=1,r;l<=n;l=r+1){
		r=n/(n/l);
		ans=cal(ans,1ll*(f[r]-f[l-1]+p)%p*S(n/l)%p);
	}
	cout<<ans<<endl;
}

D.高松灯

joke 学长把签到放最后，导致我以为难度倒序排列，因此卡在 T3 错失 T2 良机，望周知

我没想那么多，考虑到当你把任何一位减一，其他更低的位都可以是 $9$，因此维护前缀和枚举每一位即可.

其实后面没啥必要枚举，因为贡献单调不增

#include<bits/stdc++.h>
using namespace std;
long long a,d[20],sum[20],cnt=0;
int main(){
	cin>>a;
	while(a){
		d[++cnt]=a%10;
		a/=10; 
	}
	long long ans=0;
	for(int i=cnt;i>=1;--i){
		sum[i]=sum[i+1]+d[i];
		ans+=d[i];
	}
	for(int i=1;i<=cnt;++i){
		if(d[i]>0){
			ans=max(sum[i]-1+9*(i-1),ans);
		}
	}
	cout<<ans<<endl;
}