P2567 [SCOI2010]幸運(yùn)數(shù)字

題意：

我們規(guī)定只含6或8的數(shù)字為幸運(yùn)號(hào)碼，而幸運(yùn)號(hào)碼的倍數(shù)我們也認(rèn)為是幸運(yùn)號(hào)碼，問[l,r]中有多少個(gè)幸運(yùn)號(hào)碼？

題解：

第一反應(yīng)是數(shù)位dp，但其實(shí)不是，我們可以打表觀察下，幸運(yùn)數(shù)字其實(shí)很少，也就1000多個(gè)，所有我門全可以用容斥來去做。
先去搞定所有范圍內(nèi)的合法數(shù)字，然后直接容斥
對(duì)于有個(gè)幸運(yùn)數(shù)字X，在[L,R]中X的倍數(shù)的個(gè)數(shù)，就是R/X-(L-1)/X,但是兩個(gè)幸運(yùn)數(shù)字有可能會(huì)存在交集，即重復(fù)計(jì)算，所有要用到容斥，選一個(gè)幸運(yùn)數(shù)字-選兩個(gè)幸運(yùn)數(shù)字的lcm+選3個(gè)幸運(yùn)數(shù)字的lcm-…(容斥原理)
但是直接這樣做也會(huì)超時(shí)，需要剪枝：

對(duì)于兩個(gè)幸運(yùn)數(shù)字a，b，如果b%a==0，那么如果b/x就一定包含了a/x，因此這樣的b沒有必要，也就是幸運(yùn)數(shù)字的倍數(shù)我們就不用管了，這樣就剩下943個(gè)幸運(yùn)號(hào)碼

如果當(dāng)然的lcm大于r，就沒必要繼續(xù)搜索了，因?yàn)橐呀?jīng)超出范圍

我們將幸運(yùn)數(shù)字從大到小排序，可以使得lcm更快的超越上屆r

代碼：

#include<bits/stdc++.h> #define debug(a,b) printf("%s = %d\n",a,b); using namespace std; typedef long long ll; typedef pair<int, int> PII; clock_t startTime, endTime; //Fe~Jozky const ll INF_ll=1e18; const int INF_int=0x3f3f3f3f; inline ll read(){ll s=0,w=1ll;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')w=-1ll;ch=getchar();}while(ch>='0'&&ch<='9') s=s*10ll+((ch-'0')*1ll),ch=getchar();//s=(s<<3)+(s<<1)+(ch^48);return s*w; } void rd_test(){#ifdef ONLINE_JUDGE#elsestartTime = clock(); //計(jì)時(shí)開始freopen("2567.in","r",stdin);#endif } void Time_test(){#ifdef ONLINE_JUDGE#elseendTime = clock(); //計(jì)時(shí)結(jié)束printf("\n運(yùn)行時(shí)間為:%lfs\n",(double)(endTime - startTime) / CLOCKS_PER_SEC);#endif } const int N=20005; const int dig1=6; const int dig2=8; ll l,r,a[N],ans,tot; inline bool cmp(ll x,ll y){return x>y;} void init(ll x){if(x>r)return;if(x>0)a[++tot]=x;init(x*10+dig1);init(x*10+dig2); } ll gcd(ll x,ll y){return y?gcd(y,x%y):x; } void dfs(int now,int sum,ll k){if(now>tot){if(!sum)return;if(sum&1)ans+=(r/k-(l-1)/k);else ans-=(r/k-(l-1)/k);return;}dfs(now+1,sum,k);if((double long)k/gcd(k,a[now])<=(double long)r/a[now])dfs(now+1,sum+1,k*a[now]/gcd(k,a[now])); } int main(){scanf("%lld%lld",&l,&r);init(0);sort(a+1,a+tot+1,cmp);int tmp=0;for(int i=1;i<=tot;i++){for(int j=i+1;j<=tot&&a[i];j++){if(!a[j])continue;if(a[i]%a[j]==0)a[i]=0;}if(a[i])a[++tmp]=a[i];}tot=tmp;dfs(1,0,1);printf("%lld\n",ans);return 0; }

總結(jié)

以上是生活随笔為你收集整理的P2567 [SCOI2010]幸运数字的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

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