ZROI week4

于 2019-01-06 20:31:00 发布
阅读量145 收藏
点赞数
原文链接：http://www.cnblogs.com/akoasm/p/10230248.html
版权

考试

前言

起晚了，大概10点才开始看T1，被别人问了T2有点懵逼。

和

这题看了就A掉了，感觉很像原题的样子，是我的错觉吗？？

串串

某神仙有个\(O(n)\)做法问了我一下，我当时也没怎么想（因为我觉得出题人不会把On的题给到\(n^2\) )

结果人家是对的

这个写了个大讨论，然后讨论崩了，成功0分，WTF

写完这个已经11点半了，下次再也不在这种讨论题还拿不了分的题上浪费时间。。。

简单函数

这个题挺有意思，首先打表算了一下他的前几位，发现都是\(2 \times n - (n的因数和)\)

用线性筛晒出约数和直接计算即可。

其实分析一下式子也可以得出一个结论，每一个都是\(n - (非n的因数和)\)，可以等于\(n + n - (有n的因数和)\)

。。。为啥我要打表？

复杂度\(O(n)\)

写完+拍完已经\(11:50\)

简单MST

感觉好像一个原题啊，以前绝对做过的那种。

大概就是一个边取模的完全图生成树。

把边权枚举出来转到一个相同系下，做\(Kruskal\)即可。

没时间写了，写完不能一边对就打了个暴力，\(30pts\)。

总结

100 + 0 + 100 + 30

吐槽

事后\(Logey\)跟我吐槽了一下某个\(A\)掉\(T4\)的人代码和\(std\)非常像，一些没有必要的东西也在上面，拜托某些人抄代码能不能有点水平...

还有就是为啥\(std\)会让学生看到？？

upd:今天来打了一遍\(T4\)，一遍A掉了QAQ

作业

HNOI 2015 菜肴制作

解法

贪心不好做，为什么？

如果我们要取出1，那么就要取出1的所有前驱，然后一路取，复杂度\(O(n^2)\)

但是我们发现虽然不能快速直到第一个数字该填什么，但是我们直到最后一个填什么，所有没有出边的点就是有可能的最后一个数字，而根据之前的原则，选最大的。

倒着用堆维护这个序列即可。

复杂度\(O(n \; log \; n)\)

#include <bits/stdc++.h>
using namespace std;
const int MAXN = 1e5 + 10;

priority_queue < int > q;
int d[MAXN];
int res[MAXN];
int st[MAXN];

vector <int> e[MAXN + 10];

set <int> s[MAXN + 10];

int n,m,T;
int cnt;
int x,y;
int tag;
int read () {
    int q=0,f=1;char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-')f=-1;ch=getchar();
    }
    while(isdigit(ch)){
        q=q*10+ch-'0';ch=getchar();
    }
    return q*f;
}

int main () {
    cin >> T;
    while(T--) {
        tag = 0;
        n = read(),m = read();
        for(int i = 0;i < MAXN; ++i) {
            e[i].clear();
            s[i].clear();
        }
        memset(d,0,sizeof d);
        for(int i = 0;i < m; ++i) {
            x = read(),y = read();
            --x,--y;
            if(s[x].find(y) == s[x].end()) {
                s[x].insert(y);
                ++d[x];
                e[y].push_back(x);
            }
        }
        cnt = 0;
        for(int i = 0;i < n; ++i) {
            if(!d[i]) {
                q.push(i);
            }
        }
        while(!q.empty()) {
            int now = q.top();
            q.pop();
            res[cnt++] = now;
            for(int j = 0;j < e[now].size(); ++j) {
                if(!(--d[e[now][j]])) q.push(e[now][j]);
            }
        }
        if(cnt < n) {
            cout<<"Impossible!"<<endl;
        }
        else {
            for(int i = n - 1;i >= 0; --i) {
                cout<<res[i] + 1<<' ';
            }
            cout<<endl;
        }
    }
    return 0;
}

bzoj 1601 灌水

解法

建造最小生成树

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=400;
const int maxm=1000000;
int n,m;
struct edge{int u,v,w;} e[maxm];
bool cmp(edge a,edge b)
{
    return a.w<b.w;
}
bool flag[maxn][maxn];
int fa[maxn];
int find(int x)
{
    if(x==fa[x]) return x;
    else return fa[x]=find(fa[x]);
}
int kruscal()
{
    int ans=0;
    sort(e+1,e+1+m,cmp);
    for(int i=1;i<=n;++i) fa[i]=i;
    int cnt=n+1;
    for(int i=1;i<=m;++i)
    {
        int t1=find(e[i].u);
        int t2=find(e[i].v);
        if(t1!=t2)
        {
            if(cnt==1) break;
            fa[t1]=t2;
            ans+=e[i].w;
            cnt--;
        }
    }
    return ans;
}
int main()
{
    cin>>n;
    int tmp=n; 
    for(int i=1;i<=n;++i)
    {
        cin>>e[i].w;
        e[i].u=0;
        e[i].v=i;
    } 
    for(int i=1;i<=n;++i)
    for(int j=1;j<=n;++j)
    {
        int a;
        cin>>a;
        if(i!=j&&!flag[i][j]) e[++tmp].u=i,e[tmp].v=j,e[tmp].w=a,flag[i][j]=1,flag[j][i]=1;
    }
    m=tmp;
    cout<<kruscal()<<endl;
    return 0;
}

