jyzy noip模擬賽5.22-1,3

這樣的兩道題

樣例輸入

5
1 1 1 1 4

樣例輸出

2

樣例輸入

5 10

樣例輸出

8

顯然的水題

關於第一題

貪心一下，顯然是最左邊最右邊放一個度爲一的，所有的度大於1的在中間排成鏈

$ans=n-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）
*/