School of Software Engineering, SCUT, 2021 记录

题目链接:Dashboard - SCUT_SE 21-22下半学年第七周训练div.2 - Codeforces

题解链接:文叔叔 - 传文件,找文叔叔(永不限速) 


C. Simple AniPop

类似洛谷【能量项链】那题,但是加了最后剩下的那个数字,所以稍微改改就可以。

#include <bits/stdc++.h>
using namespace std;
#define FOR(i, a, b) for (int i = (a); i <= (b); i++)
#define ROF(i, a, b) for (int i = (a); i >= (b); i--)
#define int long long
const int maxn = 1010;
int n; //珠子数量
int a[maxn], f[maxn][maxn];
//f[i][j]表示把[i+1,j-1]部分所有数字删掉,留下a[i]和a[j],所得的最大分数
//那么f[i][i+n]就表示刚好剩下a[i]的时候,消除得到的最大分数,再加上a[i]就是最终分数,所以f[i][i+n]+a[i]的最大值就是最终答案
signed main(){
    int maxx=0;
    cin >> n;
    for(int i=1; i<=n; i++){
        cin >> a[i]; maxx=max(maxx,a[i]);
        a[i+n] = a[i]; //环形记录
    }
    FOR(len,3,n+1){
        for(int i=1,j=i+len-1; j<=2*n; i++,j++)
        for(int k=i+1; k<=j-1; k++)
            f[i][j] = max(f[i][j], f[i][k]+f[k][j]+a[i]*a[k]*a[j]);
    }
    int ans = 0;
    FOR(i,1,n) ans=max(ans,f[i][i+n]+a[i]);
    cout<<ans;
}

E. Eigen Substring 

用后缀自动机,每加一个字符输出一次结果,结果是fail树中叶子节点对应的minlen的最小值。

用multiset可以维护,只需要维护对应的minlen即可,不用维护结点编号。

#include <bits/stdc++.h>
#define debug cout<<"!!!"<<endl;
#define FOR(i, a, b) for (int i = (a); i <= (b); i++)
#define ROF(i, a, b) for (int i = (a); i >= (b); i--)
#define pii pair<int,int>
#define ls p<<1
#define rs p<<1|1
using namespace std;
const int N = 2e6+5, inf=0x3f3f3f3f;
int n; string s;
multiset<int> st;
//SAM
int tot=1, las=1;
int len[N], fa[N], sz[N];
int tim[N]; //1表示有,0表示无
int to[N][26];

inline void ins(int c){
    int cur = ++tot, p = las; las = cur;
	tim[cur]=1;
    len[cur] = len[p]+1; sz[cur] = 1;
    for(; p && !to[p][c]; p=fa[p])
        to[p][c] = cur;
    if(!p) fa[cur] = 1;
    else {
        int q = to[p][c];
        if(len[q] == len[p]+1) {
			fa[cur] = q;
			if(tim[q]) tim[q]=0, st.erase(st.find(len[fa[q]]+1)); //删掉q的长度(注意语法,这里是用find是删除对应的迭代器,这样才能实现有相同数字时【只删其中一个】)
		}
        else{
            int cl= ++tot; len[cl] = len[p]+1;
			if(tim[q]) st.erase(st.find(len[fa[q]]+1)); //删掉原本的q的长度(注意语法,这里是用find是删除对应的迭代器,这样才能实现有相同数字时【只删其中一个】)
            memcpy(to[cl], to[q], sizeof(to[q]));
            fa[cl] = fa[q];
            fa[q] = fa[cur] = cl;
            for(; p && to[p][c]==q; p=fa[p]) to[p][c] = cl;
			if(tim[q]) st.insert(len[fa[q]]+1); //加入q'的长度
        }
    }
    st.insert(len[fa[cur]]+1); //加入当前点的长度
}
inline void solve(){
	cin>>n>>s;
	FOR(i,1,n) ins(s[i-1]-'a'), cout<<*st.begin()<<'\n';
}
signed main(){
    ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    int T; T=1; while(T--) solve();
}

