学习打卡第D（16）天

在解决了一堆憨憨问题后终于可以写今天的博客啦
今天依旧很贪心！
先看题目吧
德国科学家总是对非洲野猴的抵抗力感到惊奇，因为他们发现在没有医疗条件的情况下，非洲野猴总是比其他所有野生动物少生病。最近的研究有了新发现，科学家Dr.Smith从非洲野猴身上发现了一种罕见的抗体，他猜测可能正是该罕见的抗体在帮助非洲野猴抵抗外来病毒的侵害。

Dr.Smith就立刻展开了对该抗体的研究。在初始的观察中Dr.Smith发现该抗体没什么特别，而且非常简单，因为抗体的每组基因只有一对基元（Dr.Smith把一组基因看成由若干对基元组成）。但是当Dr.Smith把病毒植入抗体所在的培养液后，奇迹出现了！那些简单的基因组通过不断地聚合（每个基因组两两合并生成新的基因组），最终所有的基因组合并成了一个非常庞大的基因组，而正是这庞大的基因组，因为聚合了所有原始基因组的优点，这庞大的基因组才可以慢慢地、逐个地去吃掉那些植入培养液的病毒。

下面是Dr.Smith在高倍显微镜下看到的抗体基因组基元聚合的大致过程：

图3－1 Dr.Smith在显微镜下看到的开始时有3个基因组的抗体聚合过程

Dr.Smith通过进一步观察和研究发现，抗体基因在聚合过程中似乎总是按照某个方法在进行，该方法能保证最终产生的基因组的基元对数量最多（每个基元对的存在能产生1u单位的生物能量），而在每次两个基因组聚合后所得到的新基因组的总的基元对由下面两部分相加组成：

1、每两个基因组一发生聚合，就产生一个额外的、未知的基元对。

2、当两个基因组聚合时，每个基因组中的每对基元都会与另一个基因组中的每对基元两两聚合产生一个新的基元对。

Dr.Smith还发现，抗体在每个时刻总是只有两个基因组会发生聚合，也就是说，每两个基因组的聚合都是依次进行的，而不是同时进行的。

虽然观察到聚合原理，但Dr.Smith即使在知道一开始基因组个数的前提下，还是无法统计最终聚合产生的那个庞大的基因组所具有的总能量有多大。现在他想请你编程来统计一下。
原题链接：https://www.luogu.com.cn/problem/P1974
我的答案：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
using namespace std;
int n; 
struct BigInt{
    int arr[3000];
    int len;
    BigInt(){
        memset(arr,0,sizeof(arr));
        len=0;
    }
}a[10100];
int que[50100],top,dep;
BigInt operator *(BigInt s,BigInt t){
    BigInt d;
    for(int i=1;i<=s.len;i++)
        for(int j=1;j<=t.len;j++){
            d.arr[i+j-1]+=s.arr[i]*t.arr[j];
            if(d.arr[i+j-1]>=10){
                d.arr[i+j]+=d.arr[i+j-1]/10;
                d.arr[i+j-1]%=10;
            }
        }
    d.len=s.len+t.len-1;
    for(int i=1;i<=d.len;i++){
        if(d.arr[i]>=10){
            d.arr[i+1]+=d.arr[i]/10;
            d.arr[i]%=10;
        }
        if(d.arr[d.len+1]>0)
            d.len++;
    }
    return d;
}
BigInt operator +(BigInt s,int t){
    s.arr[1]+=t;
    for(int i=1;i<=s.len;i++){
        if(s.arr[i]>=10){
            s.arr[i+1]+=s.arr[i]/10;
            s.arr[i]%=10;
        }
        if(s.arr[s.len+1]>0)
            s.len++;
    }
    return s;
}
void print(int x){
    for(int i=a[que[x]].len;i>=1;i--)
        printf("%d",a[que[x]].arr[i]);
}
int main(){
    freopen("test.in","r",stdin);
    freopen("test.out","w",stdout);
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        a[i].arr[1]=1;
        a[i].len=1;
        que[i]=i;
    }
    top=1,dep=n;
    while(top<dep){
        a[que[top]]=a[que[top]]*a[que[top+1]];
        a[que[top]]=a[que[top]]+1;
        que[++dep]=que[top];
        top+=2;
    }
    print(top);
    puts("");
    return 0;