限于本人DP處于背包--水平，不組聲（就是對神奇的數(shù)位DP，找出j個1的有多少個，然后j的多少個次方就行）

1 #include<bits/stdc++.h> 2 #define N 100005 3 #define LL long long 4 #define inf 0x3f3f3f3f 5 using namespace std; 6 inline LL ra() 7 { 8 LL x=0,f=1; char ch=getchar(); 9 while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();} 10 while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();} 11 return x*f; 12 } 13 const LL mod=10000007LL; 14 int l,wei[100]; 15 LL c[100][100],n,ans; 16 void pre() 17 { 18 for (int i=0; i<=60; i++) c[i][0]=1; 19 for (int i=1; i<=60; i++) 20 for (int j=1; j<=i; j++) 21 c[i][j]=c[i-1][j-1]+c[i-1][j]; 22 } 23 LL ksm(LL a, LL b) 24 { 25 LL sum=1; a%=mod; 26 for(;b;b>>=1,a=a*a%mod) 27 if (b&1LL) sum=sum*a%mod; 28 return sum; 29 } 30 LL solve(int x) 31 { 32 LL sum=0; 33 for (int i=l; i>=1; i--) 34 { 35 if (wei[i]==1) 36 sum+=c[i-1][x--];//計算的時候默認(rèn)i位是0，計算完后默認(rèn)i位是1，，應(yīng)該是這樣子的，，不會dp不會dp，，， 37 //這樣求出來的應(yīng)該就是最小值了233 38 if (x<0) break; 39 } 40 return sum; 41 } 42 int main() 43 { 44 pre(); n=ra()+1; 45 while (n) 46 { 47 wei[++l]=n&1; 48 n>>=1; 49 } 50 ans=1LL; 51 for (int i=1; i<=l; i++) 52 ans=ans*ksm(i,solve(i))%mod; 53 printf("%lld\n",ans); 54 return 0; 55 }

?

轉(zhuǎn)載于:https://www.cnblogs.com/ccd2333/p/6435349.html

總結(jié)

以上是生活随笔為你收集整理的bzoj 3209: 花神的数论题的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。