2019.11.4模拟赛

T1 奇因数之和

定义\(F(n)\)为\(n\)的最大奇因数,例如\(F(1) = 1\),\(F(6) = 3\),\(F(12) = 3\)。
输入\(m\),求\(\sum\limits_{i = 1}^{m}F(i)\)。
\(m \leq 10^{100}\)
分析一下我们发现,如果\(i\)是偶数,\(F(i) = F(i / 2)\),这样的话就可以分治求了。
这题难度在高精,谁考场上写大高精啊

r = open('sigma.in','r')
w = open('sigma.out','w')
ans = 0
def sove(*t):
    global ans
    x = int(t[0])
    if x == 0:
        return None
    if x % 2 == 1:
        ans += (x + 1) * (x + 1) // 4
    else:
        ans += x *x // 4
    sove(x // 2)
n = r.read()
sove(n)
w.write(str(ans))

当然可以不递归求

r = open('sigma.in', 'r')
w = open('sigma.out', 'w')
ans=  0
n = int(r.read())
while(n != 0):
    if(n & 1):
        ans += (n + 1) * (n + 1) // 4
    else:
        ans += n * n // 4
    n //= 2
w.write(str(ans))

以下是c++的代码。高精类过大没有放。
高精类模板

signed main()
{
    freopen("sigma.in", "r", stdin);
    freopen("sigma.out", "w", stdout);
    GNUM n;
    GNUM zero(0);
    cin >> n;
    GNUM ans(0);
    while (!(n == zero))
    {
        if (n.opd())
            ans = ans + (n + 1) * (n + 1) / 4;
        else
            ans = ans + n * n / 4;
        n = n / 2;
    }
    cout << ans << endl;
    return 0;
}

T2 激光塔阵

小强的激光塔阵共有\(P\)座激光塔,它们被建造在一个$N \times M \(的矩形网格上。每座激光塔都有自己的坐标\)(x,y)\(,其中\)x(1 \leq x \leq n)\(和\)y(1 \leq y \leq m)\(均为整数。当两座激光塔的\)x\(坐标和\)y\(坐标至少一个相同时,这两座塔就可以产生联系。若激光塔\)A\(与\)B\(有联系,\)B\(与\)C\(也有联系,则可以认为\)A\(与\)C$也有联系。但是现在激光塔并不一定两两都有联系,所以他想新建一些新的激光塔使得激光塔两两都有联系。
你要计算并输出小强最少需要建造多少个激光塔。
显然并查集维护一下有几个连通块就可以了。

int main()
{
    freopen("laser.in", "r", stdin);
    freopen("laser.out", "w", stdout);
    poread(n), poread(m), poread(p);
    for(register int i = 1; i <= p; ++i)
        poread(data[i].x), poread(data[i].y), data[i].id = i;
    for(register int i = 1; i <= p; ++i)
        fa[i] = i;
    sort(data + 1, data + 1 + p, cmp1);
    for(register int i = 2; i <= p; ++i)
    {
        if(data[i - 1].x == data[i].x)
        {
            register int x = find(data[i - 1].id), y = find(data[i].id);
            if(x == y)
                continue;
            fa[y] = x;
        }
    }
    sort(data + 1, data + 1 + p, cmp2);
    for(register int i = 2; i <= p; ++i)
    {
        if(data[i - 1].y == data[i].y)
        {
            register int x = find(data[i - 1].id), y = find(data[i].id);
            if(x == y)
                continue;
            fa[y] = x;
        }
    }
    register int ans = 0;
    for(register int i = 1; i <= p; ++i)
    {
        if(find(i) == i)
        {
            ++ans;

        }
    }
    printf("%d", ans - 1);
}

T3 深入虎穴

这题是个语文题不放题面了。大概就是每次求次短路都从次短路进行转移。
只放垃圾跑的贼慢的\(spfa\)了。

inline void sspfa()
{
    queue<int> q;
    memset(d, 0x3f, sizeof(d));
    memset(md, 0x3f, sizeof(md));
    for(register int i = 1; i <= k; ++i)
    {
        d[ek[i]] = md[ek[i]] = 0;
        g[ek[i]] = mg[ek[i]] = ek[i];
        q.push(ek[i]);
        v[ek[i]] = 1;
    }
    while(q.size())
    {
        register int x = q.front();
        q.pop();
        v[x] = 0;
        for(register int i = head[x]; i; i = e[i].nxt)
        {
            register int y = e[i].ver;
            if(d[y] > md[x] + e[i].edge)
            {
                if(g[y] != x)
                    md[y] = d[y], mg[y] = g[y];
                d[y] = md[x] + e[i].edge;
                g[y] = x;
                if(!v[y])
                    q.push(y), v[y] =  1;
            }
            else if(md[y] > md[x] + e[i].edge && g[y] != x)
            {
                md[y] = md[x] + e[i].edge;
                mg[y] = x;
                if(!v[y])
                    q.push(y), v[y] =  1;
            }
        }
    }
}

原文地址:https://www.cnblogs.com/Shiina-Rikka/p/11819212.html

时间: 2024-09-30 04:32:03

2019.11.4模拟赛的相关文章

2017.6.11 校内模拟赛

题面及数据及std(有本人的也有原来的) :2017.6.11 校内模拟赛 T1 自己在纸上模拟一下后就会发现 可以用栈来搞一搞事情 受了上次zsq 讲的双栈排序的启发.. 具体就是将原盘子大小copy一下排个序 用两个指针维护两个数组(原数据 和 排序后的数据), 即分为1数据和2数组 将小于1指针指向的数据的2数组中的数据全部压入栈中 后进行消除, 将栈栈顶元素与当前1数组中的1指针指向的元素进行比较 相同则消除 后重复过程 直至指针超过N 后判断一下是否两个指针都超过了N... #incl

2019.11.30训练赛总结

2019.11.30训练赛总结 Codeforces Round #499 (Div. 2) 总的来说是一场很不愉快的比赛.漏洞百出. 对于A题,其实没有什么技术含量,只是写的时候忘记了边界的情况,导致出现错误. B题,一定程度上考验了思维,既然从正面做不行,那么我可以反着来,既然求不可以正向求出答案,那我可以把答案枚举带进去看是否符合条件啊,如果数据范围再大点的话还可以二分. D题,也反应出了自己的一个漏洞,就是懒得一步一步去推数据,代数据进去.当出现bug的时候,最直接,最有效的办法,就是把

11.27 模拟赛

并没有人做的模拟赛... 出题人hx,,, T1:就是上一道矩阵乘法 数学题 T2: 一个数列中 一个区间满足,存在一个k(L <= k <= R),并且对于任意的i (L <= i <= R),ai都能被ak整除 这样的一个特殊区间 [L, R]价值为R - L 想知道序列中所有特殊区间的最大价值是多少,而有多少个这样的区间呢 这些区间又分别是哪些呢 输出每个区间的L 思路: 用两个ST表分别求一段区间的gcd和最小值 然后可以二分答案 check的时候枚举左端点,判断在这段区间

2017 9 11 noip模拟赛T2

#include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int N=205; int map[N][N]; int d[N],tag[N],book[N],f[N]; int n,m; void work(int x) { memset(d,63,sizeof(d)); memset(book,0,sizeof(book)); memset(f,0,sizeof(

2019.10.24模拟赛赛后总结

本文原创,如果有不到位的地方欢迎通过右下角的按钮私信我! A.Icow Player 题目描述 被无止境的农活压榨得筋疲力尽后,Farmer John打算用他在MP3播放器市场新买的iCow来听些音乐,放松一下.FJ的iCow里存了N(1 <= N <= 1,000)首曲子,按1..N依次编号.至于曲子播放的顺序,则是按一个Farmer John自己设计的算法来决定: * 第i首曲子有一个初始权值R_i(1 <= R_i <= 10,000). * 当一首曲子播放完毕,接下来播放的

11.4 模拟赛

终于AK了,虽然第三题主要是搞月想出来的 T1: n个1*1的小方块,把这些小方块拼成一个图形,使这个图形周长最小 思路: 枚举拼成长方形的长为i,宽为n/i 可得面积 (i+n/i+(bool)(n%i))*2 1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 #include<cstdlib> 5 #include<cmath> 6 #include<cstr

2017/11/3模拟赛

块(block)[问题描述]拼图达人小 C 手里有 n 个 1*1 的正方形方块, 他希望把这些方块拼在一起, 使得拼出的图形周长最小, 要求方块不能重叠. 擅长拼图的小 C 一下就求出了这个周长, 顺便他想考考你会不会求.[输入格式]多组数据, 第一行一个正整数 T, 表示数据组数.接下来 T 行, 每行一个正整数 n, 表示方块数.[输出格式]输出 T 行, 每行一个正整数, 表示答案.[样例输入]3 4 1122[样例输出]81420[数据范围]对于 20%的数据, n<=20:对于 40

2017/11/1模拟赛

磁星(magnet)[题目描述]在 B 城呆惯了的小 H 决定去外太空溜达一圈.人类现已发现并开发的星球(包括小 H 所在的星球)有 n 个,并且在这 n 个星球之中,人们发现了 m 对两个星球的关系.关系"xy"表示的是星球 x 对星球 y 有 1 一个单位的引导力,由于引导力还具有传递性,如果星球 x 对星球 y 能有恰好 a 个单位的引导力,星球y 对星球 z 能有恰好 b 个单位的引导力,那么星球 x 对星球 z 就能有恰好 a+b 个单位的引导力. 换言之,星球 x 对星球

11.2 模拟赛

T1: tarjan裸不能再裸了 1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 #include<cstring> 5 #include<cstdlib> 6 #include<cmath> 7 #include<vector> 8 #include<queue> 9 #define inf 2147483611 10 #defi