• BZOJ [FJOI2007]轮状病毒


    分析:

    我其实想到链状怎么做了,没有往下多想。一直想用组合做。

    最后看到CLJ的题解,才发现我已经想对了一半。。。

    对于一个链来说,设Dp[i]为长度为i的链和一个“中心”的生成树数量,可以Dp之。。
    然后考虑环,枚举1节点所在的环的长度,设为Len,那么这个环的位置有Len种情况,这个环与中心连接也有Len种情况,那么这个节点对答案的贡献就是Len^2*Dp[n-Len]。。

    转自:http://hi.baidu.com/wjbzbmr/item/3f5e90f02b913ed66225d201

    View Code
     1 #include <cstring>
     2 #include <cstdlib>
     3 #include <cstdio>
     4 #include <algorithm>
     5 #include <iostream>
     6 
     7 #define N 50
     8 
     9 using namespace std;
    10 
    11 struct BIGN
    12 {
    13     int a[N];
    14 }dp[110],ans;
    15 
    16 int n;
    17 
    18 inline BIGN operator +(BIGN a,BIGN b)
    19 {
    20     BIGN c; memset(c.a,0,sizeof c.a);
    21     c.a[0]=max(a.a[0],b.a[0]);
    22     int jin=0;
    23     for(int i=1;i<=c.a[0];i++)
    24     {
    25         jin+=a.a[i]+b.a[i];
    26         c.a[i]=jin%10000;
    27         jin/=10000;
    28     }
    29     while(jin) c.a[++c.a[0]]=jin%10000,jin/=10000;
    30     return c;
    31 }
    32 
    33 inline BIGN operator *(BIGN a,int b)
    34 {
    35     BIGN c; memset(c.a,0,sizeof c.a);
    36     int jin=0;
    37     for(int i=1;i<=a.a[0];i++)
    38     {
    39         jin+=c.a[i]+a.a[i]*b;
    40         c.a[i]=jin%10000;
    41         jin/=10000;
    42     }
    43     c.a[0]=a.a[0];
    44     while(jin) c.a[++c.a[0]]=jin%10000,jin/=10000;
    45     return c;
    46 }
    47 
    48 void prev()
    49 {
    50     dp[0].a[0]=dp[0].a[1]=1;
    51     dp[1].a[1]=dp[1].a[0]=1;
    52     for(int i=2;i<=100;i++)
    53     {
    54         dp[i].a[1]=0;dp[i].a[0]=1;
    55         for(int j=1;j<=i;j++)
    56             dp[i]=dp[i]+dp[i-j]*j;
    57     }
    58     
    59 }
    60 
    61 void prt(const BIGN &a)
    62 {
    63     printf("%d",a.a[a.a[0]]);
    64     for(int i=a.a[0]-1;i>=1;i--) printf("%.4d",a.a[i]);
    65     puts("");
    66 }
    67 
    68 int main()
    69 {
    70     prev();
    71     scanf("%d",&n);
    72     ans.a[1]=0; ans.a[0]=1;
    73     for(int i=1;i<=n;i++)
    74         ans=ans+dp[n-i]*(i*i);
    75     prt(ans);
    76     return 0;
    77 }
    没有人能阻止我前进的步伐,除了我自己!
  • 相关阅读:
    VirtualBox 使用技巧
    ThreadPoolExecutor 线程池任务队列分析 与 利特尔法则(Little's law)
    AQS 与 LockSupport
    Matrix
    Fire Net
    Travelling
    Cannon
    N皇后问题
    Safecracker
    #include <algorithm>中sort的一般用法
  • 原文地址:https://www.cnblogs.com/proverbs/p/2743952.html
Copyright © 2020-2023  润新知