hdu 4466 Triangle

题目链接：点击打开链接

题目大意：将长度为n的截为1.全部全等的三角形2.都是完整的。问不同的方法有多少种

题目分析：参考大神的方法：推荐->点击打开链接

                 

首先我们考虑周长为m的三角形有多少个

我们假设 a<=b<=c

这里分为两种情况 b=c

那么对于c的上界，c最大取floor((x-1)/2)，因为a最少为1

对于c的下界，c最小取ceil(x/3)

则所有情况为floor((x-1)/2)-ceil(x/3)+1

第二种情况，b!=c

我们考虑(a,b,c)，由于a+b>c  ，b<=c-1所以(a,b,c-1)必定也为三角形

另外需要除掉a+b==c+1的情况 ，其中a+b+c=x，可以发现x必定为奇数，数目为floor((x-(x-1)/2)/2)

则得到代码中的递推式

我们还需要把三角形全部化为连长互质的。这里通过线性筛选就可以完成

解释一下这样做的目的。

接下来我们把n分为若干 个部分，每一个部分等长，而且都能组成若干个连长互质的三角形。由于题目要求互质。

我们把其中几个组成在一起，就是和本质三角形相似的大三角形。

假设有i个本质三角形，题目要求顺序不同视为不同，插板即可解决，i-1个空，每个空都可以选择隔开或者不隔开

2^(i-1)种情况，在一起的就组成大三角形

题目注意：dp[i], 表示周长为i的有多少个三角形

#include<iostream>
#include<cstdio>
#include<map>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<string>
#define inf 100000005
#define M 20005
#define N 5000005
#define maxn 300005
#define eps 1e-10
#define zero(a) fabs(a)<eps
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define __int64 long long
#define MOD 1000000007

using namespace std;
int dp[N],permi[N];
void Init()
{
    int i,j;
    memset(dp,0,sizeof(dp));
    memset(permi,0,sizeof(permi));
    dp[3]=1;
    for(i=4;i<N;i++)
    {
        dp[i]=dp[i-1]+floor((i-1)/2.0)-ceil(i/3.0)+1;//1.b!=c的情况和b==c的情况
        if(!(1&i)) dp[i]-=floor((i-(i-1)/2.0)/2.0);//dp【i】为偶数，则dp【i-1】为奇数，那么减去上面提到的a+b==c+1，的情况
        if(dp[i]>=MOD) dp[i]-=MOD;
        if(dp[i<0])dp[i]+=MOD;//一定要注意小于0的地方
    }
    permi[1]=1,permi[2]=2;
    for(i=3;i<N;i++)
    {
        permi[i]=permi[i-1]*2;
        if(permi[i]>=MOD) permi[i]-=MOD;
        for(j=2;j*i<N;j++)
        {
           dp[j*i]-=dp[i];
           if(dp[j*i]<0) dp[j*i]+=MOD;
        }
    }
}
int main()
{
    Init();
    int t=1,n; __int64 ans;
    while(~scanf("%d",&n))
    {
        ans=0;
        for(int i=1;i*i<=n;i++)
        {
            if(n%i) continue;
            ans=(ans+(__int64)dp[i]*permi[n/i])%MOD;
            if(i*i!=n) ans=(ans+(__int64)dp[n/i]*permi[i])%MOD;
        }
        printf("Case %d: %lld\n",t++,ans);
    }
    return 0;
}