SDUACM暑期集训周赛（五）

最新推荐文章于 2023-06-22 07:54:13 发布
最新推荐文章于 2023-06-22 07:54:13 发布
阅读量295 收藏
点赞数
分类专栏： SDU_ACM_Practice 文章标签： gcd phi
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Yancy_/article/details/77435729
版权

BZOJ 2186 沙拉公主的困惑

[Problem]

询问n!中与m!互质的个数（n > m）

[Solution]

ans = n! / m ! * phi(m!)

显然[1, m!] 此区间中有phi(m!)个与m!互质的数

[m! + 1, m! + m!] 同样有phi(m!)个 gcd(m! + i, m !) = gcd(i, m!)

[2 * m! + 1, 3 * m!] phi(m!) 个
…….

[Code]

#include <cstdio>  
#define ll long long  
#define maxn 10000000  
#define maxk 1000000  
using namespace std;  
ll p[maxk], R, inv[maxn+100], N, M, tot, fac[maxn+100],  
    ans[maxn+100], T;  
bool no[maxn+100];  

void getinv()  
{  
    int i, m=(R<maxn?R:maxn); inv[1]=1;  
    for(i=2;i<m;i++)inv[i]=(ll)(R-R/i)*inv[R%i]%R;  
}  

void shai()  
{  
    int i, j;  
    for(i=2;i<=maxn+1;i++)  
    {  
        if(!no[i])p[++tot]=i;  
        for(j=1;i*p[j]<=maxn;j++)  
        {  
            no[i*p[j]]=true;  
            if(i%p[j]==0)break;  
        }  
    }  
}  

void getfac()  
{  
    int i;  
    fac[1]=1;  
    for(i=2;i<=maxn;i++)fac[i]=fac[i-1]*i%R;  
}  

void getans()  
{  
    int i, j;  
    ll t;  
    ans[1]=1;  
    for(i=1;i<=tot;i++)  
    {  
        t=(ll)ans[p[i]-1]*(p[i]-1)%R*inv[p[i]%R]%R;  
        for(j=p[i];j<p[i+1];j++)ans[j]=t;  
    }  
}  

int main()  
{  
    int i;  
    scanf("%d%d",&T,&R);  
    shai();  
    getinv();  
    getfac();  
    getans();  
    while(T--)  
    {  
        scanf("%d%d",&N,&M);  
        printf("%lld\n",fac[N]*ans[M]%R);  
    }  
    return 0;  
}

BZOJ 2186 POJ 2992 Divisors

[Solution]
把每个数字分解质因数，然后把结果存到数组中，前缀和累加一下，对于计算cnm时用幂进

行运算

[Code]

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define N 1000
int p[500], top = 0, ans[500];
int num[N][N], n, m, sum[N][N];
bool isp[N];
long long f[N][N];
void cal(int n, int m)
{
    for(int i = 1; i <= top; i++)  ans[i] = 0;
    if (n - m < m) m = n - m;
    if (f[n][m]) return;
    for(int j = 1; j <= top; j++) {
        ans[j] = sum[n][j] - sum[m][j] - sum[n - m][j];
    }
}
void init()
{
    for(int i = 2; i <= 431; i++) {
        if (isp[i])  continue;
        p[++top] = i;
        for(int j = i; j <= 431; j += i) isp[j] = true;
    }
    for(int i = 2; i <= 431; i++)  {
        int x = i;
        for(int j = 1; j <= top; j++){
            while (x % p[j] == 0) {
                num[i][j]++;
                x /= p[j];
            }
        }
    }
    for(int i = 2; i <= 431; i++) {
        for(int j = 1; j <= top; j++) {
            sum[i][j] += sum[i - 1][j] + num[i][j];
        }
    }
}
int main()
{
   // freopen("b.in", "r", stdin);
    init();
    while(~scanf("%d%d", &n, &m)) {
        cal(n, m);
        if (f[n][m]) {
            cout << f[n][m] << endl;
            continue;
        }
        long long  cur = 1LL;
        for(int i = 1; i <= top; i++) cur *= (ans[i] + 1);
        f[n][m] = cur;
        f[n][n - m] = cur;
        cout << cur << endl;
    }
    return 0;
}

POJ 1696 Space Ant

[Solution]

每次向左偏转的角度尽可能地小，判断点是不是在左边通过叉积计算

[Code]

#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
#define N 2000
const double eps = 1e-9;

int n, top;
bool f[N];
struct Point {
    double x, y;
    int num;
} a[N];
int ans[N];
bool cmp(Point a, Point b)
{
    if (a.y == b.y)   return a.x < b.x;
    return a.y < b.y;
}

int Cross(Point p0, Point p1, Point p2) {
    double tot =  (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);
    if (fabs(tot) < eps)  return 0;
    return tot < 0 ? -1 : 1;
}

