NOIp2021题目解析与总结

最新推荐文章于 2023-07-07 13:56:58 发布

入坑信奥的L同学

最新推荐文章于 2023-07-07 13:56:58 发布

阅读量552

点赞数 1

分类专栏：练习题总结文章标签： NOIp真题筛法动态规划

本文链接：https://blog.csdn.net/m0_55341679/article/details/127436718

版权

练习题总结专栏收录该内容

2 篇文章 0 订阅

订阅专栏

T1.报数
这题没什么好说的，就是暴力筛，筛就完了。但是这题需要惊人的胆识（1s跑1e7的量也敢写），所以本蒟蒻只拿了50pts。代码如下：

#include<iostream>
#include<cstdio>
using namespace std;
const int N=1e7+1005;
bool st[N];
int ans[N];
bool check(int x){
    while(x){
        if(x%10==7) return true;
        x/=10;
    }
    return false;
}
void init(){
    for(int i=1;i<=10001000;i++){
        if(check(i)){
            for(int j=1;j*i<=10001000;j++){
                st[j*i]=true;
            }
        }
    }
    int back=0;
    for(int i=10001000;i>=1;i--){
        if(!st[i]){
            ans[i]=back;
            back=i;
        }
    }
}
int main(){
    //freopen("number.in","r",stdin);
    //freopen("number.out","w",stdout);
    init();
    int T;
    cin>>T;
    while(T--){
        int x;
        cin>>x;
        if(st[x]){
            cout<<-1<<endl;
        }
        else{
            cout<<ans[x]<<endl;
        }
    }
    return 0;
}

T2.数列
这题的难度就起来了，是一个计数类的DP。看到n，m，k的值都不大，考虑设计多维状态：用dp_i,j,k,p表示已经确定了S的前i位，且a数组中已经确定了j个数（这些数不一定连续），S中有k个1，当前位置（第i-1位，从0位算起的第i位）进到第i位的进位数值大小为p。这里考虑使用向后转移的方法，即由dp_i,j,k,p推到后面的状态：dp_{i+1,j+t,k+(t+p)%2,(t+p)/2}，这里t的意思是枚举a中有多少个数的值为i，那么别的状态也可以由此导出。那么dp_i,j,k,p对于后面状态的贡献就是dp_{i+1,j+t,k+(t+p)%2,(t+p)/2} * v_i^t * C^t_n-j，那么由此可以得出代码：

#include<iostream>
using namespace std;
const int N=35,M=105,mod=998244353;
typedef long long ll;
int n,m,K;
int v[M];
ll pv[M][N];
ll C[N][N];
ll f[M][N][N][N];
void init(){
    for(int i=0;i<=m;i++){
        pv[i][0]=1;
        for(int j=1;j<=n;j++){
            pv[i][j]=pv[i][j-1]*v[i]%mod;
        }
    }
    for(int i=0;i<=30;i++){
        for(int j=0;j<=i;j++){
            if(!j) C[i][j]=1;
            else C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
        }
    }
}
inline int popcnt(int x){
    int res=0;
    while(x){
        if(x&1) res++;
        x>>=1;
    }
    return res;
}
int main(){
    cin>>n>>m>>K;
    for(int i=0;i<=m;i++){
        cin>>v[i];
    }
    init();
    f[0][0][0][0]=1;
    for(int i=0;i<=m;i++){
        for(int j=0;j<=n;j++){
            for(int k=0;k<=K;k++){
                for(int p=0;p<=(n>>1);p++){
                    for(int t=0;t<=n-j;t++){
                        f[i+1][j+t][k+(t+p)%2][(t+p)/2]=(f[i+1][j+t][k+(t+p)%2][(t+p)/2]+f[i][j][k][p]*pv[i][t]%mod*C[n-j][t]%mod)%mod;
                    }
                }
            }
        }
    }
    ll res=0;
    for(int k=0;k<=K;k++){
        for(int p=0;p<=(n>>1);p++){
            if(k+popcnt(p)<=K) res=(res+f[m+1][n][k][p])%mod;
            //这里注意加上最后一位的进位
        }
    }
    cout<<res<<endl;
}