题目链接：
http://acm.hdu.edu.cn/showproblem.php?pid=4466
题目意思：
一根长为N的木棒，长度分成若干个三角形，使得任意两个三角形都相似。对应顺序三角形全部全等的为同一种分法，求总共有多少种分法。
解题思路：
数学题。
先考虑分成一个三角形的情况。
不妨设a<=b<=c;
1、当b=c时，a至少为1，所以c<=(n-1)/2
而a<=b 所以n-2*c<=c =>c>=n/3; 故共有(n-1)/2-(n/3)+(n/3?0:1)种。
2、当b<c时，肯定有b<=c-1
此时若a+b>c 则a+b>c-1
如果a,b,c能构成三角形，则a,b,c-1也一定能够构成三角形。
反过来，如果a,b,c-1能够构成三角形，也即a+b>c-1 当a+b！=c时，一定能使a,b,c构成三角形。故可以通过dp[n-1]递推过来，然后减去a+b=c+1的情况。
此时n=2*c+1,c=(n-1)/2 只有当n为奇数的时候才有可能，而a+b=(n-(n-1)/2)，a<=b 所以这种情况共有 （n-(n-1)/2)/2种，化简得（n+1)/4
在考虑有多个三角形的情况。
假设N=a*b 则把b作为一个基本三角形，a个1，作为a个部分，中间有a-1个隔板，每个隔板可选可不选，一共有2^（a-1）种情况。
当所有的隔板都不选的话，就是一个大三角形，所以此时要把，刚才求得的一个三角形中三边不互质的数量减掉。而每个约数，对应的该约数能拆分成一个互质三角形的种数的不互质的三角形，乘以倍数就是大的不互质的三角形。
详见代码：
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<string>
#include<cstring>
#include<algorithm>
#include<vector>
#include<map>
#include<set>
#include<stack>
#include<list>
#include<queue>
#define eps 1e-6
#define INF 0x1f1f1f1f
#define PI acos(-1.0)
#define ll __int64
#define lson l,m,(rt<<1)
#define rson m+1,r,(rt<<1)|1
//#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
#define M 1000000007
#define N (5000000+10) /*
freopen("data.in","r",stdin);
freopen("data.out","w",stdout);
*/
int dp[N],bin[N]; //dp[i]表示一个三角形的情况，并且是互质的情况 void init()
{
dp[3]=1;
for(int i=4;i<N;i++)
{
dp[i]=dp[i-1]+(i-1)/2-i/3+(i%3?0:1); //b==c的情况
if(!(i&1)) //当前一个为奇数
dp[i]-=i/4; //减去，递推中不满足要求的部分
dp[i]=dp[i]%M;
if(dp[i]<0)
dp[i]+=M;
}
bin[1]=1;
for(int i=2;i<N;i++)
{
bin[i]=(bin[i-1]<<1)%M;
for(int j=2;i*j<N;j++)
{
dp[i*j]-=dp[i]; //减去一个三角形中三边不互质的部分,,
if(dp[i*j]<0) //也就是对应的约数的dp[i],对于每个互质的小的a,b,c 对应一个j*a,j*b,j*c
dp[i*j]+=M;
}
}
return ;
} int main()
{
int n,ca=0; init();
//printf("%I64d\n",dp[3]); while(scanf("%d",&n)!=EOF)
{
int ans=0;
for(int i=1;i*i<=n;i++)
{
if(n%i)
continue;
ans=(ans+1LL*dp[i]*bin[n/i])%M; //中间乘法可能会溢int
if(i*i!=n)
ans=(ans+1LL*dp[n/i]*bin[i])%M;
}
printf("Case %d: %I64d\n",++ca,ans);
}
return 0;
}