bzoj 2882 工艺

解法

用类似于后缀自动机的办法把序列搞到\(parent\) 树上求出循环最小序列即可。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int MAXN = 6e5 + 10;
int ntot;
int lst;
int fa[MAXN << 1] = {-1};
int len[MAXN << 1];
int siz[MAXN << 1];
int a[MAXN];
int n;
map <int,int> mp[MAXN << 1];
int read () {
    int q=0,f=1;char ch=getchar();
    while(!isdigit(ch)) {
        if(ch=='-')f=-1;ch=getchar();
    }
    while(isdigit(ch)){
        q=q*10+ch-'0';ch=getchar();
    }
    return q*f;
}

void extend(int c) {
    int np = ++ntot;
    len[np] = len[lst] + 1;
    siz[np] = 1;
    int now;
    for(now = lst;now != -1 and !mp[now].count(c);now = fa[now]) {
        mp[now][c] = ntot;
    }
    if(now == -1) {
        fa[np] = 0;
    }
    else {
        int q = mp[now][c];
        if(len[q] == len[now] + 1) {
            fa[np] = q;
        }
        else {
            int nq = ++ntot;
            len[nq] = len[now] + 1;
            fa[nq] = fa[q];
            mp[nq] = mp[q];
            fa[q] = nq;
            fa[np] = nq;
            for(int j = now;j != -1 and mp[j][c] == q;j = fa[j]) {
                mp[j][c] = nq;
            }
        }
    }lst = np;
}

int main () {
    //freopen("1.in","r",stdin);
    //freopen("ams.out","w",stdout);
    n = read();
    for(int i = 1;i <= n; ++i) {
        a[i] = read();
        extend(a[i]);
    }
    for(int i = 1;i <= n; ++i) {
        extend(a[i]);
    }
    int p;
    p = 0;
    for(int i = 1;i <= n; ++i) {
        cout<<mp[p].begin() -> first<<' ';
        p = mp[p][mp[p].begin() -> first];
    }
    return 0;
}

bzoj 3732 Network

解法

和货车运输基本一样，原题做法也可以。

现在考虑一个新的东西——最小生成树重构树。

建出虚点后查询\(LCA\)的权值即可。

#include <bits/stdc++.h>
using namespace std;
const int MAXN = 1e5 + 10;
const int INF = 0x3f3f3f3f;
#define LOG 20

int read () {
    int q=0,f=1;char ch=getchar();
    while(!isdigit(ch)) {
        if(ch=='-')f=-1;ch=getchar();
    }
    while(isdigit(ch)){
        q=q*10+ch-'0';ch=getchar();
    }
    return q*f;
}

int fa[MAXN * 2][22];
struct edge {
    int u,v,w;
}e[MAXN << 1];
int cnt;
int tot;
int Cnt;
int f[MAXN << 1];
int dep[MAXN << 1];
int son[MAXN][2];
int v[MAXN];
int n,m,k;
int x,y;
#define C continue
bool cmp(edge a,edge b) {
    return a.w < b.w;
}

int find(int x) {
    return x == f[x] ? x : f[x] = find(f[x]);
}

void kruskal() {
    sort(e + 1,e + Cnt + 1,cmp);
    for(int i = 1;i <= m; ++i) {
        int x = e[i].u;
        int y = e[i].v;
        int l = find(x);
        int r = find(y);
        if(l == r) C;
        son[++cnt][0] = l; //方点 
        son[cnt][1] = r;
        f[f[x]] = f[f[y]] = fa[f[x]][0] = fa[f[y]][0] = cnt;
        v[cnt] = e[i].w;
    }
}

void dfs(int now) {
    if(!son[now][0] and !son[now][1]) return;
    dep[son[now][0]] = dep[son[now][1]] = dep[now] + 1;
    dfs(son[now][0]);
    dfs(son[now][1]);
}

int find(int x,int y) {
    if(dep[x] < dep[y]) swap(x,y);
    for(int i = LOG;i >= 0; --i) {
        if(dep[fa[x][i]] >= dep[y]) {
            x = fa[x][i];
        }
    }
    if(x == y) return x;
    for(int j = LOG;j >= 0; --j) {
        if(fa[x][j] != fa[y][j]) {
            x = fa[x][j];
            y = fa[y][j];
        }
    }
    return fa[x][0];
}

int main () {
    n = read(),m = read(),k = read();
    cnt = n;
    for(int i = 1;i <= 2 * n + 10; ++i) f[i] = i;
    for(int i = 1;i <= m; ++i) {
        e[++Cnt].u = read(),e[Cnt].v = read();
        e[Cnt].w = read();
    }
    kruskal();
    dep[cnt] = 1;
    dfs(cnt);
    for(int j = 1;j <= LOG; ++j) {
        for(int i = 1;i <= (n << 1); ++i) {
            fa[i][j] = fa[fa[i][j - 1]][j - 1];
        }
    }
    while(k--) {
        x = read(),y = read();
        //cout<<find(x,y)<<endl;
        printf("%d\n",v[find(x,y)]);
    }
    return 0;
}

转载于:https://www.cnblogs.com/akoasm/p/10230248.html

aijiu8213
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值