CCPC-Wannafly Winter Camp Day8 (Div2, onsite)

目录

 

A.Aqours

B.玖凛两开花

D.吉良吉影的奇妙计划

E.Souls-like Game

G.穗乃果的考试

J.去音乐会


A.Aqours

记一个数组dp[u]表示点u到叶子节点距离的最小值。按照顺序对每个叶子节点暴力往上搜,因为题目给出的树保证了深的点一定比浅的大,也就保证了距离近的点一定会先更新。搜的和回溯的时候顺便更新dp数组就行了。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N=5e6+7;
int p[N],dp[N],ans,n;
bool is[N];
int dfs(int u,int stp){
    if(u==0){ans=-1;return N;}
    if(~dp[u]){ans=dp[u]+stp;return dp[u];}
    int tmp=dfs(p[u],stp+1);
    dp[u]=min(stp,tmp+1);
    return dp[u];
}
int main()
{
    scanf("%d",&n);
    is[1]=true;p[1]=0;dp[1]=-1;
    for(int i=2;i<=n;i++){
        is[i]=true;dp[i]=-1;
        scanf("%d",&p[i]);
        is[p[i]]=false;
    }
    for(int i=1;i<=n;i++){
        if(is[i]){
            printf("%d ",i);
            dfs(i,0);
            printf("%d\n",ans);
        }
    }
}

B.玖凛两开花

看错好久的题。

假设答案是x,也就是说0~x-1的点都一定都要与x~n-1的点相连,注意x~n-1的点可以不与任何点相连,但0~x-1的点必须都相连。这个可以做二分匹配来判断。所以只要二分答案x就行了。

也可以不二分,直接从小到大枚举答案,如果当前点被连接了,就让连接的点重连,当前点也要连。这样就少了一个log的复杂度。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e4+7;
struct Edge{
    int v,nxt;
}e[N*2];
int p[N],edn;
void add(int u,int v){
    e[++edn]=(Edge){v,p[u]};p[u]=edn;
}
int mt[N];
bool vis[N];
bool dfs(int u,int md){
    for(int i=p[u];~i;i=e[i].nxt){
        int v=e[i].v;
        if(vis[v]||v<md) continue;
        vis[v]=true;
        if(mt[v]==-1||dfs(mt[v],md)){
            mt[v]=u;
            return true;
        }
    }
    return false;
}
int main()
{
    int n,m,u,v;
    scanf("%d%d",&n,&m);
    memset(p,-1,sizeof(p));edn=-1;
    for(int i=1;i<=m;i++){
        scanf("%d%d",&u,&v);
        if(u>v) swap(u,v);
        add(u,v);
    }
    memset(mt,-1,sizeof(mt));
    int ans=0;
    for(int i=0;i<n;i++){
        if(mt[i]!=-1){
            memset(vis,false,sizeof(vis));
            if(!dfs(mt[i],i+1)) break;
        }
        memset(vis,false,sizeof(vis));
        if(!dfs(i,i+1)) break;
        ans=i+1;
    }
    printf("%d\n",ans);
}

C.御坂妹妹

三分可以卡过,正解求上下凸包,留坑。


D.吉良吉影的奇妙计划

暴力打表或者dp

dp[i][j][p][q]表示第i位,有了j个左,第i-1位是p,第i位是q的答案。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=998244353;
ll dp[50][50][5][5];
int main()
{
    int n;
    scanf("%d",&n);
    dp[2][2][1][1]=1;
    dp[2][0][2][2]=1;
    dp[2][1][0][0]=1;
    for(int i=3;i<=n*2;i++){
        int hi=min(i,n);
        int lo=max(0,i-n);
        for(int j=lo;j<=hi;j++){
            if(j>=1){
                dp[i][j][0][0]=dp[i-1][j-1][2][2];
                dp[i][j][0][1]=dp[i-1][j-1][0][0];
                dp[i][j][1][1]=(dp[i-1][j-1][1][1]+dp[i-1][j-1][0][1])%mod;
            }
            if(i-1>=j){
                dp[i][j][0][0]=(dp[i][j][0][0]+dp[i-1][j][1][1])%mod;
                dp[i][j][0][2]=dp[i-1][j][0][0];
                dp[i][j][2][2]=(dp[i-1][j][2][2]+dp[i-1][j][0][2])%mod;
            }
        }
    }
    int ans=0;
    for(int i=0;i<3;i++){
        for(int j=0;j<3;j++){
            ans=(ans+dp[n*2][n][i][j])%mod;
        }
    }
    printf("%d\n",ans);
}

E.Souls-like Game

矩阵求区间修改和区间积。

