当前位置: 代码迷 >> 综合 >> HDU 4651 Partition (水+五边形定理)
  详细解决方案

HDU 4651 Partition (水+五边形定理)

热度:78   发布时间:2023-11-15 12:55:15.0

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4651

#include<bits/stdc++.h>
using namespace std;#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<int,int>
#define fi first
#define se second
#define mk(x,y) make_pair(x,y)
const int mod=1e9+7;
const int maxn=1e5+4;
const int maxm=1e5+100;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题目意思很简单,用不大于n的数组成n的方案数。
刚学完生成函数这道题应该不难,
利用五边形定理及其递推式,即
p(n) = sigma p(n-qi),其中qi是第i个五边形数,公式是
(3*x*x-x)/2。
详细证明我也不大懂。。QAQ。
*/
ll p[maxn];
int f1[300],f2[300];
int main(){for(int i=1;;i++){f1[i]=(3*i*i-i)>>1,f2[i]=(3*i*i+i)>>1;if(f1[i]>=maxn) break;if(f2[i]>=maxn) break;}mst(p,0),p[0]=1;rep(i,1,maxn) {for(int j=1,idx=1;;j++,idx*=-1){if(f1[j]<=i) p[i]+=idx*p[i-f1[j]];else break;if(f2[j]<=i) p[i]+=idx*p[i-f2[j]];else break;}p[i]%=mod,p[i]=(p[i]+mod)%mod;}int t;scanf("%d",&t);while(t--){int n;scanf("%d",&n);printf("%lld\n",p[n]);}return 0;
}

 

  相关解决方案