2018 计蒜之道 初赛 第五场

Posted hl666

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了2018 计蒜之道 初赛 第五场相关的知识,希望对你有一定的参考价值。

这次的比赛没有现场打,而是等到了今天才来补。

主要是因为那时候和HHHOJ上的比赛冲突了,所以就没写。

这次前三题的难度都比较低,但是就是一个T4要莫比乌斯反演。又是不可食用的。

好了我们开始看题。

A. 贝壳找房搬家

这道题刚开始看的时候没看懂题意,觉得T1就是这种立体几何的题目,有种想的感觉。

因为我认为这个方块可以不规则地想怎么放就怎么放的,其实题目中有一句话:

我们可以把这堆箱子看成一个\(x \times y \times z\) 的长方体。

什么?刚开始只能是长方体吗?好吧好像还是不可做的样子。

然后一看时限:5000ms,那不就……暴力一下?

我们枚举最短边\(a(1<=a<=\sqrt[3]n)\),然后是次短边\(b(1<=b<=\sqrt n)\)

然后我们只需要算出第三边\(h(h为整数)\),然后讨论一下谁是长,谁是宽,谁是高。

然后直接更新最大最小值即可。注意这里长宽都是原来+2的结果,但高是+1。

CODE

#include<cstdio>
using namespace std;
int t,n,h,v,MAX,MIN;
inline char tc(void)
{
    static char fl[100000],*A=fl,*B=fl;
    return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
    x=0; char ch=tc();
    while (ch<‘0‘||ch>‘9‘) ch=tc();
    while (ch>=‘0‘&&ch<=‘9‘) x=x*10+ch-‘0‘,ch=tc();
}
inline void write(int x)
{
    if (x/10) write(x/10);
    putchar(x%10+‘0‘);
}
inline int min(int a,int b)
{
    return a<b?a:b;
}
inline int max(int a,int b)
{
    return a>b?a:b;
}
int main()
{
    //freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
    register int a,b; read(t);
    while (t--)
    {
        read(n); MAX=-1; MIN=1e9;
        for (a=1;a*a*a<=n;++a)
        for (b=1;b*b<=n;++b)
        {
            if (n%(a*b)) continue; h=n/a/b; 
            MIN=min(MIN,(a+2)*(b+2)*(h+1)-n); MAX=max(MAX,(a+2)*(b+2)*(h+1)-n);
            MIN=min(MIN,(a+2)*(b+1)*(h+2)-n); MAX=max(MAX,(a+2)*(b+1)*(h+2)-n);
            MIN=min(MIN,(a+1)*(b+2)*(h+2)-n); MAX=max(MAX,(a+1)*(b+2)*(h+2)-n);
        }
        write(MIN); putchar(‘ ‘); write(MAX); putchar(\n);
    }
    return 0;
}

B. 贝壳找房算数(简单)

这还是比较水的,直接\(O(n^2logn)\)判断即可。

因为分解以及gcd都是log的,所以不会超时。

CODE

#include<cstdio>
using namespace std;
int n,a,b,ans;
long long k;
inline int f(int x)
{
    int tot=1;
    while (x) tot*=x%10,x/=10;
    return tot;
}
inline int gcd(int m,int n)
{
    return n?gcd(n,m%n):m;
}
int main()
{
    register int i,j; scanf("%d%lld",&n,&k);
    for (i=1;i<=n;++i)
    for (j=1;j<=n;++j)
    {
        int a=f(i),b=f(j);
        if (!a||!b) continue;
        if (gcd(a,b)<=k) ++ans;
    }
    printf("%d",ans);
    return 0;
}

C. 贝壳找房算数(中等)

这个也是很Easy的,因为我们注意到像124,142和241等等这些数它们的贡献是一样的。

因此我们发现有用的\(f(i)\)很少,所以我们统计一下所有\(f(i)=x\)的数分别有多少。

这里由于数的分布不均匀,因此我们上map大法

这里主意一下对于map的遍历,需要用指针和一个奇奇怪怪的STL来完成

这遍历非人哉,我才不会告诉你我是从网上找的方法

CODE

#include<cstdio>
#include<map>
using namespace std;
const int mod=998244353;
map <int,int> t;
map <int,int>::iterator a,b;
int n,w;
long long k,ans;
inline int f(int x)
{
    int tot=1;
    while (x) tot*=x%10,x/=10;
    return tot;
}
inline int gcd(int m,int n)
{
    return n?gcd(n,m%n):m;
}
int main()
{
    register int i; scanf("%d%lld",&n,&k);
    for (i=1;i<=n;++i)
    ++t[f(i)];
    for (a=t.begin();a!=t.end();++a)
    for (b=a;b!=t.end();++b)
    {
        int x1=a->first,y1=a->second,x2=b->first,y2=b->second;
        if (!x1||!x2) continue;
        if (gcd(x1,x2)>k) continue;
        if (x1^x2) ans=(ans+(long long)y1*y2*2)%mod; else ans=(ans+y1*y2)%mod;
    }
    printf("%lld",ans);
    return 0;
}

以上是关于2018 计蒜之道 初赛 第五场的主要内容,如果未能解决你的问题,请参考以下文章

2017 计蒜之道 初赛 第五场 B. UCloud 的安全秘钥(简单)

2017 计蒜之道 初赛 第五场 C. UCloud 的安全秘钥(中等)

2017 计蒜之道 初赛 第五场 UCloud 的安全秘钥(中等)

2017 计蒜之道 初赛 第五场 D. UCloud 的安全秘钥(困难)

计蒜之道第五场C题

2018 计蒜之道 初赛 第三场