jyzy noip模拟赛5.22-1,3

最新推荐文章于 2020-06-20 18:04:18 发布

startaidou

最新推荐文章于 2020-06-20 18:04:18 发布

阅读量192

点赞数

分类专栏：数学文章标签：数学

本文链接：https://blog.csdn.net/startaidou/article/details/90473345

版权

数学专栏收录该内容

10 篇文章 0 订阅

订阅专栏

这样的两道题

在这里插入图片描述样例输入

5
1 1 1 1 4

样例输出

在这里插入图片描述样例输入

5 10

样例输出

显然的水题

关于第一题

贪心一下，显然是最左边最右边放一个度为一的，所有的度大于1的在中间排成链

$\sum_{i=1}^{n}(max(2,d_{i})-2)$

不过考场上我不是这么写的。。。。。不过思路差不多

#include<bits/stdc++.h>
using namespace std;
int n,d[25],len=0,a[25];
int main()
{
	freopen("tree.in","r",stdin);
	freopen("tree.out","w",stdout);
	scanf("%d",&n);
	if(n<0) cout<<"hsmtql"; 
	int i=1;
	while(i<=n)
	{
		scanf("%d",&d[i]);
		++i;
	}
	sort(d+1,d+n+1);
	i=1;
	while(i<=n)
	{
		if(d[i]!=1) break;
		++i;
	}
	--i;
	if((i<2)&&(d[n]!=1))
	{
		printf("-1\n");
		return 0;
	}
	int j=i+1;
	while(j<=n)
	{
		a[++len]=d[j];
		++j;
	}
	j=1;
	i-=2;
	while(j<=len)
	{
		a[j]-=2;
		i-=a[j];
		++j;
	}
	if(i==0)
	{
		printf("%d\n",len+1);
	}
	else printf("-1\n");
	return 0;
}
/*
1996年：东方灵异传（TOH1）
1997年：东方封魔录（TOH2）
1997年：东方梦时空（TOH3）
1998年：东方幻想乡（TOH4）
1998年：东方怪绮谈（TOH5）
2002年：东方红魔乡（TOH6）
2003年：东方妖妖梦（TOH7）
2004年：东方萃梦想（TOH7.5）
2004年：东方永夜抄（TOH8）
2005年：东方花映冢（TOH9）
2005年：东方文花帖（TOH9.5）
2007年：东方风神录（TOH10）
2008年：东方绯想天（TOH10.5）
2008年：东方地灵殿（TOH11）
2009年：东方星莲船（TOH12）
2009年：东方非想天则（TOH12.3）
2010年：东方文花帖DS（TOH12.5）
2010年：东方三月精（TOH12.8）
2011年：东方神灵庙（TOH13）
2013年：东方心绮楼（TOH13.5）
2013年：东方辉针城（TOH14）
2014年：弹幕天邪鬼（TOH14.3）
2014年：东方深秘录（TOH14.5）
2015年：东方绀珠传（TOH15）
2017年：东方凭依华（TOH15.5）
2017年：东方天空璋（TOH16）
*/

关于第二题

估计所有人都会这道题

$\begin{gathered} \begin{bmatrix} f_{1} \\ f_{0} \end{bmatrix} \begin{pmatrix} \begin{bmatrix} 1 & 1 \\ 1 & 0 \end{bmatrix} \end{pmatrix}^{n}= \begin{bmatrix} f_{n+1} \\ f_{n} \end{bmatrix} \end{gathered}$

所以快速幂一下

不过我用了其他的小玩意

