首页 > Codeforces 900D Unusual Sequences:记忆化搜索

Codeforces 900D Unusual Sequences:记忆化搜索

题目链接:http://codeforces.com/problemset/problem/900/D

题意:

  给定x,y,问你有多少个数列a满足gcd(a[i]) = x 且 ∑(a[i]) = y。

 

题解:

  由于gcd(a[i]) = x,所以y一定是x的倍数,否则无解。

  那么原题就等价于:问你有多少个数列a满足gcd(a[i]) = 1 且 ∑(a[i]) = y/x。

 

  设f(k)为gcd(a[i]) = 1 且 ∑(a[i]) = k时的答案。

  只满足条件∑(a[i]) = k的数列共有2^(k-1)种(隔板法)

  然后要从中去掉gcd不为1的数列。

  每个和为k且gcd不为1的数列a1,对应着一个和为k的因数且gcd为1的数列a2。

  因为a1可以由a2整体放大而来。

  那么也就是f(k) = 2^(k-1) - ∑ f(p),其中p为k的因数(p != k)。

  搜索 + map记忆化即可。

 

  由于需要用到的f(k),k均为y/x的因数,最多sqrt(y/x)个。

  加上map的log复杂度,所以总复杂度为O(sqrt(n)*log(n))。

 

AC Code:

 1 #include 
 2 #include 
 3 #include <string.h>
 4 #include 
 5 #define MOD 1000000007
 6 
 7 using namespace std;
 8 
 9 int x,y;
10 map<int,int> mp;
11 
12 long long quick_pow(long long n,long long k)
13 {
14     long long ans=1;
15     while(k>0)
16     {
17         if(k&1) ans=(ans*n)%MOD;
18         n=n*n%MOD;
19         k>>=1;
20     }
21     return ans;
22 }
23 
24 long long dfs(int i)
25 {
26     if(i==1) return 1;
27     if(mp.count(i)) return mp[i];
28     long long ans=quick_pow(2,i-1);
29     for(int j=2;j*j<=i;j++)
30     {
31         if(i%j==0)
32         {
33             ans=((ans-dfs(j))%MOD+MOD)%MOD;
34             if(i/j!=j) ans=((ans-dfs(i/j))%MOD+MOD)%MOD;
35         }
36     }
37     ans=((ans-1)%MOD+MOD)%MOD;
38     return mp[i]=ans;
39 }
40 
41 int main()
42 {
43     cin>>x>>y;
44     cout<<(y%x==0 ? dfs(y/x) : 0)<<endl;
45 }

 

转载于:https://www.cnblogs.com/Leohh/p/8464069.html

更多相关:

  • Gym - 102082Ghttps://vjudge.net/problem/2198225/origin对于数列中任意一个数,要么从最左边到它不递减,要么从最右边到到它不递减,为了满足这个条件,就要移动,而移动的最少步数就是逆序对数。所以这个数要么往左移动,要么往右移动,所以两个取最小就好了 #include

  • 雪花算法根据时间戳生成有序的 64 bit 的 Long 类型的唯一 ID 各 bit 含义: 1 bit: 符号位,0 是正数 1 是负数, ID 为正数,所以恒取 041 bit: 时间差,我们可以选择一个参考点,用它来计算与当前时间的时间差 (毫秒数),41 bit 存储时间差,足够使用 69 年10 bit: 机器码,能编...

  •   这是一篇很水的blog 扫雷 link 一道很水的dp,考虑上一上,这一行,与下一行是否有雷即可 #include #include #include #include using namespace std; inline long long read...

  • P2429 制杖题 题目描述 求不大于 m 的、 质因数集与给定质数集有交集的自然数之和。 输入输出格式 输入格式:第一行二个整数 n,m。 第二行 n 个整数,表示质数集内的元素 p[i]。 输出格式:一个整数,表示答案,对 376544743 取模。 输入输出样例 输入样例#1:2 15 3 5 输出样例#1:60 说...

  • 把字符串看作是26进制的数,从后往前翻译,那么就可以把两个串变成对应的26进制的数字,那么只要把两个数加起来除以二就得到中间的串对应的数了,同理再转化回来就行了。但是这样会有一个问题就是串的长度有2e5,26的2e5次方显然不可求,所以需要对每一位进行手动的加和运算。对于两个串,我们假设a串从后往前的每一位对应的数值为a0, a1,...

  • 中国剩余定理说白了就是小学时候的韩信点兵的完全版。给定一系列数,给定条件是一个数MOD这一些列数的结果,问你最后这个数最少为多少。 抽象出来就是N个同余方程,利用扩展GCD就可以求得这一结果,本题给定的数都是互质的,因此处理起来就简单了。 代码如下: #include #include #inc...