[41] (CSP 集训) CSP-S 模拟 9

A.邻面合并

观察到 \(m\) 很小,支持我们 \(mn2^{2m}\)

状压枚举二进制状态, \(f_{i,j,k}\) 表示到第 \(i\) 位的状态为 \(j\),上一位状态为 \(j\) (\(i,j\) 为状压位) 的方案数

考虑转移的时候加了什么

判断这一行与上一行的联通情况,如果这一行的某一个连通块和上一行正好对上了,那么就可以直接扩展矩形,答案不增加,否则对每个连通块,答案加一

我对这个联通情况的设计是不咋好的,还是说一下

\(a\) 是原矩阵第 \(i\) 行,\(b\) 是枚举到的状压状态,我判断连通块是根据

  • 连通块内不存在 \(a_j=0\)\(j\)
  • 连通块内不存在 \(b_j=0\)\(j\)

这么一做就会出很多冗余状态,但是我也不太会设计更好的了,所以就只能这样了

#include<bits/stdc++.h>
using namespace std;
int n,m;
int a[101][9];
int f[2][(1<<8)+1][(1<<8)+1];
int minn[101][(1<<8)+1];
vector<pair<int,int>>thi,las;
inline void work(int i,int j,vector<pair<int,int>>&ts){
    if(i==0) return;
    ts.clear();
    for(int k=1;k<=m;++k){
        int tmp=(j>>(k-1))&1;
        int tmplast=(k==1?tmp:(j>>(k-2))&1);
        if(a[i][k]==0){
            if(!ts.empty() and ts.back().second==0){
                ts.back().second=k-1;
            }
        }
        else if(k==1 or tmp!=tmplast or a[i][k-1]==0){
            if(!ts.empty() and ts.back().second==0){
                ts.back().second=k-1;
            }
            ts.push_back({k,0});
        }
    }
    if(!ts.empty() and ts.back().second==0){
        ts.back().second=m;
    }
}
int main(){
    freopen("merging.in","r",stdin);
    freopen("merging.out","w",stdout);
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;++i){
        for(int j=1;j<=m;++j){
            scanf("%d",&a[i][j]);
        }
    }
    memset(minn,0x3f,sizeof minn);
    memset(minn[0],0,sizeof minn[0]);
    for(int i=1;i<=n;++i){
        memset(f[i&1],0x3f,sizeof f[i&1]);
        for(int j=0;j<=(1<<m)-1;++j){
            for(int k=0;k<=(1<<m)-1;++k){
                work(i,j,thi);
                work(i-1,k,las);
                int res=0;
                for(auto p:thi){
                    res++;
                    for(auto q:las){
                        if(p==q){
                            res--;
                            break;
                        }
                    }
                }
                f[i&1][j][k]=min(f[i&1][j][k],minn[i-1][k]+res);
                minn[i][j]=min(minn[i][j],f[i&1][j][k]);
            }
        }
    }
    int ans=0x7fffffff;
    for(int j=0;j<=(1<<m)-1;++j){
        ans=min(ans,minn[n][j]);
    }
    cout<<ans<<'\n';
}

B.光线追踪

好题,有意思

转化成角度,那么就相当于是在做区间问题

首先射线只会碰到一个矩形的下或左边,因此剩下两个边可以直接忽略

假如射线会遇到多个矩形,那么我们肯定是选横/纵坐标最小的那个(哪个不变就看哪个,对于横线肯定是看纵坐标,竖线就是看横坐标)

所以问题就转化成了区间修改与单点查询最小值

因为询问的是矩形编号,所以应该是单点查询最小值编号

然后是一些小点,比如如果有一个是 \(0\) 应该返回另一个,或者两个最值相等应该返回编号大的(后来的会覆盖先来的),还有当给定的向量存在 \(0\) 时的处理

这个题还卡精度,卡精度我就直接搬有理数类套 map 了,离散化一下还是简单的,注意求斜率分母为 \(0\) 的情况

#include<bits/stdc++.h>
#define int long long
#include"../include/hdk/frac.h"
//https://www.cnblogs.com/HaneDaCafe/articles/18439223

