转载请注明出处,谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents by---cxlove
题目:给出长度为n的铁丝,将铁丝分为若干 部分,每个部分折成三角形,要求所有三角形相似。
三角形三边相等视为相等,三角形顺序不同视为不同
http://acm.hdu.edu.cn/showproblem.php?pid=4466
首先我们考虑周长为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)种情况,在一起的就组成大三角形
#include<iostream>
#include<cstdio>
#include<map>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
#include<set>
#include<string>
#include<queue>
#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 pb(a) push_back(a)
#define mp(a,b) make_pair(a,b)
#define mem(a,b) memset(a,b,sizeof(a))
#define LL long long
#define MOD 1000000007
#define lson step<<1
#define rson step<<1|1
#define sqr(a) ((a)*(a))
#define Key_value ch[ch[root][1]][0]
//#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
int dp[N],fac[N];
void Init()
{
dp[3]=1;
for(int i=4;i<N;i++)
{
dp[i]=dp[i-1]-ceil(i/3.0)+floor((i-1)/2.0)+1;
if(!(i&1)) dp[i]-=floor((i-(i-1)/2.0)/2.0);
if(dp[i]>=MOD) dp[i]-=MOD;
if(dp[i]<0) dp[i]+=MOD;
}
fac[1]=1;fac[2]=2;
for(int i=3;i<N;i++)
{
fac[i]=fac[i-1]*2;
if(fac[i]>=MOD) fac[i]-=MOD;
for(int j=2;i*j<N;j++)
{
dp[j*i]-=dp[i];
if(dp[j*i]<0) dp[j*i]+=MOD;
}
}
}
int main()
{
Init();
int cas=0,n;
while(cin>>n)
{
LL ans=0;
for(int i=1;i*i<=n;i++)
{
if(n%i) continue;
ans=(ans+(LL)dp[i]*fac[n/i])%MOD;
if(i*i!=n) ans=(ans+(LL)dp[n/i]*fac[i])%MOD;
}
cout<<"Case "<<++cas<<": "<<ans<<endl;
}
return 0;
}