BZOJ 1002([FJOI2007]轮状病毒-基尔霍夫矩阵)

1002: [FJOI2007]轮状病毒

Time Limit: 1 Sec   Memory Limit: 162 MB
Submit: 1402   Solved: 758
[ Submit][ Status][ Discuss]

Description

 给定n(N<=100),编程计算有多少个不同的n轮状病毒。

Input

第一行有1个正整数n。

Output

将编程计算出的不同的n轮状病毒数输出

Sample Input

3

Sample Output

16

HINT

Source

[ Submit][ Status][ Discuss]




不想吐槽了,kUbuntu各种坑人……

各种卡机子不说,好像Google的词库也没法用。

Codeblock根本进不了Win盘(没权限估计)

话说看PPT那种换一页等3s真不是人受的

话说我到底有没有写题解啊!(1G内存党伤不起)

至于基尔霍夫矩阵什么的我才不知道呢

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<functional>
#include<cmath>
#include<cctype>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define RepD(i,n) for(int i=n;i>=0;i--)
#define MAXN (100+10)
#define F (10000)
int n;
struct Highn
{
    int a[10000],len;
    Highn(){len=0;memset(a,0,sizeof(a));}
    Highn(int b)
    {
        len=0;memset(a,0,sizeof(a));
        while (b) a[++len]=b%F,b/=F;
        if (len==0) len=1;
    }
    friend Highn operator*(int b,Highn a)
    {
        For(i,a.len) a.a[i]*=b;
        For(i,a.len) a.a[i+1]+=a.a[i]/F,a.a[i]%=F;
        if (a.a[a.len+1]) a.len++;
        return a;
    }
    friend Highn operator+(Highn b,Highn a)
    {
        Highn c;
        c.len=max(a.len,b.len);
        For(i,c.len)
        {
            c.a[i]+=a.a[i]+b.a[i];
            c.a[i+1]+=c.a[i]/F;
            c.a[i]%=F;
        }
        c.len++;
        while (!c.a[c.len]) c.len--;
        return c;
    }
    friend Highn operator-(Highn a,Highn b)
    {
        Highn c;
        c.len=max(a.len,b.len);
        For(i,c.len)
        {
            c.a[i]+=a.a[i]-b.a[i];
            if (c.a[i]<0) c.a[i]+=F,c.a[i+1]--;
        //  c.a[i+1]+=c.a[i]/F;
        //  c.a[i]%=F;
        }
        while (!c.a[c.len]) c.len--;
        return c;
    }

    void print()
    {
        printf("%d",a[len]);
        ForD(i,len-1)
        {
            printf("%04d",a[i]);
        }
        puts("");
    }
}f[MAXN];
int main()
{
//	freopen(".in","r",stdin);
//	freopen(".out","w",stdout);
    scanf("%d",&n);
    f[1]=1,f[2]=5;
    if (n<=2) f[n].print();
    else
    {
        Fork(i,3,n) f[i]=3*f[i-1]-f[i-2]+2;
        f[n].print();
    }


	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值