另外,赛时使用的线段树维护,线段树底层结点中,非inf值是叶子节点的minlen,inf值是非叶子结点的值,用inf表示非法值。

#include <bits/stdc++.h>
#define debug cout<<"!!!"<<endl;
#define FOR(i, a, b) for (int i = (a); i <= (b); i++)
#define ROF(i, a, b) for (int i = (a); i >= (b); i--)
#define pii pair<int,int>
#define ls p<<1
#define rs p<<1|1
using namespace std;
const int N = 2e6+5, inf=0x3f3f3f3f;
int n; string s;
int mn[N<<2];
//SAM
int tot=1, las=1;
int len[N], fa[N], sz[N];
int leaf[N]; //1表示是叶子,0表示不是叶子
int to[N][26];

void build(int p,int l,int r){
    mn[p] = inf;
	if(l==r) return;
    int mid = (l+r)>>1;
    build(ls,l,mid);
    build(rs,mid+1,r);
}
void upd(int p,int l,int r,int cur,int val){
    if(l==r && l==cur){mn[p]=val; return;}
    int mid = (l+r)>>1;
    if(cur<=mid) upd(ls,l,mid,cur,val);
    else upd(rs,mid+1,r,cur,val);
    mn[p] = min(mn[ls],mn[rs]);
}
inline void ins(int c){
    int cur = ++tot, p = las; las = cur;
    leaf[cur] = 1;
    len[cur] = len[p]+1; sz[cur] = 1;
    for(; p && !to[p][c]; p=fa[p])
        to[p][c] = cur;
    if(!p) fa[cur] = 1;
    else {
        int q = to[p][c];
        if(len[q] == len[p]+1) {
			fa[cur] = q;
			if(leaf[q]) leaf[q]=0, upd(1,1,(int)2e6+2,q,inf);
		}
        else{
            int cl= ++tot; len[cl] = len[p]+1;
            // to[cl] = to[q];
            memcpy(to[cl], to[q], sizeof(to[q]));
            fa[cl] = fa[q];
            fa[q] = fa[cur] = cl;
            for(; p && to[p][c]==q; p=fa[p]) to[p][c] = cl;
			if(leaf[q]) upd(1,1,(int)2e6+2,q,len[fa[q]]+1);
        }
    }
    upd(1,1,(int)2e6+2,cur,len[fa[cur]]+1);
}
inline void solve(){
	cin>>n>>s;
    build(1,1,(int)2e6+2);
	FOR(i,1,n) ins(s[i-1]-'a'), cout<<mn[1]<<'\n';
}
signed main(){
    ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    int T; T=1; while(T--) solve();
}

G. Nim plus 

博弈题,可以用dp混过去。(思路稍微有点难懂,我还没有完全理解)

#include <bits/stdc++.h>
using namespace std;
#define FOR(i, a, b) for (int i = (a); i <= (b); i++)
#define ROF(i, a, b) for (int i = (a); i >= (b); i--)
const int N = 5010;
int n,m,a[2][N],f[2][N];

signed main(){
    ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    cin>>n>>m;
    FOR(j,0,1) FOR(i,1,m) cin>>a[j][i];
    FOR(i,1,n){
        for(int j=1; j<=m && a[0][j]<=i; j++)
            if(f[1][i-a[0][j]]==0) f[0][i]=1; //假如此时对面没办法赢,那我就能赢了
        for(int j=1; j<=m && a[1][j]<=i; j++)
            if(f[0][i-a[1][j]]==0) f[1][i]=1; //同理
    }
    if(f[0][n]==1) cout<<"Long Long nb!";
    else cout<<"Mao Mao nb!";
}

H. Treasure Hunt 

费用流问题,要用到最小费用可行流算法,求最大费用可行流。

第一次学费用流,说下可行流和最大流的区别。

