Educational Codeforces Round 101 (Rated for Div. 2) D. Ceil Divisions 思维 + 根号数-CSDN博客

本文链接：https://blog.csdn.net/m0_51068403/article/details/113975904

这篇博客探讨了一种针对数组操作的算法，目标是通过不超过n+5次操作，将所有元素变为1和2的组合。作者提出以n的平方根为分母进行上取整操作，有效地减少数组中的大数值，并给出了具体的步骤和操作次数计算。文章还包含了一个完整的C++代码示例，展示了如何实现这一策略。

摘要由CSDN通过智能技术生成

传送门

题意： 给一个数组 $a_i=i$ ，每次可以进行操作 $a_x=\left \lceil \frac{a_x}{a_y} \right \rceil$ ，操作不能超过 $n + 5$ 次，最终需要把数组中的数变成 $n - 1$ 个 $1$ 和一个 $2$ 。

思路： 比较显然的是数组中原来的 $1$ 和 $2$ 是不需要操作的，可以把 $2$ 当做最终的 $2$ ，让后可以以 $a_n$ 为分母，让前面 $> 2$ 的数都除 $a_n$ ，但是这个样子的话， $a_n$ 就需要用 $2$ 来消除，但是操作明显不够。那么我们考虑一下 $\sqrt{n}$ ,显然我们可以用 $\sqrt{n}$ 当做分母， $a_n=\left \lceil \frac{a_n}{a_{\sqrt{n}}} \right \rceil$ 进行两次即可将 $a_n$ 化成 $1$ 。而 $\sqrt{n}=447$ ，用 $2$ 来消除的话还是有点不够，我们可以继续开根号，我们发现从 $2 e 5$ 开根号开到 $> 2$ 的最后一个，最多只有5个根号数，假设为a b c d e 。对于e我们可以 $a_e=\left \lceil \frac{a_e}{a_d} \right \rceil$ ，进行两次。对于每个根号数都这样操作，最多需要10次，让后对于不是根号数的直接以 $a_n$ 作为分母，一步到位。当根号数为5的时候，操作数为 $10 + n - 2 - 5 = n + 3$ ，满足题意。
还有注意是上取整，开根号的时候注意判断是否为平方数，我为了方便直接把每个平方根都+1了。

//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;

//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }

typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;

const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;

int n;
int s[20],tot;
vector<PII>v;
bool st[N];

int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);

    int _; scanf("%d",&_);
    while(_--)
    {
        scanf("%d",&n);
        int nn=n; v.clear(); tot=0;
        while(nn>2) st[nn]=true,s[++tot]=nn,nn=(int)sqrt(nn)+1;
        for(int i=3;i<=n;i++)
        {
            if(st[i]) continue;
            v.pb({i,n});
        }
        for(int i=1;i<tot;i++)
        {
            v.pb({s[i],s[i+1]});
            v.pb({s[i],s[i+1]});
        }
        v.pb({s[tot],2});
        v.pb({s[tot],2});
        printf("%d\n",(int)v.size());
        for(int i=0;i<(int)v.size();i++) printf("%d %d\n",v[i].X,v[i].Y);
        for(int i=1;i<=tot;i++) st[s[i]]=false;
    }

	return 0;
}
/*

*/