2023湖南省赛游记/题解

省赛拖了大哥们的后腿,感觉随便补个正常一队水平的人,我们一队肯定能AK。只能说自己真的菜,全程帮不上什么忙,还负贡献,真的想笑

B
暴力sg

#include 
#define ll long long
#define ull unsigned long long
constexpr int N=1e6;
int sg[N+5];
int vis[N+5];
void getsg(){
    for (int i=1;i<=N;i++){
        // vis[sg[i-1]]=i;
        int res=sqrt(1.0*i);
        if (i==res*res){
            for (int k=1;k<=res;k++){
                vis[sg[i-(i-res*res+k*res)]]=i;
            }
        }else{
            for (int k=0;k<=res;k++){
                vis[sg[i-(i-res*res+k*res)]]=i;
            }
        }
        while (vis[sg[i]]==i){
            sg[i]++;
        }
    }
}
void yrzr(){
    getsg();
    int n;
    std::cin>>n;
    int sum=0;
    for (int i=1;i<=n;i++){
        int x;
        std::cin>>x;
        sum^=sg[x];
    }
    if (sum){
        std::cout<<"First\n";
    }else{
        std::cout<<"Second\n";
    }
}

D
大哥赛后说是HN省选dp原题(见状压dp P3188 [HNOI2007] 梦幻岛宝珠),然后这个题比那个题范围还要大,看完题解发现是贪心题,但是还是按位分组的思想去解决问题
首先对于一个目标物品来说,它的需求可以表示为几个二进制位的组合,然后我们分解每一个目标物品到位上,然后将其按位合并,最后可以表示为一个二进制串(用数组表示)
然后对于每一件物品,重量小的物品可以合并成大物品,但大物品不能拆解成小物品,所以我们按位从小到大枚举目标物品数组,贪心地用小物品去满足当前位的数量,然后将剩下的物品贪心地两两合并成下一位的物品。

简要的说就是:目标物品总体可以表示为一个数组,然后拥有的物品贪心从小到大去填补目标重量

#include 
#define ll long long
#define ull unsigned long long
constexpr int inf=1e9;
std::priority_queue<int,std::vector<int>,std::greater<int>> q[5105];
void yrzr(){
    int n;
    std::cin>>n;
    for (int i=1;i<=n;i++){
        int x,y;
        std::cin>>x>>y;
        q[x].push(y);
    }
    
    int m;
    std::cin>>m;
    std::vector<int> sum(5105);
    while (m--){
        int t,h;
        std::cin>>t>>h;
        while (h){
            if (h&1) sum[t]++;
            h>>=1;
            t++;
        }
    }
    for (int i=0;i<=5100;i++){
        if (sum[i]>=2){
            int temp=sum[i],p=i;
            sum[i]=0;
            while (temp){
                if (temp&1) sum[p]++;
                temp>>=1;
                p++;
            }
        }
    }

    int ans=0;
    for (int i=0;i<=5100;i++){
        if (q[i].size()<sum[i]){
            std::cout<<"-1\n";
            return;
        }

        for (int j=1;j<=sum[i];j++){
            ans+=q[i].top();
            q[i].pop();
        }
        while (q[i].size()>=2){
            int temp=q[i].top();
            q[i].pop();
            temp+=q[i].top();
            q[i].pop();
            q[i+1].push(temp);
        }
    }
    std::cout<<ans;
}

F
图论题,首先Floyed求出一个字母变成另一个字母的最小花费,然后对于两个不同的字母,暴力预处理出他们变成同样字母的最小花费。最后A串不动,枚举B串位置,求最小即可

可能有重边,输入时 d i s dis dis取最小

#include 
#define int long long
#define ull unsigned long long
constexpr int inf=1e16;
void yrzr(){
    int n,m;
    std::cin>>n>>m;
    std::vector<int> a(n),b(n);
    for (int i=0;i<n;i++){
        std::cin>>a[i];
    }
    for (int i=0;i<n;i++){
        std::cin>>b[i];
    }

    std::vector<std::vector<int>> dis(405,std::vector<int>(405,inf)),c(405,std::vector<int>(405,inf));
    for (int i=1;i<=m;i++){
        int x,y,z;
        std::cin>>x>>y>>z;
        dis[x][y]=std::min(dis[x][y],z);
    }

    for (int i=1;i<=400;i++){
        dis[i][i]=c[i][i]=0;
    }

    for (int k=1;k<=400;k++){
        for (int i=1;i<=400;i++){
            for (int j=1;j<=400;j++){
                dis[i][j]=std::min(dis[i][j],dis[i][k]+dis[k][j]);
            }
        }
    }


    for (int i=1;i<=400;i++){
        for (int j=1;j<=400;j++){
            for (int k=1;k<=400;k++){
                c[i][j]=std::min(c[i][j],dis[i][k]+dis[j][k]);
            }
        }
    }

    int ans=inf;
    for (int i=0;i<n;i++){
        int sum=0;
        bool ok=1;
        for (int j=0;j<n;j++){
            if (c[a[j]][b[(j+i)%n]]==inf){
                ok=0;
                break;
            }else{
                sum+=c[a[j]][b[(j+i)%n]];
            }
        }
        if (ok){
            ans=std::min(ans,sum);
        }
    }
    if (ans==inf){
        std::cout<<"-1\n";
    }else{
        std::cout<<ans<<"\n";
    }
}

