f(1)=1? f(2x)=f(x)? f(2x+1)=f(x)+f(x+1)?給出n<=10^6，求:所有的i滿足f(i)=n

第一道類歐算法

我們考慮一個(gè)性質(zhì) f(2x+1)=f(x)+f(x+1)=f(2x)+f(2x+2)

所以，顯然有f(2x+1)>f(2x)?f(2x+1)>f(2x+2)

那么現(xiàn)在我們知道了f(2x+1),自然考慮枚舉一個(gè)f(2x)

可以按照以下形式轉(zhuǎn)移:f(2x+1),f(2x)->f(x+1),f(x) （f(x+1)=f(2x+1)-f(2x)）

發(fā)現(xiàn)轉(zhuǎn)以后的兩者也是相鄰的，但是不知道奇偶性，較大的那個(gè)肯定是基數(shù)

那么就可以進(jìn)行輾轉(zhuǎn)相減，邊界就是f(x+1)=f(x)=1

注意這里可以借鑒輾轉(zhuǎn)相除的方法，不然會(huì)Tle

#pragma GCC optimize("O3") #pragma G++ optimize("O3") #include<vector> #include<stdio.h> #include<string.h> #include<algorithm> #define LL long long #define M 998244353 using namespace std; vector<LL> A; int p[1000010]; LL gcd(LL x,LL y,LL& a,LL& b){if(x==y) return x>1?0:b=((a&=0)+=1)++;if(x>y){if(!gcd(--x%y+1,y,a,b)) return 0;b=b*p[x/y]%M; a=(b+1)%M; return 1;} else {if(!gcd(x,--y%x+1,a,b)) return 0;a=a*p[y/x]%M; b=(a-1+M)%M; return 1;} } int main(){freopen("func.in","r",stdin); freopen("func.out","w",stdout);LL n; scanf("%lld",&n); p[0]=1;for(int i=1;i<=n;++i) p[i]=(p[i-1]<<1)%M;for(LL x,y,i=1;i<=n;++i){if(gcd(n,i,x,y)) A.push_back(x);if(x<0) printf("%d %d\n",x,i);}sort(A.begin(),A.end());for(LL i=0;i<A.size();++i) printf("%lld\n",A[i]); }

轉(zhuǎn)載于:https://www.cnblogs.com/Extended-Ash/p/9477142.html

總結(jié)

以上是生活随笔為你收集整理的Jzoj5317 Func的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

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