#include<bits/stdc++.h>
using namespace std;
char a[30003000];
long long p;
int len=0;
struct ma
{
    long long a[3][3];
    ma()
    {
        memset(a,0,sizeof(a));
    }
    ma operator*(ma x)
    {
        ma s;
        register int i=1,j=1,k=1;
        while(i<=2)
        {
            j=1;
            while(j<=2)
            {
                k=1;
                while(k<=2)
                {
                    s.a[i][j]=(s.a[i][j]+a[i][k]*x.a[k][j])%p;
                    ++k;
                }
                ++j;
            }
            ++i;
        }
        return s;
    }
    ma operator^(register long long k)
    {
        ma s;
        ma e;
        s.a[1][1]=a[1][1];
        s.a[1][2]=a[1][2];
        s.a[2][1]=a[2][1];
        s.a[2][2]=a[2][2];
        e.a[1][1]=e.a[2][2]=1;
        while(k)
        {
            if(k&1) e=e*s;
            k>>=1;
            s=s*s;
        }
        return e;
    }
};
inline long long gcd(register long long n,register long long m)
{
    while(m^=n^=m^=n%=m);
    return n;
}
inline long long lcm(register long long n,register long long m)
{
    return n/gcd(n,m)*m;
}
inline long long power(register long long n)
{
    ma a,ans;
    a.a[1][1]=a.a[1][2]=a.a[2][1]=1;
    a.a[2][2]=0;
    ans=a^n;
    return ans.a[1][2];
}
long long pi[10001000],k[10001000];
inline long long geet(register long long p)
{
    register long long s=sqrt(p);
    register int i=2,tot=0;
    while(i<=s)
    {
    	if(!(p%i))
    	{
        	pi[++tot]=i;
            k[tot]=1;
        	while(!(p%i))
            {
                p/=i;
                k[tot]*=i;
            }
    	}
    	++i;
    }
    i=1;
    while(i<=tot)
    {
        k[i]/=pi[i];
        ++i;
    }
    if(p!=1)
    {
        k[++tot]=1;
        pi[tot]=p;
    }
    i=1;
    while(i<=tot)
    {
        if(pi[i]==2) k[i]*=3;
        else if(pi[i]==3) k[i]*=5;
        else if(pi[i]==5) k[i]*=20;
        else if((pi[i]%5==1)||(pi[i]%5==4)) k[i]*=pi[i]-1;
        else k[i]*=(pi[i]+1)<<1;
        ++i;
    }
    register long long ans=k[1];
    i=2;
    while(i<=tot)
    {
        ans=lcm(ans,k[i]);
        ++i;
    }
    return ans;
}
int main()
{
	freopen("fibonacci.in","r",stdin);
	freopen("fibonacci.out","w",stdout);
    register long long m=0,n=0;
    scanf("%s%lld",a+1,&p);
    len=strlen(a+1);
    m=0;
    n=0;
    if(p==1)
    {
        printf("%lld\n",n);
        return 0;
    }
    m=geet(p);
    register int i=1;
    while(i<=len)
    {
        n=((n<<3)+(n<<1)+(a[i]^48))%m;
        ++i;
    }
    ++n;
    n%=m;
    if(!n) printf("%lld\n",n);
    else if(n==1) printf("%lld\n",n);
    else printf("%lld\n",power(n));
    return 0;
}
/*
1996年：东方灵异传（TOH1）
1997年：东方封魔录（TOH2）
1997年：东方梦时空（TOH3）
1998年：东方幻想乡（TOH4）
1998年：东方怪绮谈（TOH5）
2002年：东方红魔乡（TOH6）
2003年：东方妖妖梦（TOH7）
2004年：东方萃梦想（TOH7.5）
2004年：东方永夜抄（TOH8）
2005年：东方花映冢（TOH9）
2005年：东方文花帖（TOH9.5）
2007年：东方风神录（TOH10）
2008年：东方绯想天（TOH10.5）
2008年：东方地灵殿（TOH11）
2009年：东方星莲船（TOH12）
2009年：东方非想天则（TOH12.3）
2010年：东方文花帖DS（TOH12.5）
2010年：东方三月精（TOH12.8）
2011年：东方神灵庙（TOH13）
2013年：东方心绮楼（TOH13.5）
2013年：东方辉针城（TOH14）
2014年：弹幕天邪鬼（TOH14.3）
2014年：东方深秘录（TOH14.5）
2015年：东方绀珠传（TOH15）
2017年：东方凭依华（TOH15.5）
2017年：东方天空璋（TOH16）
*/

startaidou

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
jyzy noip模拟赛5.22-1,3

这样的两道题样例输入51 1 1 1 4样例输出2样例输入5 10样例输出8显然的水题关于第一题贪心一下，显然是最左边最右边放一个度为一的，所有的度大于1的在中间排成链ans=n−1−∑i=1n(max(2,di)−2)ans=n-1- \sum_{i=1}^{n}(max(2,d_{i})-2)ans=n−1−i=1∑n(max(2,di)−2)不过 ...
复制链接

扫一扫