I
数位dp,不知道省赛脑子犯了什么抽,受前几天一道数位dp,直接把状态数组压入map暴力记忆化的影响,想当然的觉得这题也是这么做,大哥直接提出时间复杂度是不够的。后面大哥们说是数位dp加个容斥。

补题的时候也一直在想数位dp+容斥,发现自己容斥不出来(容斥见识的得太少了),看了题解后发现只要一个数位dp就可以了

突然发现这题的思想其实跟前几天做的数位dp是一样的,甚至状态更简单,我们最重要的状态其实就只有当前位第几位,然后还要知道现在有多少个不同的数,这就足够了。因为对于这10个数字来说,后面的数字它只有两种身份,之前选了/没选,所以当前位相同,之前不同的数字个数相同,剩下的方案数就相同,就可以用记忆化。(赛时真的蠢傻了,这种套路思想学数位dp的时候见过很多了)

#include 
#define ll long long
#define ull unsigned long long
constexpr ll mod=1e9+7;
ll f[200005][12][2][2];
int now[12];
void yrzr(){
    int n;
    std::cin>>n;
    std::vector<int> a(n+1),b(n+1);
    for (int i=1;i<=n;i++){
        char ch;
        std::cin>>ch;
        a[i]=ch-'0';
    }
    for (int i=1;i<=n;i++){
        char ch;
        std::cin>>ch;
        b[i]=ch-'0';
    }
    reverse(a.begin()+1,a.end());
    reverse(b.begin()+1,b.end());

    int A;
    std::cin>>A;
    std::function<ll(int,int,int,int)> dfs=[&](int len,int sum,int dif1,int dif2){
        if (len==0){
            return 1LL*(sum==A);
        }
        if (sum>A){
            return 0LL;
        }
        if (dif1&&dif2&&f[len][sum][dif1][dif2]){
            return f[len][sum][dif1][dif2];
        }

        int x=(dif1?0:a[len]),y=(dif2?9:b[len]);
        ll res=0;
        for (int i=x;i<=y;i++){
            if (now[i]){
                now[i]++;
                res=(res+dfs(len-1,sum,dif1|(a[len]!=i),dif2|(b[len]!=i)))%mod;
                now[i]--;
            }else{
                now[i]++;
                res=(res+dfs(len-1,sum+1,dif1|(a[len]!=i),dif2|(b[len]!=i)))%mod;
                now[i]--;
            }
        }
        if (dif1&&dif2) f[len][sum][dif1][dif2]=res;
        return res;
    };
    std::cout<<dfs(n,0,0,0);

}

J
分别枚举圆心在三个坐标轴,然后二分圆的半径,接下来考虑怎么check,对于一个点来说,我们去计算出它在球内时,圆心所在的区间,存下区间左右端点和权值(表示是左端点还是右端点),遍历一遍看是否有点满足题意即可

#include 
#define ll long long
#define ull unsigned long long
constexpr int N=1e5+5;
int n;
ll x[N],y[N],z[N];
struct node{
    double pos;
    int val;
}a[N<<1];
bool check(double r){
    int cnt=0;
    for (int i=1;i<=n;i++){
        double res=r*r-y[i]*y[i]-z[i]*z[i];
        if (res<0) continue;
        res=sqrt(res);
        a[++cnt]={x[i]-res,1};
        a[++cnt]={x[i]+res,-1};
    }
    std::sort(a+1,a+1+cnt,[&](node i,node j){
        return i.pos<j.pos;
    });
    int sum=0;
    for (int i=1;i<=cnt;i++){
        sum+=a[i].val;
        if (sum>=n/2){
            return 1;
        }
    }

    cnt=0;
    for (int i=1;i<=n;i++){
        double res=r*r-x[i]*x[i]-z[i]*z[i];
        if (res<0) continue;
        res=sqrt(res);
        a[++cnt]={y[i]-res,1};
        a[++cnt]={y[i]+res,-1};
    }
    std::sort(a+1,a+1+cnt,[&](node i,node j){
        return i.pos<j.pos;
    });
    sum=0;
    for (int i=1;i<=cnt;i++){
        sum+=a[i].val;
        if (sum>=n/2){
            return 1;
        }
    }

    cnt=0;
    for (int i=1;i<=n;i++){
        double res=r*r-x[i]*x[i]-y[i]*y[i];
        if (res<0) continue;
        res=sqrt(res);
        a[++cnt]={z[i]-res,1};
        a[++cnt]={z[i]+res,-1};
    }
    std::sort(a+1,a+1+cnt,[&](node i,node j){
        return i.pos<j.pos;
    });
    sum=0;
    for (int i=1;i<=cnt;i++){
        sum+=a[i].val;
        if (sum>=n/2){
            return 1;
        }
    }

    return 0;
}
void yrzr(){
    std::cin>>n;
    for (int i=1;i<=n;i++){
        std::cin>>x[i]>>y[i]>>z[i];
    }

    double l=0,r=1e6;
    for (int i=0;i<=50;i++){
        double mid=(l+r)/2;
        if (check(mid)){
            r=mid;
        }else{
            l=mid;
        }
    }
    std::cout<<std::fixed<<std::setprecision(7)<<l;
}