在网络流中,随便一种流都叫做可行流,而总流量最大的流叫做最大流。不管是最大流还是可行流,都有可能有多种,他们对应的总费用可能不同。最小费用可行流就是随便一个流,求最小总费用;最小费用最大流,就是所有最大流中取最小总费用。而实际上,前者只需要在后者板子上改一行就可以了,所以板子几乎通用。

回到本题,每到达一个点都能获得对应金币,那么我们可以把这个金币看成费用,指向该点的所有边的费用就是这个金币数(这里转化比较巧妙),然后边的流量就是边允许走的次数(这里显然),跑一个最大费用可行流即可。

#include <bits/stdc++.h>
using namespace std;
#define debug cout<<"!!!"<<endl;
#define FOR(i, a, b) for (int i = (a); i <= (b); i++)
#define ROF(i, a, b) for (int i = (a); i >= (b); i--)
const int inf = 1<<30;
const int N = 110;
int n,m,a[N],val[N][N],cost[N][N],s,t, mincost=0;
int pre[N],flow[N],dis[N],vis[N];
vector<int> g[N];

int spfa(){
    FOR(i,1,n) dis[i] = inf; //离s点距离,初始是inf
    FOR(i,1,n) pre[i] = 0; //记录前驱,0表示一开始没有前驱
    dis[s] = 0; flow[s] = inf; vis[s] = 1; //超级原点的初始设置
    queue<int> q; q.push(s); //队列中放入原点
    while(q.size()){
        int u = q.front(); q.pop();
        vis[u] = 0; //出队,vis标记为0
        for(int v:g[u]){
            if(val[u][v] <= 0) continue; //流量<=0,实际上不连通
            if(dis[v] <= dis[u]+cost[u][v]) continue; //松弛失败,不做任何操作
            dis[v] = dis[u]+cost[u][v]; //进行松弛,更新距离
            pre[v] = u; //记录前驱(为了服务于mcmf函数,要记录路径)
            flow[v] = min(flow[u], val[u][v]); //记录可行的最大流量(val的前缀最小)
            if(!vis[v]) q.push(v), vis[v]=1; //当前点进入队列
        }
    }
	if(dis[t] == inf || dis[t]>=0) return 0; //如果无法到终点,或者到终点不赚钱或亏钱了,那以后就都不继续走了(证明略)
    //重点!!!dis[t]>=0就是【最小费用可行流】比起【最小费用最大流】需要加的唯一一句代码!!!
	else return 1; //能到终点,还能赚钱,这次走比不走好
}
void mcmf(){
    while(spfa()){
        int now = t;
        //maxflow += flow[t]; 记录最大流(这里其实跑的不是最大流,这个值没有意义,仅作为板子参考)
        mincost += flow[t]*dis[t]; //flow是本次流量,dis是这一条路上的【单位流量的费用】的总和
        while(now!=s){ //从终点一路往前,更新流量的信息
			val[pre[now]][now] -= flow[t]; //正向路,流量减少
			val[now][pre[now]] += flow[t]; //反向路,流量增加
            now = pre[now]; //往前走
        }
    }
}
signed main(){
    ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    cin>>n>>m; FOR(i,1,n) cin>>a[i]; //输入每个点赚的钱
    s=0; t=n; //超级原点和终点
    g[0].push_back(1); g[1].push_back(0); //超级原点和起点建立双向连接
    val[0][1] = inf; // 0->1流量无限
    cost[0][1] = -1*a[1]; //注意为了求最大费用流,把权值全部倒转了,正向是-1倍
    cost[1][0] = a[1]; //反向则是1倍

    FOR(i,1,m){
        int u,v,w; cin>>u>>v>>w;
        g[u].push_back(v); g[v].push_back(u); //双向建边
        val[u][v] = w; //注意只有单向流量
        cost[u][v] = -1*a[v]; //注意为了求最大费用流,把权值全部倒转了,正向是-1倍
        cost[v][u] = a[v]; //反向则是1倍
    }
    mcmf(); //跑mcmf
    cout<<-1*mincost; //输出最大费用可行流大小,注意乘上-1(因为代码里跑的是最小费用可行流)
}