double dis(Point u, Point v)
{
    return sqrt((u.x - v.x) * (u.x - v.x) + (u.y - v.y) * (u.y - v.y));
}
bool lcmp(Point u, Point v1, Point v2)
{
    int d = Cross(u, v2, v1);
    if (!d)  return dis(u, v2) < dis(u, v1);
    return d > 0;
}

int main()
{
   //freopen("b.in", "r", stdin);
    int T;
    scanf("%d", &T);
    while(T--)  {
        scanf("%d", &n);
        for(int i = 1; i <= n; i++) scanf("%d%lf%lf", &a[i].num, &a[i].x, &a[i].y);
        sort(a + 1, a + n + 1, cmp);
        for(int i = 1; i <= n; i++)  f[i] = false;
        ans[1] = a[1].num;  f[1] = true;
        Point lastu;
        lastu.x = 0; lastu.y = a[1].y;
        Point u = a[1];
        int top = 1;
        for(int i = 2; i <= n; i++) {
            int minj = -1;
            for(int j = 2; j <= n; j++) {
                if (!f[j]  && (Cross(lastu, u, a[j]) >= 0))  {
                    if (minj < 0 || lcmp(u, a[minj], a[j]))  minj = j;
                }
            }
            if (minj < 0)  break;
            ans[++top] = a[minj].num;
            f[minj] = true;
            lastu = u;
            u = a[minj];
        }
        printf("%d", top);
        for(int i = 1; i <= top; i++)  printf(" %d", ans[i]);
        printf("\n");
    }

    return 0;
}

POJ 1319 Pipe Fitters

[Solution]

对于第一种放置方式，直接是(int)a * (int)b
对于第二种放置方式，最后一行可以防止(int)b, 倒数第二层取决于b的小数部分
需要高度与实际高度的比例可以推出
注意第二种放置方式始终有一层的高度是无法优化的
注意第二种放置方式把a当行一次，把b当行一次

[Code]

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
int cal(double a, double b)
{
        int d1 = (int)b, d2 = (int) (b + 0.5);
        d2--;
        int k = (int)((a - 1) * 2 / sqrt(3.0)); // important
        if (a >= 1) k++;
        int ans2 = (k / 2) * (d1 + d2);
        if (k % 2) ans2 += d1;
        return ans2;
}
int main()
{
   // freopen("b.in", "r", stdin);
    double a, b;
    while (~scanf("%lf%lf", &a, &b))  {
        int ans1 = (int)a * (int)b;
        int ans2 = max(cal(a, b), cal(b, a));
        if (ans1 >= ans2)
            printf("%d grid\n", ans1);
        else
            printf("%d skew\n", ans2);
    }
    return 0;
}

HDU 3746 Cyclic Nacklace

[Solution]

kmp判断

首先考虑虽然我们可以添加到左边，也可以添加到右边，但实际上等效于添加到一边，因为最终会连成一个环。

因此我们只思考在字符串的右面加上一些字符，考虑最后一位的next值，len-next[len]为循环节长度d， next[len] % d 便是最后多的那一部分

[Code]

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define N 200500
char c[N];
int nxt[N];
void get_next(char* s, int* nxt)
{
    int n = strlen(s);
    int j = 0;
    nxt[0] = nxt[1] = 0;
    for(int i = 1; i < n; i++){
        while (j > 0 && s[i] != s[j]) j = nxt[j];
        if (s[i] == s[j]) j++;
        nxt[i + 1] = j;
    }
}
int main()
{
  //  freopen("b.in", "r", stdin);
    int T;
    scanf("%d", &T);
    while (T--)  {
        scanf("%s", c);
        get_next(c, nxt);
        int len = strlen(c);
        int ans = len - nxt[len];
        int t = nxt[len] % ans;
        if (!t && nxt[len]) t = ans;
        printf("%d\n", ans - t);
    }
    return 0;
}

BZOJ 2186 沙拉公主的困惑

[Solution]

