习题日常第十六练

1麦森数（题目链接）

这个题是个快速幂和高精度以及一点点数论的组合。我一开始发现根据高精度判断位数是不可能的，于是看了下别人的题解，可以通过log函数来计算位数。剩下的用高精度乘法加快速幂就解决了。代码如下。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
#include<sstream>
#define ll long long int
const int mod=100003;
using namespace std;
int a[1001],b[1001],c[1001],d[1001];
void cheng()
{
    memset(c,0,sizeof(c));
    ll i,j,k;
    for(i=0;i<500;i++)
        for(j=0;j<500;j++)
    {
        c[i+j]+=b[j]*a[i];
    }
    for(i=0;i<500;i++)
    {
        c[i+1]+=c[i]/10;
        c[i]%=10;
    }
    for(i=0;i<500;i++)
        b[i]=c[i];
    return;
}
void bian()
{
    memset(c,0,sizeof(c));
    ll i,j,k;
    for(i=0;i<500;i++)
        for(j=0;j<500;j++)
    {
        c[i+j]+=a[i]*a[j];
    }
    for(i=0;i<500;i++)
    {
        c[i+1]+=c[i]/10;
        c[i]%=10;
    }
    for(i=0;i<500;i++)
        a[i]=c[i];
    return;
}
void mypow(ll n)
{
    while(n)
    {
        if(n&1)
        {
           cheng();
        }
        bian();
        n>>=1;
    }
}
int main()
{
    ios::sync_with_stdio(false);
    ll i,j,l,k,n,m;
    cin>>n;
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    memset(c,0,sizeof(c));
    m=(ll)(n*log10(2)+1);
    a[0]=2;
    b[0]=1;
    cout<<m;
    mypow(n);
    b[0]--;
    for(i=499;i>=0;i--)
    {
        if((i+1)%50==0)
            cout<<endl<<b[i];
        else
            cout<<b[i];
    }
  return 0;
}

2御阪美琴（题目链接）

这个题用到点贪心的思路吧。首先对于这些数，顺序实际不重要，因为每一堆实际是由自己调整的。对于每一个大的数，肯定是由一个相对大的数分出来的更好，因为这个大的数被分以后可以确保所需要剩下的数中最大的数有可能被分出来。用优先队列进行模拟。代码如下。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
#include<sstream>
#define ll long long int
const int mod=100003;
using namespace std;
ll a[100100],b[100100],c[100100],d[100100];
priority_queue<ll>p1,p2;
int main()
{
    ios::sync_with_stdio(false);
    ll i,j,l,k,n,m,ans=0;
    cin>>n>>m;
    for(i=0;i<m;i++)
    {
        cin>>a[i];
        ans+=a[i];
        p2.push(a[i]);
    }
    if(ans!=n)
        cout<<"ham";
    else
    {
        ll flat=0;
        p1.push(n);
        while(!p2.empty())
        {
            ll x=p1.top();
            ll y=p2.top();
            if(p1.size()>p2.size())
            {
                flat=1;
                break;
            }
            if(x<y)
            {
                flat=1;
                break;
            }
            else if(x>y)
            {
                p1.pop();
                p1.push(x/2);
                p1.push(x-x/2);
            }
            else
            {
                p1.pop();
                p2.pop();
            }
        }
        if(p1.empty()&&p2.empty())
        {
            cout<<"misaka";
        }
        else
            cout<<"ham";
    }
  return 0;
}

3子串（题目链接）

这个题就当我复习kmp算法的一个板子题吧。对于2到16每一个进制进行转换为字符串，然后用kmp算法查看是否为子串即可。另外，直接用string里的find函数可能会更方便。具体代码如下。

#include<cmath>
#include <iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<cstring>
#include<math.h>
#include<stack>
#include<algorithm>
#include<queue>
#include<bitset>
#include<sstream>
#define ll long long int
const int mod=100003;
using namespace std;
int a[1050100];
string s,s1;
string chp(ll n,ll k)
{
    string sp;
   int i=0,lp;
   while(n)
   {
       a[i++]=n%k;
       n/=k;
   }
   lp=i;
   for(i=lp-1;i>=0;i--)
   {
       if(a[i]>=10)
        sp+=a[i]-10+'A';
       else
        sp+=a[i]+'0';
   }
   return sp;
}
bool kmp(string x,string y)
{
    memset(a,-1,sizeof(a));
    int i,j,k;
    a[0]=-1;
    i=0;
    k=-1;
    while(i<x.length()-1)
    {
        if(k==-1||x[i]==x[k])
        {
            i++;
            k++;
            a[i]=k;
        }
        else
            k=a[k];
    }
    i=0;
    j=0;
    int len1=x.length();
    int len2=y.length();
    while(j<len2)
    {
        if(x[i]==y[j]||i==-1)
        {
            i++;
            j++;
        }
        else
        {
            i=a[i];
        }
        if(i==len1)
            return 1;
    }
    return 0;
}
int main()
{
    ios::sync_with_stdio(false);
    int i,j,l,k,n,m,ans=0;
    cin>>n>>s;
    for(k=2;k<=16;k++)
    {
        string s2;
        for(i=1;i<=n;i++)
        {
            s2+=chp(i,k);
        }
       if(s2.length()>=s.length())
       {
           if(kmp(s,s2))
           {
               ans=1;
               break;
           }
       }
    }
    if(ans==1)
        cout<<"yes";
    else
        cout<<"no";
  return 0;
}