K. Number Puzzle 

待更....


L. Long Long Wanna Buy 

大致题意:给定n行m列个(si,pi)对,在每行中选一个(si,pi)对,求使得\frac{\sum p}{s(min)}最小时的\sum p

注意到s的范围很小,只有1000,所以可以枚举s(min),把所有所有si >= s(min)对应的最小pi加起来即可,记录一下后缀最小就很容易处理。

#include <bits/stdc++.h>
using namespace std;
#define int long long
#define debug cout<<"!!!"<<endl;
#define FOR(i, a, b) for (int i = (a); i <= (b); i++)
#define ROF(i, a, b) for (int i = (a); i >= (b); i--)
const int N = 1010;
int n,m,s[N][N],p[N][N],mnp[N][N];
int pos[N];

signed main(){
    ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    cin>>n>>m;
    FOR(i,1,n) FOR(j,1,m) cin>>s[i][j]>>p[i][j]; //输入信息
    FOR(i,1,n) mnp[i][m+1]=0x3f3f3f3f; //mnp初始化
    FOR(i,1,n) ROF(j,m,1) mnp[i][j]=min(mnp[i][j+1],p[i][j]); //求p的后缀最小值

    double res=1e9; int ans=0;
    FOR(smn,1,1000){
        FOR(i,1,n) while(s[i][pos[i]]<smn && pos[i]<=m) pos[i]++; //移动指针
        FOR(i,1,n) if(pos[i]==m+1) break; //检查越界
        int tmp=0;
        FOR(i,1,n) tmp+=mnp[i][pos[i]]; //当前的sum_p
        if((double)tmp/smn < res) {res = (double)tmp/smn; ans=tmp;}
    }
    cout<<ans;
}

M. Camouflage 

按题意模拟即可,代码难度不高,思维难度为0

#include <bits/stdc++.h>
using namespace std;
#define FOR(i, a, b) for (int i = (a); i <= (b); i++)
#define ROF(i, a, b) for (int i = (a); i >= (b); i--)
const int dx[4]={1,0,-1,0}, dy[4]={0,1,0,-1};
const int N = 1005, inf=0x3f3f3f3f;
int n,m,d,mp[N][N],ans[N][N];
int vis[N][N], cnt=0;

void dfs(int x,int y){
    // if(cnt==-1) return;
    if(mp[x][y]==-1) return; //遇到-1结束
    if(x==1 || x==n || y==1 || y==m) {cnt=-1;} //到边界,是-1
    if(vis[x][y]) return;
    vis[x][y]=1; //标记遍历过
    if(mp[x][y]==1){
        if(cnt!=-1) cnt++;
        FOR(i,0,3){
            int xi=x+dx[i], yi=y+dy[i];
            if(!vis[xi][yi]) dfs(xi,yi);
        }
    }
}
void dfs2(int x,int y){
    if(mp[x][y]==-1) return;
    mp[x][y]=-1; //修改mp
    FOR(i,0,3){
        int xi=x+dx[i], yi=y+dy[i];
        if(xi<1 || xi>n || yi<1 || yi>m || mp[xi][yi]==-1) continue;
        dfs2(xi,yi);
    }
}
signed main(){
    ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    cin>>n>>m>>d; char ch;
    FOR(i,1,n) FOR(j,1,m){
        cin>>ch;
        if(ch=='.') mp[i][j]=1;
        else mp[i][j]=-1;
    }
    FOR(i,1,n) FOR(j,1,m){
        dfs(i,j);
        if(mp[i][j]==1) ans[i][j]=cnt;
        cnt=0;
    }
    FOR(i,1,n) FOR(j,1,m){
        if(0<ans[i][j] && ans[i][j]<d) dfs2(i,j);
    }
    FOR(i,1,n){
        FOR(j,1,m){
            if(mp[i][j]==-1) cout<<'#';
            else cout<<'.';
        }
        cout<<'\n';
    }
}

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值