2018 计蒜之道 初赛 第五场

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

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

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

好了我们开始看题。

A. 贝壳找房搬家

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

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

我们可以把这堆箱子看成一个(x imes y imes 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('
');
	}
	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;
}