using namespace std;
int n;
struct ques{
    int op;
    int xa,ya,xb,yb;
}q[100001];
hdk::frac a[400001];
int cnt,ree=1;
map<hdk::frac,int>mp;
bool cmp(const hdk::frac &A,const hdk::frac &B){
    return B<A;
}
inline int maped(pair<int,int> x){
    if(x.second==0) return 1;
    return mp[hdk::frac(x.first,x.second)];
}
struct __ans{
    int ans,minn;
};
struct stree{
    struct tree{
        int l,r;
        int pos,minn;
    }t[1000001];
    #define tol (id*2)
    #define tor (id*2+1)
    #define mid(l,r) mid=((l)+(r))/2
    void build(int id,int l,int r){
        t[id].l=l;t[id].r=r;
        t[id].pos=0;t[id].minn=0;
        if(l==r) return;
        int mid(l,r);
        build(tol,l,mid);
        build(tor,mid+1,r);
    }
    void change(int id,int l,int r,int val,int pos){
        if(l<=t[id].l and t[id].r<=r){
            if(t[id].pos==0 or t[id].minn>=val){
                t[id].pos=pos;
                t[id].minn=val;
            }
            return;
        }
        if(r<=t[tol].r) change(tol,l,r,val,pos);
        else if(l>=t[tor].l) change(tor,l,r,val,pos);
        else{
            change(tol,l,t[tol].r,val,pos);
            change(tor,t[tor].l,r,val,pos);
        }
    }
    __ans ask(int id,int pos){
        if(t[id].l==t[id].r){
            return {t[id].pos,t[id].minn};
        }
        __ans res={};
        if(pos<=t[tol].r) res=ask(tol,pos);
        else res=ask(tor,pos);
        if(t[id].pos==0) return res;
        if(res.ans==0 or t[id].minn<res.minn) return {t[id].pos,t[id].minn};
        else if(t[id].minn==res.minn and t[id].pos>res.ans) return {t[id].pos,t[id].minn};
        return res;
    }
};
stree x,y;
int cal(__ans x,__ans y,int xx,int yy){
    if(x.ans==0) return y.ans;
    if(y.ans==0) return x.ans;
    if(xx==0){
        if(q[x.ans].ya<q[y.ans].ya) return x.ans;
        return y.ans;
    }
    if(yy==0){
        if(q[x.ans].xa<q[y.ans].xa) return x.ans;
        return y.ans;
    }
    if(x.minn*xx==y.minn*yy){
        if(x.ans>y.ans) return x.ans;
        return y.ans;
    }
    if(x.minn*xx<y.minn*yy){
        return x.ans;
    }
    return y.ans;
}
signed main(){
    freopen("raytracing.in","r",stdin);
    freopen("raytracing.out","w",stdout);
    scanf("%lld",&n);
    for(int i=1;i<=n;++i){
        scanf("%lld",&q[i].op);
        if(q[i].op==1){
            scanf("%lld %lld %lld %lld",&q[i].xa,&q[i].ya,&q[i].xb,&q[i].yb);
            if(q[i].xa) a[++cnt]=hdk::frac(q[i].ya,q[i].xa);
            if(q[i].xb) a[++cnt]=hdk::frac(q[i].ya,q[i].xb);
            if(q[i].xa) a[++cnt]=hdk::frac(q[i].yb,q[i].xa);
        }
        else{
            scanf("%lld %lld",&q[i].xa,&q[i].ya);
            if(q[i].xa) a[++cnt]=hdk::frac(q[i].ya,q[i].xa);
        }
    }
    sort(a+1,a+cnt+1,cmp);
    for(int i=1;i<=cnt;++i){
        if(i==1 or !(a[i]==a[i-1])){
            mp[a[i]]=++ree;
        }
    }
    x.build(1,1,ree);
    y.build(1,1,ree);
    for(int i=1;i<=n;++i){
        if(q[i].op==1){
            int r1=maped({q[i].yb,q[i].xa});
            int r2=maped({q[i].ya,q[i].xb});
            int mid=maped({q[i].ya,q[i].xa});
            x.change(1,min(r2,mid),max(r2,mid),q[i].ya,i);
            y.change(1,min(r1,mid),max(r1,mid),q[i].xa,i);
        }
        else{
            __ans tmp1=x.ask(1,maped({q[i].ya,q[i].xa}));
            __ans tmp2=y.ask(1,maped({q[i].ya,q[i].xa}));
            cout<<cal(tmp1,tmp2,q[i].xa,q[i].ya)<<'\n';
        }
    }
}
posted @ 2024-10-07 21:13  HaneDaniko  阅读(15)  评论(0编辑  收藏  举报