K
补的时候相成状压了,然后状压+矩阵加速???,发现不会写
然后看了题解知道是容斥了,所以以后对于这种小数据除了状压也有可能是容斥

想到容斥就很好做了,每次直接暴力二进制枚举哪些城市不去,在转移的矩阵和答案矩阵上去掉那些转移点,然后dp矩阵加速即可

#include 
#define ll long long
#define ull unsigned long long
constexpr ll mod=1e9+9;
constexpr int N=25;
struct Matrix{
    int lenn,lenm;
    ll a[N][N];
    Matrix() {memset(a,0,sizeof(a));}
 
    Matrix operator*(const Matrix &b)const{
        if (lenm==b.lenn){
            Matrix res;
            res.lenn=lenn;
            res.lenm=b.lenm;
            for (int i=1;i<=lenn;i++){
                for (int j=1;j<=b.lenm;j++){
                    for (int l=1;l<=lenm;l++){
                        res.a[i][j]=(res.a[i][j]+a[i][l]*b.a[l][j]%mod)%mod;
                    }
                }
            }
            return res;
        }
    }

    Matrix operator+(const Matrix &b)const{
        if (lenn==b.lenn&&lenm==b.lenm){
            Matrix res;
            res.lenn=lenn;
            res.lenm=lenm;
            for (int i=1;i<=lenn;i++){
                for (int j=1;j<=lenm;j++){
                    res.a[i][j]=a[i][j]+b.a[i][j];
                }
            }
            return res;
        }
    }

    Matrix operator-(const Matrix &b)const{
        if (lenn==b.lenn&&lenm==b.lenm){
            Matrix res;
            res.lenn=lenn;
            res.lenm=lenm;
            for (int i=1;i<=lenn;i++){
                for (int j=1;j<=lenm;j++){
                    res.a[i][j]=a[i][j]-b.a[i][j];
                }
            }
            return res;
        }
    }

};
int s[10],n,m,k,d;
int a[25][25];
ll solve(int mask){
    Matrix e,f,ans;
    ans.lenn=1;
    ans.lenm=f.lenn=f.lenm=e.lenn=e.lenm=n;
    std::vector<int> vis(n+1);
    
    for (int i=1;i<=n;i++){
        for (int j=1;j<=n;j++){
            e.a[i][j]=a[i][j];
        }
        f.a[i][i]=1;
        ans.a[1][i]=1;
    }
    
    for (int i=1;i<=k;i++){
        if (mask&(1<<(i-1))){
            int x=s[i];
            for (int j=1;j<=n;j++){
                e.a[x][j]=e.a[j][x]=0;
            }
            vis[x]=1;
            ans.a[1][x]=0;
        }
    }

    int temp=d;
    while (temp){
        if (temp&1) f=f*e;
        e=e*e;
        temp>>=1;
    }
    ans=ans*f;
    ll sum=0;
    for (int i=1;i<=n;i++){
        if (vis[i]) continue;
        sum=(sum+ans.a[1][i])%mod;
    }
    return sum;
}
void yrzr(){
    std::cin>>n>>m>>k>>d;
    d--;
    for (int i=1;i<=k;i++){
        std::cin>>s[i];
    }
    for (int i=1;i<=m;i++){
        int x,y;
        std::cin>>x>>y;
        a[x][y]=a[y][x]=1;
    }
    ll ans=0;
    for (int i=0;i<(1<<k);i++){
        if (__builtin_popcount(i)&1){
            ans=(ans-solve(i)+mod)%mod;
        }else{
            ans=(ans+solve(i))%mod;
        }
    }
    std::cout<<ans;
}

你可能感兴趣的:(算法)