div2好像纯暴力也可以,div1线段树加预处理矩阵的幂次。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N=1e4+7;
const ll mod=998244353;
struct rec{
    ll p[3][3];
    ll val;
    void read(){
        for(int i=0;i<3;i++) for(int j=0;j<3;j++) scanf("%lld",&p[i][j]);
    }
    ll getv(){
        val=0; for(int i=0;i<3;i++) for(int j=0;j<3;j++) val=(val+p[i][j])%mod;
        return val;
    }
    void operator = (const rec k){
        for(int i=0;i<3;i++) for(int j=0;j<3;j++) p[i][j]=k.p[i][j];
    }
    rec operator * (const rec k)const{
        rec res;
        for(int i=0;i<3;i++){
            for(int j=0;j<3;j++){
                res.p[i][j]=0;
                for(int t=0;t<3;t++){
                    res.p[i][j]=(res.p[i][j]+p[i][t]*k.p[t][j]%mod)%mod;
                }
            }
        }
        return res;
    }
    rec operator ^ (int k)const{
        rec res,a;
        for(int i=0;i<3;i++)
            for(int j=0;j<3;j++)
                res.p[i][j]=(i==j),a.p[i][j]=p[i][j];
        while(k){
            if(k&1) res=res*a;
            a=a*a;k>>=1;
        }
        return res;
    }
}t[N<<2],lz[N<<2];
void pu(int rt){
    t[rt]=t[rt<<1]*t[rt<<1|1];
}
void pd(int rt,int l,int r){
    if(lz[rt].val){
        int m=l+r>>1;
        t[rt<<1]=lz[rt]^(m-l+1);
        t[rt<<1|1]=lz[rt]^(r-m);
        lz[rt<<1].val=lz[rt<<1|1].val=1;
        lz[rt<<1]=lz[rt];lz[rt<<1|1]=lz[rt];
        lz[rt].val=0;
    }
}
void bd(int rt,int l,int r){
    lz[rt].val=0;
    if(l==r){
        t[rt].read();
        return;
    }
    int m=l+r>>1;
    bd(rt<<1,l,m);
    bd(rt<<1|1,m+1,r);
    pu(rt);
}
void upd(int rt,int l,int r,int L,int R,rec val){
    if(L<=l&&r<=R){
        t[rt]=val^(r-l+1);
        lz[rt]=val;
        lz[rt].val=1;
        return;
    }
    pd(rt,l,r);
    int m=l+r>>1;
    if(L<=m) upd(rt<<1,l,m,L,R,val);
    if(m<R) upd(rt<<1|1,m+1,r,L,R,val);
    pu(rt);
}
rec que(int rt,int l,int r,int L,int R){
    if(L<=l&&r<=R) return t[rt];
    pd(rt,l,r);
    int m=l+r>>1;
    rec res;
    for(int i=0;i<3;i++) for(int j=0;j<3;j++) res.p[i][j]=(i==j);
    if(L<=m) res=res*que(rt<<1,l,m,L,R);
    if(m<R) res=res*que(rt<<1|1,m+1,r,L,R);
    return res;
}
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);n--;
    bd(1,1,n);
    while(m--){
        int t,l,r;
        rec tmp;
        scanf("%d%d%d",&t,&l,&r);
        if(t==1){
            tmp.read();
            upd(1,1,n,l,r,tmp);
        }
        else{
            printf("%lld\n",que(1,1,n,l,r-1).getv());
        }
    }
}

G.穗乃果的考试

统计每个1对答案的贡献,队友代码

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const ll mod=998244353;

char s[2005];
int main()
{
    ll n,m;
    scanf("%lld%lld",&n,&m);
    ll ans=0;
    for(ll i=1;i<=n;i++)
    {
        scanf("%s",s+1);
        for(ll j=1;j<=m;j++)
            if(s[j]=='1')
                ans=(ans+(i*j)%mod*(n-i+1)%mod*(m-j+1)%mod)%mod;
    }
    printf("%lld\n",ans%mod);
    return 0;
}

J.去音乐会

两个范围其实都只有1e3,所以可以n^2暴力,然后就是一个拓展欧几里得解不定方程。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
void egcd(int a,int b,int &x,int &y){
    if(!b){x=1;y=0;}
    else{egcd(b,a%b,y,x);y-=x*(a/b);}
}
int main()
{
    int t,a,b,c,d,x,y;
    ll n,ans;
    scanf("%d",&t);
    while(t--){
        scanf("%d%d%d%d%lld",&a,&b,&c,&d,&n);
        ans=0;
        for(int i=a+1;i<=b;i++){
            for(int j=c+1;j<=d;j++){
                int aa=-b,bb=d,cc=i-j;
                int t=__gcd(aa,bb);
                if(cc%t) continue;
                aa/=t;bb/=t;cc/=t;
                egcd(aa,bb,x,y);
                if(bb<0) bb=-bb;
                x=((1ll*cc*x)%bb+bb)%bb;
                if(1ll*b*x+i<=n) ans+=((n-i)/b-x)/(bb)+1;
            }
        }
        printf("%lld\n",ans);
    }
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值