只想说尽可能不会再把数字转成字符串处理了
[Code]

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define N 3010
#define M 2
#define ll long long
int cnt = 0;
int id[N * 34], ch[N * 34][2], tim[N * 34];
int pow[33];
int n, m;
int a[N];
void init()
{
    cnt = 0;
    for(int i = 0; i < M; i++) ch[0][i] = 0;
    id[0] = 0;
    tim[0] = -1;
}
void insert(int x, int k)
{
    int u = 0;
    for(int i = 30; i >= 0; i--)  {
        int v = (x & (1 << i)) > 0 ? 1 : 0;
        if (!ch[u][v]) {
            ch[u][v] = ++cnt;
            for(int j = 0; j < M; j++) ch[cnt][j] = 0;
            id[cnt] = 0;
            tim[cnt] = 0;
        }
        u = ch[u][v];
        tim[u]++;
    }
    id[u] = k;
}
void update(int x, int t)
{
    int u = 0;
    for(int i = 30; i >= 0; i--)  {
        int v = (x & (1 << i)) > 0 ? 1 : 0;
        u = ch[u][v];
        tim[u] += t;
    }
}
int cal(int x)
{
    int u = 0;
    for(int i = 30; i >= 0; i--)  {
        int v = (x & (1 << i)) > 0 ? 1 : 0;
        if (ch[u][1 - v] && tim[ch[u][1 - v]]) u = ch[u][1 - v];
        else u = ch[u][v];
    }
    return id[u];
}
int main()
{
   // freopen("b.in", "r", stdin);
    int T;
    scanf("%d", &T);
    while(T--)  {
        init();
        scanf("%d", &n);
        for(int i = 1; i <= n; i++)  {
            scanf("%d", a + i);
            insert(a[i], i);
        }
        int ans = 0;
        for(int i = 1; i <= n; i++)  {
            update(a[i], -1);
            for(int j = i + 1; j <= n; j++) {
                update(a[j], -1);
                int k = cal(a[i] + a[j]);
                ans = max(ans, ((a[i] + a[j]) ^ a[k]));
                //printf("%d %d %d\n", i, j, k);
                update(a[j], 1);
            }
            update(a[i], 1);
        }
        printf("%d\n", ans);
    }
    return 0;
}
YuHsin
关注
专栏目录
SDUACM暑期集训周赛（二）
Yancy_的博客
08-08 285
POJ 3613 Cow Relays
2016年暑假集训周赛#1题解
qq_33094929的博客
07-28 513
2016年暑假集训周赛#1        周赛链接 Problem A：跑跑卡丁车系列之游戏下载 题目大意：给出n个数,然后Q个询问,每次询问输入一个数字,输出n个数里面有多少个数小于等于它(总的减去大于这个数的) 思路：先排序,然后用upper_bound返回大于的第一元素 #include #include #include using namespace std; in
【FJNUMCS暑假集训周赛#1】
这太Imba了
07-29 396
2016年暑假集训周赛#1Problem A: 跑跑卡丁车系列之游戏下载Time Limit: 1000 MS Memory Limit: 257792 KB 64-bit interger IO format: %lld Java class name: Main某天，金斌在实验室玩跑跑卡丁车，众多群巨也一起加入其中。蛋蛋发现这是一个虐爆群巨的好机会，所以他也准备去
第五次周赛总结
yeyuluo的博客
11-28 106
第五次周赛总结 A.防AK题目——超难系列 B.Fenoix超厌恶xxx C.lwm学姐的木屋 D.Fenoix的趣事 E.Fenoix的签到题 F.cgy学长爱吃糖果 G.来呀，贪心呀 H.说是栈你会做么？ I.我太难了 J.我也太难了 K.Fenoix的针对 L.来自Fenoix的刁难 ...
周赛补题【2022年湘潭大学暑假集训选拔赛2022.7.2】
WSY444的博客
07-08 355
jiayouya
2016暑期集训---周赛(水题)
Code_J_xer的博客
07-21 495
今天被这个题前前后后折磨了好几次。其实最开始就想到可能需要用数组来存数。但是做题过程中却偏偏只是一步一步地放大储存容量。我感觉这可能是程序员的一个小小职业病。由于计算机的存储空间以及计算性能是有限的，所以一直精打细算时间跟空间，一种“不浪费时间和空间的思想”让我在这题上栽倒了好几次。第一次，用int存数。显然，要挂。其实这次WA是不意外的。 第二次，于是把所有的数都换成了long long。结果。
CSDN周赛回顾
文盲的专栏
05-25 1090
得之我幸，失之我命
集训两次周赛总结
大神养成中.....
07-20 1253
集训过去一半了，进行了两次周赛。成绩还算满意，但感觉还有待提高。 第一次比赛：http://acm.hdu.edu.cn/webcontest/contest_show.php?cid=5131 这次比赛被老大坑了，大概之做了最后的一个半小时，水出来两道题目。。。后面下来做了一下，都是简单题目，还写了一个 题解： 第二次比赛：http://acm.hdu.edu.cn/webcontest
CSDN第十五期周赛
weixin_50089904的博客
12-19 340
第三题有点考察思维了，既然说不能成环那我们为了保证建立最多的边，只要使两方边尽可能相等就可以了，(n/2) * (n - n/2)，这就是答案了。也可能是第一次使用类似的oj，不太适应他的评测机制，总感觉自测比较麻烦。不过评测机倒是挺快的。第四题是之前给新生培训的一道原题，dp+贪心的解法就可以做出来。第一次在csdn周赛中取得满分的成绩，先来说一些做题体验吧。总的来说这次体验还是蛮不错的，也希望自己能取得更好的成绩。第二题只需记录最大的数字，最后输出最大数+1即可。这次的题感觉比以往几期稍微简单一点。
【蓝桥杯集训·周赛】AcWing 第96场周赛
你的博客
03-26 1025
（2）因为能够被2520整除的是2520的倍数，所以当前n是2520的多少倍（下取整），即存在多少个能够被2520整除的数。,dn 以及三个整数 k,a,b。有一个容量为 n 的背包和 m+1 种物品，每种物品都有无限多个。接下来 m 行，每行包含四个整数 li,hi,vi,wi。一个整数，表示 [1,n] 范围内完美数的个数。第 i 种物品的体积为 vi，价值为 wi。第一行包含 5 个整数 n,k,a,b,q。第一行包含四个整数 n,m,v0,w0。注意审题，不要将 n,m 的含义弄混。
leetcode周赛难度-leetcode:leetcode
07-06
周赛难度LeetCode 题解 一题目列表 细绳 # 标题 标签 困难 解决方案 描述 1 细绳 中等的 2 细绳 中等的 3 细绳 中等的 4 细绳 中等的 DP # 标题 标签 困难 解决方案 描述 1 DP 中等的 2 DP 中等的 3 DP 中等的 4 DP ...
leetcode周赛得分3-mstcBlog:微博
07-07
标题 "leetcode周赛得分3-mstcBlog:微博" 暗示了这是一个与LeetCode周赛相关的项目，其中可能包含作者在解决LeetCode编程挑战时的思路、代码实现以及可能的成绩分享。LeetCode是一个在线平台，提供各种算法题目来...
leetcode周赛积分-leetcode:leetcode刷题记录
07-07
leetcode周赛积分 leetcode 推荐 No. Problems Solutions Tag 32 DP, Stack 周赛 (weekly-content) No. Solutions Completion Date 2020.07.05 刷题记录 No. Problems Solutions Completion Date 329 2020.07.26 16 ...
CSDN 周赛 59 期
最新发布
文盲的专栏
06-22 497
一如既往的填空题，以及新鲜出炉的带坑编程题
BZOJ 刷题总结
Yancy_的博客
01-29 1478
背景 辣鸡的人总要想法自救，便产生了寒假学些新算法、在Bzoj刷些题的想法。一来为明年省赛做准备…压力不小；二来寒假也可以有些事情做。 PS 1.不定时更新做题的思路和吐槽 2.按照hzw刷题顺序训练，具体依照BZOJ题表 3.希望寒假能够至少刷够100道题目吧，在此立个FLAG，希望别被青岛的妖风吹跑 Problem &amp; Solution BZOJ 1003 ...
SDU省赛选拔-ACM ICPC 2010–2011, NEERC, Northern Subregional Contest
Yancy_的博客
04-14 1252
A. Alien Communication Masterclass【题目】给定序列A、B，构造一个等式，使得对于任意的A[i]进制都满足，对于任意的B[i]都 不满足【分析i】同一个式子，不同进制的结果如果不相同，一般是在进位上，首先我们可以想到如果对于A序列的进制都满足呢，可以写成x1 * x2 * x3 *x4 = 0这种，每种进制可以构造一个等于0的式子，对于A[i]，可以构造（10 -
sduacm16级寒假训练 自测
Yancy_的博客
02-16 584
POJ 3579 Median【sol】显然，每个点与其他的点形成的差是有序的，这样相当于有n组升序的序列，询问中位数。二分中位数，然后判断比他小的数的个数，复杂度n* logn*logn#include<cstdio> #include<iostream> #include<cstring> #include<cmath> #include<algorithm> using namespace
Educational Codeforces Round 37 所有题目整理！
Yancy_的博客
02-05 548
A.Water The Garden int ct = a / p; if (a % p != 0) ++ct; 等价于 ct = (a + p - 1) / p; B. Tea Queue solution：按照题意模拟即可 C. Swap Adjacent Elements 题意：给定1-n的一个排列，存在一些i，其中1 思路：首先考虑数字1，假设其位置为pos1，必然
2017暑期集训Day 11 背包
Yancy_的博客
07-21 498
A - Bottles[Problem]有n个苏打水的瓶子，每瓶中装有一定量苏打水，给你每瓶的体积，现在像把这些苏打水装到尽量少的瓶子中，并且移动尽可能少的苏打水。[Solution]f[i]代表凑齐i体积需要瓶子的最小数量，g[i]代表盛有的苏打水体积，在满足f[i]尽可能小的情况下,g[i]最大[Code]#include<cstdio> #include<iostream> #include<
重庆八中周赛题解：编程竞赛算法分析
本次重庆八中周赛涉及了多种算法和问题解决策略，主要包括文本排列、矩阵覆盖和区间子序列处理等主题。下面分别对这些知识点进行详细解析： **P1 文本左右对齐** 此题考察的是字符串排列和贪心算法的应用。在处理...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值